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The atmosphe rically formed protective 


atmosphere from surfaces of an aluminum alloy, 


carbon and stainless steels 


fil 


ms were removed by 


two brasses, 


abrasion in 


& magnesium alloy, 


an argon 
and low- 


The resulting surfaces were 0.12 to 0.76 volt more negative with 


respect to a calomel electrode than surfaces prepared and measured under normal atmospheric 


conditions 


Appreciable changes in electrochemical solution potentials of notched specimens, 


stressed in tension, occurred at or just above stresses at which the true stress true strain 


curves deviated from the modulus lines 


These changes in potential were caused by rupturing 


of the protective films at the roots of the notches and were of the order of 0.16 to 0.70 volt 


at failure, depending on the material 


Stress corrosion is postulated to occur in corrosive 


media, at stresses sufficient to rupture the protective film, by electrolytic action between 


the filmed (cathodic) and film-free (anodic 


1. Introduction 


ress-COrrosion cracking has been defined as the 
spontaneous failure of metal by cracking under the 
ombined action of corrosion and stress, residual or 
ipplied” [1]. Service failures have frequently re- 
sulted from residual stresses set up in the material or 
structure during fabrication. These may 
iit from the spinning or deep drawing of sheet 
rass, the work-hardening of aluminum-alloy tubing, 
the pressing of oversized bushings into aluminum- 
alloy fittings, the riveting or welding of steel struc- 
tures, Or even as the result of the application of 
ssive torques in the assembling of brass-pipe 
ructures [2]. Stress-corrosion cracking may be 
ther intererystalline, as in aluminum alloys, alpha 
rass, and low-alloy steels, or transcrystalline, as is 
sually the case in magnesium-base alloys and the 
istenitic stainless steels, 

It is agreed by most workers in the field that some 
degree of stress-corrosion cracking can be produced 
1 many alloys and commercially pure metals by the 
proper combinations of stresses and corrosive media. 
lhe ASTM-AIMME Symposium on Stress-Corro- 
sion Cracking of Metals [3], published in 1944, gives 

excellent summary of work in the field to that 
lat. A number of papers have also appeared in 
ecent vears [4 to 8, 10] that contain excellent 
ibliographies. For the most part, the data in the 
terature are from investigations of the relative 
susceptibilities of various materials to stress-cor- 
rosion cracking, susceptibilities of particular materi- 
\o stress corrosion in specific corrosive media, or 
susceptibilities of particular materials to 
osion after various mechanical or thermal 
treatments. 

lhe purpose of this investigation was to study the 
itechanism of stress corrosion and not to determine 
the relative susceptibilities of various materials to 

ss corrosion. Recently the various theories of 
inechanism of stress corrosion were summarized 
Harwood [10]. Of these theories the “General- 
theory of stress corrosion”? proposed by Mears, 
wo, and Dix [11], the “Mechanical theory” of 


stresses 


stress 


col 


in brackets indicate the literature references at the end of this paper 


areas, 


Keating [5], and the “Film rupture theory’’, or better, 
“Film rupture mechanism,” are pertinent to the 
present paper. It is considered in this paper that 
the “Film rupture mechanism” and the “Mechanical 
theory” are logical parts of the Generalized theory. 

In the Generalized theory, corrosion is postulated 
to occur by an electrochemical mechanism along 
localized paths that are anodic to the surrounding 
metal. There will be stress concentrations at the 
bases of localized corroded areas; the deeper the 
attack and “the smaller the radius at the base 
of the path the greater would be the stress concentra- 
tion. Such a condition would act to pull the metal 
apart along more or less continuous paths. At suf- 
ficient concentrations of stress, the metal might start 
to tear apart by mechanical action the tearing 
action described above would expose fresh metal, 
unprotected by films, to the action of the corrosive 
environment. Because this freshly exposed metal 
is more anodic, an increase in current from the base 
of the localized path to the unaffected surface would 
be expected "It is seen from this quotation 
that Mears and his coworkers [11] postulate a mech- 
anism in which protective films on metal surfaces are 
ruptured. Keating [5] suggested that each advance- 
ing section of the crack in time encounters an obstacle 
in the form of a nonmetallic inclusion, lattice dis- 
continuity, or unfavorably orientated grain boundary. 
In each case, the result may be reduction in the stress- 
concentration effect to a level at which cracking 
cannot proceed. Such a would explain 
branching of stress-corrosion cracks, sudden changes 
in direction at grain boundaries, ete. Gilbert and 
Haden [12] report the results of an investigation on 
an aluminum alloy containing 7 percent of mag- 
nesium in which they found that the cracking process 
was not continuous but stepwise 

There are few published data from which the elec- 
trochemical solution potential of the freshly torn 
metal at the root of-the crack can be estimated. 
Work at the National Bureau of Standards as early 
as 1940 showed that the solution potential of freshly 
quenched or freshly abraded 245-T aluminum-alloy 
sheet was much more electronegative (anodic) than 
that of the same material permitted to stand for 


process 





some minutes before measurements were made [13]. | cartridge brass, 61-percent-Cu_36-percent-Zn-3-) 
Similar results were recently obtained at the Bureau | cent-Pb leaded brass, AZ31f magnesium alloy, |, 
on 18-percent-Cr-8-percent-Ni austenitic stainless | carbon steel, and type 302 stainless steel. 
steels. Mears and Brown [14] report a solution po- | chemical compositions are given in table 1. 
tential greater than —3.0 v for commercially pure | aluminum alloy, the leaded brass, and the st 
aluminum from which the film was removed by | were obtained commercially. The magnesium a 
scratching the surface with a glass stylus in the | was supplied by the Dow Chemical Co., and 
electrolyte. Druet and Jacquet [15] report values | cartridge brass was made in the Bureau’s Exp 
of —0.75 v for mechanically polished specimens con- | mental Foundry from electrolytic copper and 99 )9- 
taining 99.99 percent of aluminum (presumably with | percent zine. Prior to machining of the specimens 
a protective film intact) and —1.384 v for the same | the materials were heat treated as follows: (a) The 
material with measurements made immediately after | aluminum-alloy stock was solution heat-treated at 
electrolytic polishing. | 920°+10° F, and quenched in water at approxi- 
The author is not aware of any data giving the | mately 70° F, (b) the cartridge brass was rolled to 
magnitude of changes in solution potential produced | diameter of about 0.62 in. and then annealed for | 
by rupture of protective films on the metals or alloys | hour at 1,560° F, after which the Brinell hardness 
that are particularly susceptible to stress corrosion. | number was 66 and the grain size approximately 
Such data are necessary to establish the validity of | 0.025 mm, (c) the low-carbon steel, originally in the 
the so-called film rupture process in the generalized | cold-rolled condition, was heated to 1,625° F and 
theory of stress corrosion. In the early phases of | cooled with the furnace, (d) the stainless steel was 
the investigation of the mechanism of stress corrosion, | heated to 1,950° F and quenched in cold water. The 
the following have been determined and are reported | magnesium alloy, in the form of an extruded rod 
in this paper: (a) the electrochemical solution poten- | in. in diameter, and the leaded-brass rod, \ in 
tials of several alloys free from the protective effects | in diameter, were used as received. The magnesium 
of films,’ and (b) relationships between the applied | alloy evidently had been extruded at a sufficiently 
stresses and the electrochemical solution potentials | high temperature that the metal had recrystallized 
of these alloys | following extrusion. The leaded brass had a Brinel! 
hardness number of 113 and a grain size of 0.015 to 
2. Materials and Methods of Test 0.025 mm. The tensile properties, as determined on 
subsize (0.355-in. diameter) ASTM round tensile 
specimens, are given in table 2. 
When the protective film on 248-T4 aluminum 
alloy is broken in air, it tends to repair itself quite 


It is recognized that contact of an electrolyte, at least with the aluminum and | f ‘ a be al f : 2 : : . 
magnesium alloys studied, may produce a protective film on the metal surfaces rapidly ‘ Phere fore, specimens were prepared by 


The materials used in this investigation were 245- 
T4 aluminum alloy, 70-percent-Cu-30-percent-Zn 


Tasie ! Chemical composition of metals used 


! | } 
Material Pb Mn Me 


48 aluminum alloy * 
Cartridge brass 

Leaded brass 

AZ31 magnesium alloy * 
Low-carbon stee 

Ty pe 302 stainless steel 


* Determined spectroscopically. © By difference Manufacturer's analysis 


Taste 2, Mechanical properties of metals used | placed on the prepared surface and brought into 
electric contact with a calomel electrode of the 
Ultimate viet | Rtenention saturated KCl] type by means of an agar-agar bridge. 

Material tensile . . 
strenath strength =| in 14 in Potential measurements were made with a pH meter, 
or for values in excess of 1.3 v, with the pH meter in 

1,000 thjin.? | 1,000 thin? | Percent | series with a L&N type K potentiometer. 

pa ni a3 a It was recognized that the preparation of the speci- 
Leaded brass =s —$ » men surfaces by abrasion left the surfaces to be 

AZA1f magnesium alloy 11.0 s) il : . ene 
Low-carbon stee! 4.7 1.5 “0 measured in the cold-worked condition and might 
Type 302 stainless steel m1 3. " produce phase changes in the stainless steel. In 
some earlier work at the Bureau the author reduced 
abrading them with “400 aloxite’’ metallographic | the thickness of some 245-T4 aluminum-alloy sheet 
polishing paper, and all measurements were made in approximately 5 percent without appreciable chang- 
an argon atmosphere* in a dry box. In making | Ing Its solution potential. Druet and Jacquet (15) 
measurements, drops of the desired electrolyte were showed that severely cold-worked high-purity alumi- 
, num was only about 0.04 v more negative than the 


: ' sis of the argon by the mass spectrograph method showed 0.2 percent | . hy - 
(hata . annealed material. The author has seen no data on 
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fect of cold-work on the solution potentials of 

materials studied nor no the effect of a phase 
ve on the solution potential of the stainless steel. 
ress-electrochemical solution potential relation- 
; were obtained on notched specimens made in 


rdanee with the drawing in figure 1. The 
imens were covered, except for the threaded 
:, with two coats of Glyptal red enamel. After 


enamel had dried it was cut through by means of 
zor blade so as to expose an extremely narrow 
th of metal around the specimen at the root of 
noteh The specimens were placed in cells in a 
sile-testing machine, figure 2, with a calomel 
trode placed directly in the cell if potassium 
ride was used as the electrolyte or in a beaker 
nected to the cell by an agar-agar bridge if some 
her electrolyte was used. The specimens were 
ssed to failure in the corroding solution, and the 
trochemical solution potentials that corresponded 
given loads were determined either with a pH 
meter used as a vacuum-tube voltmeter or a Miller 
nultimeter (used asa high-resistance potentiometer). 
The free cross-head speed of the testing machine was 











held between 0.05 and 0.10 in./min. In order to 
obtain true stress-true strain and load—true stress 
be . 2 144" ——+}-———- 2 1/4" — oat 
R=008 —| THREAD 1/2-13 
Vv 
(00 (it 
—_— 
| 90° L 
i. : 5/8" 
“~ 0.20020! 
be —_— — 4 1/2"— — 000 
Ficure 1 Details of notched specimens, 





Figure 2, in place in cell 


Notched specimen (A) 
(B) in testing machine. 


trode (C) was placed directly in cell if potassium chloride was the corrodent, 
se an agar-agar bridge was used. 


relationships, three or more additional notched 
specimens of each material were broken in the 
uncorroded condition.” True stress-electrochemical 
solution potential relationships were then obtained 
by combining load-true stress and load-electro- 
chemical solution potential data. 


3. Results and Discussion 
3.1. Electrochemical Solution Potentials of the Metals 


The different electrolytes used and the electro- 
chemical solution potentials of the various materials, 
determined in air and in argon, are given in table 3. 
All values are in terms of a calomel electrode of the 
saturated KCl type and were determined at approx- 
imately 25° C. The values reported for atmospher- 
ically formed filmed surfaces were equilibrium values 
obtained by immersion of the prepared surfaces in 
their specific electrolytes; the values reported for the 
specimens prepared in argon were instantaneous 
maxima and not necessarily equilibrium values. 
These values are given inasmuch as they represent 
the conditions under which changes in electro- 
chemical solution potential values reported in a 
later section of this paper were obtained. 


TABLE 3 Electrochemical solution potentials (with respect to 
calomel electrode of saturated KCl with 


surface S prepared and measurements made im air or argon 


type) of materials 


Electrochemical 
solution 
potentials in 


Material Electrolyte 

Air Argon 

Volt Volt 
248-14 aluminum alloy KC] (saturated 1.67 1.43 
7, -Cu-30"-Zo brass 5°) NH.OH4 47 0.06 

5°) (NH),CO 
61% -Cu-36%-Zn-3%-Pb do 39 Us 
brass 

AZ31 magnesium alloy 3.5% NaCl+2% KeCr0, 1.51 1. 45 
Low-carbon steel 5% NHiNO “”) 0.72 
Type 302 stainless steel KC! (saturated a» 78 


The data indicate that, under normal atmospheric 
conditions, protective films are formed on the surfaces 
of these alloys. These filmed materials are in some 
cases as much as 0.76 v more electropositive than the 
same materials with surfaces prepared under argon. 
The largest differences in potential are found on the 
aluminum alloy, the brasses, and the stainless steel, 
all of which are relatively corrosion-resistant mate- 
rials. Those materials that corrode more readily, 
such as low-carbon steel and the magnesium alloy, 
had very much smaller differences in potential 
between filmed and unfilmed surfaces. 

In order to determine whether the difference in 
solution potential of the cartridge brass was due to 
the formation of an oxide film or to some other cause, 
the following investigation was carried out. <A 

‘Lhe use of notched specimens, of course, introduces stress concentrations of 
undetermined magnitude at the root of the notch. The maximum stress at the 
root of the notch is a function not only of the geometry of the specimen, but also 
f the applied stress. Although the radii of the notches used were small (approxi 


mately 0.005 in.), they were nevertheless large compared to those at the roots of 
atress-corrosion cracks in many materials 
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specimen was placed im the dry box, and a surface 
was prepared and measured under argon. Ap- 
proximately 20 percent of the argon was then 
replaced by dry oxygen, and additional measure- 
ments were made on the same surface. The experi- 
ment was repeated with wet oxygen, wet carbon 


dioxide, and dry ammonia gas. The results are 
given in table 4, 
Tanue 4 Electrochemical solution pole ntials of 70%-Cu 30% 


Zn hrass with surfaces prepared man argon atmosphere and 


measurements made ina mired almos phe re 


Electrolyte was NH,OU4 (NH pCO 

Electrochemical 

Atmosphere olutio 

potential 
Vad 
Argon mi 

~y argon + ar iry wt “) 
“yr weon? ar wet © st 4) 
; =r argon + ar wet CO ht x 
wy. argon + aF lry Nu “it “l 


| : 


It will be noted that in every 
potential was more positive after the addition of 
the second gas, the change being most marked in the 
case of oxygen, either dry or wet. In the atmos- 
phere containing ammonia, the surface of the alloy 
adjacent to drops of liquid used in making the meas 
urements became tarnished; its final potential was 

0.31 v compared to —0.47 v for brass measured in 
air. Measurements made on the surfaces of brass 
specunens covered with a blue corrosion product 
after exposure to an atmosphere containing ammonia, 
carbon dioxide, air, and water vapor were 0.75 ¥ more 
electropositive than film-free surfaces of the same 
alloy. 


case the solution 


3.2. Stress-Electrochemical Solution Potential 
Relationships 


The relationships between the true stresses and 
electrochemical solution potentials of the materials 
are shown in figures 3 to 8. The maximum values 
of the potentials were usually obtained just at the 
failure of the specimen as the result of the exposure 
of large areas of unfilmed material. These values 
are of course not shown on the figures but are given 
in table 5. In any measurement of electrochemical 
solution potentials the measured potentials are 


Vaxrim um change in electrochemical solution pote n- 
tial with stress 


Tasie 5, 


Change in 
electrochem 
ical solution 

potential 


Material Electrolyte 


Volt 
MS-T4 aluminum alloy Saturated KCI 0. 610 
7,-Cu- we) -Zn Brass Y) NHOH+5% (NHoO CO 700 
61% -Cu-%-Za-3%-Pb do So 

brass 

AZ31 magnesium alloy Lo, NaCl+ 2, Ker, fH) 
Low-caurbon stee “| NENO Th 
Ww? stainleas steel taturated KCI WS 
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tionships for 24S-T4 aluminum alloy specimens. 


Corroding medium, saturated KCl. Data are from seven specimens 
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Corroding medium, 5 percent NH,OH +5 percent (NH) »COs.. 
seven specimens. 
































ssarily the integrated potentials over the entire 
in contact with the electrolyte It is therefore 
« expected that potential measurements made 
notched specimens of the same material may 
somewhat from specimen to specimen, hence 
ranges of potential at given true stresses rather 
average curves are shown in the graphs 
There was an appreciable change in the solution 
ential of the aluminum alloy (fig. 3) at stresses at 
ch the true stress-true strain curve deviated 
m the modulus line, that is, at stresses at which 
tic deformation had occurred. Marked changes 
he potentials of the brass specimens (figs. 4 and 5) 
irred at stresses that were 10,000 to 25,000 Ib/in 
ove those at which the true stress-true strain 
ves deviated from the modulus line. The solution 
ential of the magnesium alloy changed appreciable 
tresses only about one-third of those necessary to 
niuce deviation from the modulus line (fig. 6). The 
eld point in the low-carbon steel was accompanied 
a marked change in potential (fig. 7). 
iflicient to produce a deviation of the true stress 
1¢ strain curve from the modulus line for stainless 
eel also produced appreciable changes in the 
olution potential (fig. 8). The electrochemical 
olution potentials (except for the low-carbon steel) 
became increasingly more negative as the stresses 
were increased. For every material (except low- 
arbon steel) stresses were reached, as is seen in 
livures 3, 4, 5, 6 and 8, at which the slope of the true 
stress-potential curves decreased markedly, and the 


Stresses 
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tionsh ps for 1Z31 magnesium alloy specimens 
Corroding medium, 3 percent NaCl +2.0 percent KoCrOy Data are from 


eight specimen 


potential increased rapidly to the failure of the 
specimen 

It is postulated that during the loading of the speci- 
mens in tension the protective film is broken over 
minute areas of the surface. If the potential of the 
unprotected area alone could be measured, it would 
presumably be of the same order of magnitude as 
that determined on a film-free surface. The area 
where the film has been ruptured is very small com- 
pared with that of the filmed surface of a specimen, 
and consequently the measured potential differs by 
only a small amount from that obtained on an un- 
broken filmed surface. As the stresses are increased 
the ratio of the film free to the filmed surface of a 
specimen increases, and the solution potential be- 
comes more negative, until at the instant of fracture 
the area of unfilmed surface may be very large com- 
pared to that of the filmed surface, and the potential 
may closely approach that of completely unfilmed 
material. 

At a constant stress the protective film tends to 
repair itself. This was demonstrated by results 
obtained on aluminum-alloy specimens. After an 
appreciable change had been obtained in the solution 
potential during loading, the applied load was held 
constant, and within | minute the solution potential 
was the same as that of the unstressed specimen. 
After the specimens broke (exposing large unfilmed 
areas) the Sastsocliouteal solution potentials rapidly 
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assumed the same values as those of the filmed and 
unstressed materials. 
Examination of 
notched specimens and the marked changes in the 
potentials of all of the materials except low-carbon 
steel at fracture indicated that the final failure had 


the fractured surfaces of the 


been almost instantaneous. In the low-carbon steel 
specimens, however, both the appearance of the frac- 
ture and the small change in potentials at failure 
indicated that prior to the final fracture, cracks had 
penetrated into the material from the roots of the 
notches and that the progression of fracture was not 
rapid and was probably discontinuous. Conse- 
quently, at failure the amount of freshly exposed 
unfilmed metal was too small to produce any marked 
change in the potential 

The data indicate that when the stresses are high 
enough to cause plastic deformation the atmospheri- 
cally formed protective film will be broken and the 
exposed metal will become electronegative with re- 
spect to the filmed material. These film-free areas 
will, in general, be very small compared to the filmed 
areas. If such areas are connected to filmed areas by 
a conducting liquid, for example, by a drop of liquid 
containing dissolved gases or salts, a short-circuited 
electrolytic cell will be set up, the circuit being com- 
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Corroding medium, 5-percent NH\NO Data are from nine specimens 
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tionships for type 302 stainless steel specimens. 


Corroding medium, saturated KCI Data are from six specimens 


pleted through the metal itself. As the film-free 
(anodic) area is generally very small compared with 
the filmed area, the current density over the anodic 
area will be high, and corrosion will proceed at a 
rapid rate antl @ readjustment of stresses permits 
se cet ws of the protective surface film. Corrosion 
will then become more general (or stress corrosion 
will proceed along different paths) until the stress 
concentration at the particular path under considera- 
tion again becomes sufficient to rupture the pro- 
tective film. The author believes that successive 
breaks in the time-potential curves reported by 
Gilbert and Haden [12] indicate that such processes 
were occurring, probably at several points, on their 
test specimens. 

The experiments reported in this paper confirm 
the general idea of the film-rupture mechanism 
suggested by Mears, Brown, and Dix. In the theory 
as first presented, it was postulated that tearing of 
the metal was necessary to expose film-free material 
Hlowever, the present data indicate that film rupture 
occurs with comparatively small amounts of plastiv 
deformation of the underlying material. When the 
film ruptures, stress corrosion may start in the film- 
free areas at stresses much below those normally 
required to tear the metal apart. Unless conditions 
permit film repair, stress corrosion may progress to 
cause complete failure in a relatively short time. If 
the film is restored and the stresses are not increased 
or readjusted, stress corrosion may cease. If, how 
ever, the stresses are readjusted by any means what- 
soever, for example, by a change of external load or 
by corrosion, the film may again be broken and 
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corrosion proceed discontinuously to failure. 
ast process may be that which occurs in the 
ry stalline type of stress corrosion. 


4. Summary 


a part of a detailed study of the mechanism of 
s corrosion, the electrochemical solution poten- 
of filmed and film-free surfaces and the true 
electrochemical solution potential relationships 
everal alloys were determined. 
he electrochemical solution potentials of surfaces, 
which the atmospherically formed films had 
removed by abrasion in an inert atmosphere, 
0.10 to 0.76 Vv more electronegative (with 
ect to a calomel electrode) than filmed surfaces 
he same materials 
\ppreciable changes in electrochemical solution 
entials accompanied plastic deformation, and just 
or to fracture the solution potentials were 0.16 
to 0.70 v (depending on the allov) more electro- 
ative than the unstressed material 
The results of the investigation indicate that at 
esses sufficiently high to produce failure of the 
protective film on metals, stress corrosion occurs in 
a corrosive medium by electrolytic action between 
film-free (anodic) and filmed (cathodic) areas. The 
process proceeds continuously or discontinously to 
failure, depending on whether or not the subsequent 
adjustment of stresses permits reforming of the 
protective film. 
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A Note on the Bounds of the Real Parts of the Characte:- 
istic Roots of a Matrix’ 


P. Stein * 


I'wo theorems are given for the bounds of the real parts of the characteristic roots of 


an n> 
is then specialized in several ways so as 
elements of the matrix 


1. Fundamental Theorems 


Let A be Let M and m be 


the maximum and minimum real parts of its char- 


(dy,) an nxn matrix 


acteristic roots 
Let u,, Ue, , Uy, be a set of positive numbers 
In this note we prove the following two 


The orem «4 1 


theorems 


. {1 - l - / 

i) M<smax;- (a,,+4,,) >‘ la a, u 
(2 au, St \ 
\1 - l - ) 

~ 

(ui) m=>mimn (a,,-+-d,,)- a a,, u 

r (2 su, — ie \ 
Theorem B. Tf the elements a,, are real and the 


_ rs are real and non-negative, then 


l —— 

): M smax (: +—S> a,,u, }. 
. “<= 
sr 


Using these theorems and giving specific values to 
the set u,, we obtain some further inequalities for M 


elements a 


' 1 
min («- +—S) a,,u, 
' Ure 


and m. These inequalities are numbered theorems 1 
to 4. Theorem B is substantially a theorem of L. 
Collatz? 

To prove theorem A(i), we write A= //-+-iA, where 
H and A are Hermitian matrices whose elements h,, 
and k,, are given by 

A.* Ay, Gor k, 7 


(@,,—a,,). (1) 


It is known that the real parts of the characteristic 
roots of A are bounded above and below by the 
maximum and minimum characteristic roots (real) 
of H. Hence to prove A(i) it is sufficient to prove 
that the maximum characteristic root of // is less 
than the right-hand side of (i), and to prove A(ii) it 


! This work was performed under National Bureau of Standards contract with 
the University of California at Los Angeles, and was sponsored in part by the 
Office of Naval Research 

* University of Natal, South Africa, 
Angeles 


and University of California at Los 


? Math, Zeit. 48, 221 (1942 





n matrix, depending on the use of an arbitrary set of n positive numbers 
to lead to theorems for the bounds in terms of the 





The set 


is sufficient to prove that the minimum characteri 
root is greater than the right-hand side of (ii). 

Let \ and uw be the maximum and minimum char- 
acteristic roots of 77. If x is a vector with n com- 
ponents, then A,, is one value of 2*//,/x*rx for all 
values of r, where z* is the transpose of the conjugate 
of x. Hence \2h,,, for all r, and uw sh,, for all r 

Let Ul’ be the diagonal matrix with elements «,, 


My, . . . , Up. Let B=U-'HU and let its elements 
be b,,. Then 
en ea ee a 2 
u, 


As B is a transform of H/, B and // have the same 
characteristic roots, so that ’ and uw are character- 
istic roots of B. If @ is any characteristic root of B, 


then a@ lies in at least one circle with center 6,, and 


radius 
n 
—" 
-_— b.. 
sel 
onwr 
Thus we have 
n" 
a—b,,|s >> |6,, 
sel 
avr 
for at least one value of r. 
As }\—6,, 20, and w—b,, S0, we have 
n 
a ‘ 
156,455 |b. 
sel 
sar 


for a least one value of r; and 


u2b,, b,, 


snr 


for at least one value of r. By (1) and (2), this com- 
pletes the proof of theorem A. 

Theorem B may be proved by a modification of the 
proof given by L. Collatz (see footnote 3), or alter- 
natively, as follows. 

Let C=A+NE, where E is the unit matrix, and 
N is a positive number so chosen that a,,4+-N>0 for 
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n- 


he 


nd 
or 








('lU’, where U is the diagonal matrix 


Let D=l 


h elements u;, Ue, « Alias 


he characteristic roots of C and PD are the same, 
ereas the characteristic roots of A are those of C 
ninished by NV 

if d.. is an element of DD). then d (1/u.)a..u 

ain the maximum positive root of a matrix with 
nnegative elements is bounded above and below 
the maximum and minimum values of the sums 
elements of a row. The maximum positive root 
such a matrix is also a root of maximum mod- 
is, and so greater than the real part of any other 
ot 

he sum of the rth row of the matrix J) is 


a 
ostr 


a 


Hence theorem B follows by subtracting N from 
ach row 


2. The Case of Real Elements 


In this section we shall suppose all the elements of 
i to be real Let 


R.=DS4n, Ci=dSe 
s= | s=| 


Theorem 1. (i) If a,.+a,2 0, r#¥s, if R.4+-C.>0. 
for all r, if ay2a,, for all r, and a,Za,,, r#i, of 
R,+0C,2 (R,+-C,) for all r and RAC R,+C,, 

p then 
R,+C R,+C, 

a + > ): 


VJ < max ( a.,4 


~ 


ii) Tf Qnm <a,, for allr and a,,Sa,,, r#m, then 


R,+C 


ae , R,+C, 
m min ( Gan” = ’ Gana 5 ) 


lo prove (i), We may suppose i= 1, and let 


Be tC, ' R, +c’ 


9 , 9 


We apply A (7), where we may leave out the mod- 
ulus sign. We take u,=o and 


rxél 


R,-+-¢ 
~— oa 
» 
then 
R,+-¢ R,+-C, , 
. ‘ 0 
») ) 
and 
max U,=o 
Hence in either case 
l aa) As) , « 
ays > : U.S 4,70 ° (3) 


Again for r#1, we have 


1 */a,,+4a l */a,,+a 
yy 5 aT < " - . ’ Xs re ar 
a,,4 ( = ju, a,,-+-max u,-—> ( = 
Urgel é ar reewl\ ~ 
avr avr 
=r, tC. (4) 


From (3) and (4), part (i) of theorem 1 follows. 
Part (ii) of theorem | follows similarly from A(i). 

Theorem 2 Let a,, be real, a,,20, r#s. (i) If 
ay2a,,, for all r, a,,2a,,, r#i, and R,=R, for all r, 
R,>R,, rAp, then M smax (a,,4+-Ry, ayy+R,). (Ud) 
If Omn =a, for all r, Ay» Sa,,, rA~Am, and Rt SR,, for 
al r, R'SR,, rp’, then M2min (a,,+, 
Gant Rt). 

‘The proof of this is similar to the proof of theorem 
1 using theorem B, and may be omitted. 

Theorem 3. If a,, is real and a,,20, r#s, then 


7 


— 
( iz cons j 
We apply theorem B and take 


n n 1 
M < max (< +)>3 a3, —>> a2,! . 


n 1 
> 2 
u a 
Ser 
where we suppose u,#0. By the Hélder-Schwartz 


inequality 


| n | n 1 n i 
s> ‘ Sia Vl w)\ 
- A,Us ; ( a... u; 
r sel le \ s=l s=1 
sx sr ser 
n n 1 
7 2 ‘ 2 2 
S u“— S a... . 
sel g= | 


ser 


The theorem follows from the definition of u,. 
Theorem 4. lf a,, is real and a,,20, r#¥s, and if 


1, 4a,,+a% R,+C r 
ay+—), = u,<a,,4 : z ‘max U,: ~- 
i 20 mi | then 
R, +0, M <max a,, 
” g, | r 

— max wU, e’: ' min}| ( Ss} 1, —max L: |*( > } max L,)*}. 

sl | ( .a@l r a= r ) 
- ; 107 
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BR We take u, 





L} and apply the Hélder-Schwartz | If z 
inequality and obtain 


at a 
= L(2 1-1) 


fr, Sth 
max ) LS L 


, y>0, r+y—a, we have 


a’ —y* >a* —(a—r)*=2ar—27’>ar 


Applying this to (5) we get 
vn 


L, ( >) L,—L, ); , <( > L, )- max L 


s=l 


From (6) and (7) and theorem B the 
follows. 


Los ANGELES, October 23, 1951 
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Chart for the TE ,, Mode Piston Attenuator 
Charles M. Allred 


A nomogram is given and described that expedites the determination of the dependence 
of attenuation on frequency, conductivity, and radius in a cylindrical waveguide, TE,, mode, 
piston attenuator. 




















a0 
110° 5 2 25 3 35 4 45 5 6 7 8 9 imoe 
A. 4. eee ee eee eee ee Peers eeeeee peer | 
' 
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gs 
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— 
ae 
— 
8 1 » 
30 9 . 
z 
, et 
os a A. yr-(2t) 2 
* i *15.99 22 f 
. ta (Ge) asa Tear) 
* + 3." 
wo e G, radius, meters 
aif O , conductivity, mhos per meter 
5 f , frequency , cycles per second 
a 2 : A, attenuation, db per rodius 
©, © 
at ew C, velocity of light, 3x10® meters per second 
a 
Ses wi” PJ, permeability of spoce , 417 x 10’ henry per meter 
vr @ x 
4 6 $ 6 7 » 9 ' i 12 13 \4 15 16 
2 — oe ee ee -” " wall a oe oe meee! — 
a SS 
6.00 99 98 .97 96 15.95 oe ee 15.90 89 688 87 86 15.85 84 83 82 Bi 5.80 79 78 77 6 S75 74 «73 72 a 1570 
A 
Figure 1. TE, mode piston attenuator chart. 
Piston attenuators (waveguide-below-cutoff) are The accompanying nomograph expedites the 


extensively used as adjustable attenuators because | determination of the attenuation constant, A, for 
of their desirable characteristics, namely, a linear | the most commonly used TE,, mode in a cylindrical 
decibel scale and potentially high accuracy.!. How- | guide. A good approximate’ functional relationship 
ever, their attenuation constant (decibels per radius | of A, the conductivity of the guide (nonmagnetic 
of travel) is a slowly changing function of frequency. | material), its radius, and the operating frequency is 
This variation is due to two factors: (1) the proximity | shown on the nomograph, and is the equation the 
of the operating frequency to the cutoff frequency, | nomograph solves. The effect of skin-depth is 
and (2) the frequency dependence of skin-depth, | negligible for the TMo, mode, and a nomograph for 
vhich alters the effective diameter of the attenuator. | this mode along with the TE,, mode for the case 
a | of infinite conductivity is already available.’ 
R. E. Grantham and J. J. Freeman, A standard of ; attenuation for microwave reRER 
rements, Trans. Am. Inst. Elec. Engrs. 67, 535 (1948 R. E. Lafferty, Piston attenuator chart. Electronics 21, No. 2, 132 (1948 
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This nomograph may be used in designing an 
attenuator to secure minimum frequency dependence 
over a given frequency range and to determine 
that dependence after its construction. To exemplify 
the latter, consider the radius a=2 em, the con- 
ductivity o=—1.510' mhos per meter, and the 
frequency f-=1» 10° eyeles per second; then ae 
3> 10° mhos, af 210° meters per second, and the 
extension to the A seale of a straight line passing 
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through the above points on the ao and af sea 
gives a value of 15.94 db per radius for A. In desig 
ing an attenuator, one simply chooses values of 
and ¢ to place the operating point preferably on t! 
right end of the ao scale and about the knee of | 
af scale. 


WasHineton, October 29, 1951. 
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able of the Zeros and Weight Factors of the First Twenty 
Hermite Polynomials’ 


Herbert E. Salzer, Ruth Zucker, and Ruth Capuano 


The chief use of this table of the zeros and weight factors of the Hermite polynomials 
is in the calculation of integrals over the interval [ ©, ©], when the integrand is either 
the product of ¢ and a polynomial, or may be closely approximated by « times a poly- 
nomial. The zeros and weight factors, z{/ and aj”, respectively, together with the aux- 
iliary quantities B)=a\” exp [(z7/”)*], which are useful in computation, are all tabulated 
here for the first twenty Hermite polynomials. The zeros z{? and 8! are given to 13 or 
more decimals, and the weight factors a,” to 13 significant figures Although other shorter 
tables have appeared, this present table will enable one to cope with problems requiring 
much higher degree and accuracy, both in problems involving direct quadratures and those 
arising in the numerical solution of linear integral equations for the range [— @ , |. Thus 
the use of this table in any direct quadrature can give exact accuracy as far as the 39th 
degree (the only inexactitude arising from the use of rounded values of z™ and at 


1. Introduction By considering e* f(r) K(y,z) as a polynomial in z, 
, and making use of (2) for y=r™,i=1,2,..., n, 
The main purpose of the table of the zeros zr; the approximate solution of (3) is reduced to the 
ind weight factors (“Christoffel numbers’) a?’ of | solution of a set of only n linear equations for the 
the Hermite polynomials is to provide the basic | (approximate values of) f(z"). These values of 
intities that occur in the quadrature formula f(x) can then be used in (2) to find f(y) for any 
, value of y. The advantage in the use of (2) instead 
e-® f(x)dzr=Dd> al f(x)+R,, 1) | of an equally spaced quadrature formula for solving 
i (3) is that here only n points are needed to give 
tied accuracy obtainable by approximating e* f(z) K(y,2) 
, as a polynomial of the (2n—1)th degree in x. Fora 
R. rj full description of the method, including examples, 

© 2" (2n)(2n—1)... (a +2) (+1) see A. Reiz [1], especially p. 4 to 10, 16 to 21. 

for some &, © &. ) [2, p. 101 102, 369].? 


When the integrand does not contain e~* explicitly, . , 
ae SP al ; Pi. wi 3. Important Properties of Hermite 

but behaves like « times a polynomial, instead of 2 

1), there is Polynomials 


F (x) dr ele? F @) dre Dales” F(ze). The Hermite polynomials may be defined by 
' Jj-e = the formula 
(2) 
: : 4 . n = d" r 
For dealing with quadratures of the type in (2), it H(z) =(—1)"« qe 6 )s (4) 
is convenient to follow the example of A. Reiz [1] " 
and to tabulate also the quantities BS’ =a’ exp | : 
r’’)*). From the form of the remainder term &,, | °F by the generating function 
1) is seen to be exact (save for the use of rounded | 
values of 2? and aS’) when f(z) is a polynomial | toe 3 (x)=. (5) 
of at most the (2n—1)th degree. P — ae | m 
2. Numerical Solution of Integral Equations | They are given explicitly by the formula 
Besides problems involving direct quadratures, ni} (—1)*n! 
there are those arising in the numerical solution of H,(x) (2z)"-*. (6) 


; f= Ai(n— 2k)! 
linear integral equations of the second kind for 


integrais over » ge D, «x amely, 
es ey eee ee : on The numerical values of the coefficients of //,(x) for 
: . the first 30 values of n are contained in H. E. Salzer 
o(yW=J(y) +r JS @)K (y,x)dz. \) | (3). The polynomials H(z) satisfy the differential 
| equation 





$ project was sponsored in part by the Office of Air Research, AMC, 
} 


’” 9 , L9 -~ 
cures in brackets indicate the literature references at the end of this paper. Hi, (r) 27H, (z) 2nil,(z) 0, ( ‘ ) 


lll 











the recurrence formula 
He (2) —2rH, (x) +-2nHT,,_ (2) =0, n=2, (8) 
and the differential-recurrence formula 
HT,’ (a 2niT,,_\(2). (9) 


The Hermite polynomials have the orthogonality 


property 


e~? H, (x) H,, (2) dz (10) 


wi2"n!lsm 


A most important consequence of (10) is (1), where 


the a!’ are defined by 


“H, (x) 


t—z;> 


l = i 


a" 
: H,'(z*°’) . 


dr. (11) 


are also expressible in the following form 


ah 
" 
I he a, 


[2, p. 344), [4, p. 311): 
}On+ly? 
gi 2"*t!n! 
ay ; (12) 
H,,'(2 
which, from (9), can be written as 
ri 2"-'(n—1)! 13 
ay 4 (1.5) 
n{H,_,(x‘"’)}? 


advancing differences. 


the latter form being more convenient for compu- | 


tation. 


4. Description of Tables 


This table gives the zeros 2;" 
* and 8{*’, for the first 20 Hermite polynomials. 


a, 
The absolute values of the zeros, |z;"’|, are tabulated 
to 15 decimals for n=1(1)13, to 14 decimals for 


n= 14(1)16, and to 13 decimals for n=17(1)20. To 
each nonvanishing |7{"’| there corresponds the two 


zeros 2x," tin’) and 2," x{|. The weight 
factors a{"’=a,"’,,, are tabulated to 13 significant 
figures. For a given n, zeros having the same ab- 


solute value are associated with the same weight 
factor. The auxiliary quantities 8; are tabulated to 
13 decimals. 


5. Method of Computation and Checking 


Since //,(z) is an odd or even poly nomial for n odd 
or even, respectively, it was convenient to change the 


variable to N=4z? and to calculate Y\{"’ =42,"", 
where the nonvanishing ,"’ are the zeros of poly- 
nomials A,,.»)(X), n odd, and B,.~)(.X), n even, of 
degree [n/2] in X. Thus the degree of the polyno- 
mials A).»)(X) and By,.,)(.X) did not exceed 10, and 
the coefficients were considerably smaller than those 
of H(z). A first approximation to ..;"’ was ob- 
tained from E. R. Smith [5], which gives the ; 


’, the weight factors 








Wha (; r/2) 


h, (x) =(—1)"e" ” 
~ dx" 


by taking twice the square of the latter. The po! 
nomials Ay» (XV) and By,» VY) were tabulated 
five or six values of \ in the neighborhood of ea: 
X;"’, at intervals of 0.0001 in = The zeros .\ 

were calculated by a formula for inverse interpolati: 
given in H. E. Salzer [6], and then were checked } 
direct interpolation in the polynomials A,,..)(.X) ar 
Bno(X), using the Gregory-Newton formula wit! 
The accuracy of the zen 
was determined from the approximate relatior 





Bi " 


. ‘ 0.0001 
‘ror AS oe (A). 
error in_LX, A ‘ 


where F(X;"’) denotes either Aj,»)(X) or 
Biuo( XY), and A is the difference between 
consecutive entries of FY) on either side of Y" 

The values of 2{ were found from 2" 


two 





l ; : 
=sv-Ay"’, and checked back by squaring to get 





a: — 
Af? m4n° Additional checks upon the 2;"’ were 
pe ‘rformed by calculating 
S » n(n—1) 
St ym)? 
— qv, » ? (14) 
i=l - 
n l ! 
' in )t as ' 
a’ 2{=(—1)"-! »\y for n odd, (l5a 
im] ((m-+ 1)/2)! j 


where x’ denotes the omission of the factor 2{" =0, 
and 

n ' 

gz i’ =(—1)*" I for n even. 15b) 


inl (Wye 
The weight factors a"? were calculated from (13), 

using 

) for n odd, 


H,, (z 3" ) B, a~—1) 9 (A ,” 


and 


H(z) =X?’ 


Ajin-y py (X,”’) for n even. 


Formula (12) was used to check the a;"’, in the form 


ee rey het ned, ist 
. ! , : 
a” ==. ——» forn odd, i#=(n+1) (12a 
; { A" e°o 2 \« XY; 2 . 2 | 
and 
$9On-3,1 
: ri 2"-8n! 
a!) = — ; > » for n even. 12b 
X”( Beak) } : 


Additional over-all checks upon a;") were performed, 


using (1) for P(xz)=1, 2’, zx‘, , at, 2°9<2n—1: 


zeros of /9 2”*' and in a few cases for P (z) gets? 2" 4 2p. 
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pens that that integral is known to be ya e7"8J,(/ 3 

where J,(t) denotes the modified Bessel function of 
P(r) —2*. The choice of these different powers order zero (see [26, p. 394]). As the integrand 
for P(r) was to ensure that every significant | e~".J)(f) is even, and as af JI.( Ss”) = a, as 

re of the a”? was covered in the checks. J,(—2x;"), from (1), whenever n is even, 

The auxiliary quanitites 8;"’ were computed by 


m+! the left member of (1) equal to r( k + : ) 


* oo 


¢ OJ) dt~ da!’ Jz’), (16) 





ap ong weight factors a" by es", where 


e” “were obtained from the NBS tables of the 
yonential function (7) and were checked by the 
ation log B{"’=log a{"’+ X;"'/4, where the loga- 
hms were obtained from the NBS tables [8]. 

Che functions 2", a{"’, and B;"’ were calculated 


/0 


where the summation covers only the positive values 
of z\". In the present example the error in (16) 
(from the term R, in (1) and the inequality 
Ji” (x)! <1) does not exceed 10>" for n as small as 10. 
about two more places than are tabulated, and | For the right member of (16), Jo(x{) was computed 
available in manuscript form. All entries given | to as many places as were sufficient to obtain 
re, which were obtained by rounding the manu- | a Jo(xi”) to 13 places. The J)(2{") was calcu- 
ript table, are correct to within about a unit in the | Jated by interpolation in the Bessel function table 
t place. | of the Harvard Computation Laboratory [27]. The 
The functions in this present table were compared | following are the results of these computations (the 
with those in three much shorter similar tables given | extra fourteenth place in a!'?Jo(2{"°) was carried to 
by R. E. Greenwood and J. J. Miller [9], A. Reiz [1], | avoid the cumulation of rounding errors) 
and Z. Kopal [10], and the few errors of more than a ‘i 








init in the last place in those three tables were noted 
and reported to the journal “Mathematical Tables ; Jo(x) a1 Jo(a..(10)) 
and Other Aids to Computation.”” Also, the values 
of X, for n ? upon the worksheets were compared 
, , ones enum 4s sneem @nnes anne 
with those in J. Barkley Rosser [11], with perfect : - eanes canon = : pen oe — 
‘ “ Oot . : . he * * vo | 
agreement. ‘Finally, the. zeros in the table of the 8 0. 36515 41522 17 0. 01236 93757 8926 | 
Harvard University Laboratory [28] were checked by mT) 0. 06452 02154 9 0. 00008 66923 1312 | 
calculating V2\x,""|, to give complete agreement. 10 0. 37050 3331 0. 00000 28308 0582 | 
6. Illustration 10 5 
So a” Jo(z;") =0.78515 05503 338. 
teb 


\s an example of the use of (1) for quadratures, | 
suppose that one wished to evaluate numerically, say | This answer differs by only a unit in the thirteenth 
place from the true value of vz e7'* J,(1/8) which 
Je was found (employing the power series for J,(t)) to 
denotes the Bessel function of order zero. It hap- | be, to 13 decimals, 0.78515 05503 339. 


@ 


to 13 places, the integral e-"Jo(t)dt, where Jo(t) 


7. Table of the Zeros and Weight Factors of the First Twenty Hermite Polynomials 


zr a" a 
n=l n=l n=l 
0. 00000 00000 O0000 1. 7724 53850 906 1.77245 38509 055 
ne=2 n=2 n=2 
n 0. 70710 67811 86548 0. 88622 69254 528 1. 46114 11826 611 
n=3 n=3 n=3 
N 0. 00000 00000 00000 1. 1816 35900 604 lL. 18163 59006 037 
1. 22474 48713 91589 0. 29540 89751 509 1. 32393 11752 136 
n=4 n-4 n 4 
) 0. 52464 76232 75290 0. 80491 40900 055 1.05996 44828 950 
1. 65068 01238 85785 0. (1)81312 83544 725 1. 24022 58176 958 
n=5 n=5 n-—5 
0. 00000 00000 00000 0. 94530 87204 829 0. 94530 87204 829 
0. 95857 24646 13819 0. 39361 93231 522 0. OS8658 09967 514 
2. 02018 28704 56086 0. (1)19953 24205 905 1. 18148 86255 360 | 
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0. 


3. 02063 
. 88972 48978 


n 6 
13007 74119 
33584 90740 
35060 49736 

n-—7 


00000 00000 


. 81628 TRS2S 


67355 16287 
65196 13568 
ns 
38118 60902 
15719 37124 
8165 67566 
03063 74202 
n-— 9 


00000 ODO00O 
72355 10187 


. 46855 32802 


26658 O5845 
19090 32017 


n- 10 


34200 13272 
03661 OS8207 
75668 36492 
53273 16742 
13615 9LISS8 


nell 
00000 00000 


65680 95668 
32655 70844 


2. 02594 SO15S8 


78329 00097 
66847 O8465 


n- 12 


31424 03762 


. 94778 83912 
. 59768 26351 


27950 70805 
70251 


nn 13 


00000 00000 
60576 38791 
22005 50365 
85310 76516 


. 51973 56856 


. 87871 
1. 47668 
2. 09518 32585 


24660 89783 
10133 75961 


n—-14 


20174 55106 
37873 
7311 


. 74847 07249 
. 46265 69336 


30444 85704 


27617 
13697 
74492 


00000 
5SU65 
67471 
35233 


07322 
16780 
5843 
57244 


00000 
52838 
16668 
31843 
S152S8 


23705 
89514 
99882 
32790 
37738 


00000 
82100 
94933 
ORT RR 


eV oo 
81652 
59583 


54359 
10164 
52605 
01060 
20800 
69782 


00000 
71060 
90748 
01512 
78238 
72410 
78640 


7256 
2040 
4114 
o772 
8540 
0227 
7363 


0. 


0. 


n 6 


0. 72462 95952 
0. 15706 73203 
2)45300 09905 


nu=7 


16175 
72526 
58281 
12450 


0. 81026 
0. 42560 
154515 
3)97178 


n-8 


0. 66114 
0. 20780 23258 
(1)17077 98300 
3)19960 40722 


n=9% 


0. 72023 52156 
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Film Dosimetry of Electrons in the Energy Range 
0.5 to 1.4 Million Electron Volts’ 


J. Fleeman and F. S. Frantz, Jr. 


Two types of films were exposed to electrons in the energy range 0.5 to 1.4 million 


electron volts 


The film response was found to be linearly proportional to the dosage re- 


ceived and independent of the energy of the incident electrons in the measured energy 


range. 
paper absorber surrounding the film. 


The photographic density is further shown to be dependent on the amount of 
A step-wedge filter is suggested that would permit 


exptrapolation of densities to zero filter thickness. 


l. Introduction 


The recent advances in radioactivity leading to 
high-activity radiation sources has made the study 


of personnel monitoring increasingly important. For 


camma rays, many of the techniques used earlier for 
X-ravs have provided development guides. How- 
ever, because high-activity electron sources have 
only recently become available, personnel monitoring 
systems for electrons are somewhat behind those for 
gamma rays. Several techniques, including photo- 
graphic films similar to .the film badges used for 
X-ravs, have been suggested. In the film technique, 
the opacity of the processed film is a measure of the 
radiation received by it. However, each type of 
film requires a calibration curve relating the darken- 
ing of the film to the radiation exposure. 

Monitoring techniques do not require a high ac- 
curacy as the short-time biological responses for low 
exposures are not appreciable, and in addition, there 
are considerable individual variations in radiation 
response. Usually an accuracy of 15 to 30 percent 
is acceptable. It is therefore entirely feasible to 
determine the ionization produced in an air cavity by 
the radiation and from this to eompute the dose de- 
livered to the surrounding tissue. The theoretical 
and experimental factors involved in this computa- 
tion are not accurately known, but it is estimated 
that the over-all accuracy of the computation is 
within 5 percent. 

Because of the high attenuation in the body of the 
radiation from commonly available electron sources, 
the dose received near the surface of the skin is most 
important. 
averaging 7 mg/em®, there is no reason for measuring 
the dosage at lesser depths. Actually, because of 
electron scattering, the dosage may increase some- 
what beyond this depth, but for moderate electron 
energies personnel monitoring dosages are measured 
at this depth. 

Density-dosage curves for two films, Minimax 
Dental X-ray Film Extra Fast and du Pont Dosimeter 
type 552 single film packet,’ when exposed to elec- 

e work described here was sponsored by the Atomic Energy Commission 
vil Defense, Radiological Monitoring Instruments Specifications (Dec. 6, 


The Minimax dental film may be obtained from the Minimax Co., Chicago, 
The du Pont film was a specially prepared film packet used by the AEC 
no longer available 


However, as the skin has a dead layer | 


| 
| 





| 


trons of 0.5 to 1.4 million electron volts (Mev) and 
also to the beta rays from uranium, are presented in 
the present paper. 

The polystyrene used in this experiment has an 
average atomic number and density close to that of 
human tissue. In addition, 7 mg/cm? of polystyrene 
is used to simulate the dead layer of tissue. Conse- 
quently, the correlation experiment reported here 
relates the darkening of the film previously men- 
tioned when exposed to electron bombardment in 
terms of the energy dissipated in 7 mg/cm? of poly- 
styrene, which is equivalent to the dead layer of 
human tissue. It must be clearly understood that 
the dosage is not the energy dissipated in the film 
but in the polystyrene. The two are related by way 
of the relative stopping powers of the materials. ~ 


2. Apparatus and Procedure 


High-energy electrons produced in an electron 
accelerator * passed through a 20 10~*-em-thick 
aluminum window into an evacuated collimator, and 
thence through a second window into either an ion- 
ization chamber or into the film badge. <A shutter 
near the first window permitted accurate control of 
the film exposure. Aluminum baffles in the collima- 
tor served to reduce the amount of scattered radiation 
incident on the film or ionization chamber and to 
make the main electron beam parallel. <A carriage 
supported both the film-badge holder and the ion- 
ization instruments so that either could be accurately 
positioned in the electron beam (fig. 1). 

The extrapolation-type ionization chamber ® (fig. 
2) was arranged so that the plates were perpendicular 
to the beam and the spacing could be varied. The 
plate upon which the electron beam was incident had 
a thickness of 7 mg/cm’ of polystyrene. The collec- 
tor plate size (1.04-cm diameter) was determined to 
be well within the area of uniform electron-beam 
intensity. An FP-54 electrometer located directly 
below the chamber measured the ionization current 
to this collector by indicating the voltage drop across 
a fixed resistor in the grid circuit. 





‘E. E. Chariton and H. 8. Hubbard, Gen. Elec. Rev. 43, 272 (1940 
1G. Pailla, Radiology 20, 22 (1937 
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Cross-seitional drawing of the 
chamber 


Ficure 2 extrapolation 


A, Accelerating electrode; B, amber bushings; C, collecting electrode 
Gi, guard ring; RK, Bakelite rings; 8, phosphor bronze spring; I, Bakelite insulator 


It has recently been recommended ° that electron 
dosages be expressed in terms of energy absorbed per 
rram of tissue at the position in question. It has 
-~_ shown ’ that the energy, &,., dloahel per gram 
of material is related to the number, J,,, of ion pairs 
produced in a small air cavity per gram of air by 


E,. sW,., 


* Recommendations of the International Commission on Radiological Units 
(London, 1990), Radiology 356, 117 (1951). See also, U. Fano and L. §. Taylor, 
Radiology 55, 743 (1050 

'L. H. Gray, Brit. J. Radiology 10, 600 (1907 


where s is the stopping power of the material relati 
to air,” and W is the energy required to produce a 
ion pair (5.2 107" erg to within 3%). J, is dete 
mined from the potential, |" (volts), across the eri 
resistor, 4.56 10° ohms; the area of the colleeto 
0.855 em?*; the density of air at normal temperatu 
and pressure, 0.00129 g/cm‘; the actual air pressur: 
P (mm of Hg); the temperature, 7 (°A); the ele: 
tronic charge, 1.60 10->" coulomb; the time. t: an 
the plate separation, d (em). Therefore, 


(1.020)(5.20)\10°')V(760)71t 
(4.56) 10°V(0.855)d(0.00129)P(273)1.6000 
vt7y 
dP 


ie 
1840 


Values of V were obtained experimentally for sey 
eral different plate separations, d. The slope of the 
resulting curve gave V/d, and this, together with the 
pressure and temperature, gave the exposure rate, 
F,/t. 

Data were obtained at each 0.1 Mev in the range 
0.5 to 1.4 Mev. The electron dosage rate was first 
determined at a given accelerating potential by the 
ionization chamber. The chamber was removed 
from the beam, and film were exposed to different 
dosages by timing the exposures. lonization meas 
urements were then repeated to determine the con- 
stancy of the dosage rate. For changes in dosage 
rates of less than 10 percent, average dosage rates 
were used; for changes larger than 10 percent, in 
terpolated dosage rates were used. 

Film exposures were also made on a large uranium 
plaque by placing the film directly on the tight 
7 mg/em’ covering of the plaque. The dosage rate 
of the plaque had previously Roaie measured * to be 
20.1 (ergs/g)/hour at the surface of the covering by 
an extrapolation chamber method. 

All films exposed at a given energy plus a control 
film were processed simultaneously. They were de- 
veloped for 5 minutes in Eastman Liquid X-ray Devel- 
oper and Replenisher at 20° C, fixed for 5 minutes at 
20° C, washed for 30 minutes, and dried. Film den- 
sity was determined on a calibrated Ansco-Sweet 
densitometer, model 11, for densities below 2.5 and 
on an Ansco color densitometer model 12 for higher 
densities. 


3. Results and Discussion 


Figures 3 and 4 show the results of the density- 
dosage determination of Minimax and du Pont 552 
film, respectively. In each figure, two sets of 
curves are presented. The upper set are the data 
for the density-dosage relationship, using mono- 
energetic electrons, and the results are plotted, 
using the various values of the incident energy of 
the electrons as a parameter; the lower curve is 
the density-dosage relationship, using the uranium 

* Computed values of s vary with electron energy, but 1.02 represents a mear 
value that is probably accurate to +2% for the energy range of interest here 


* We are indebted to the New York Operation Office of AEC for the loan « 
this standard and for its calibration constant 











i. An examination of the data for mono- 
retic electrons seems to indicate that the re- 
se of these films is independent of the energy. 
rder to establish this more accurately, the 
itv was plotted against the energy by using 
ral values of the dosage as a parameter (50, 100, 
By the method of least squares, 

as determined that to within the accuracy of 

results the density-dosage relationship is inde- 

dent of the energy in the energy interval 0.5 

1 Mev. These results can understood in 
following qualitative manner: The dosage to 
ch a film has been exposed is essentially a meas- 
of the amount of energy that has been expended 
the electron in its traversal through a film. The 

‘in cause of this energy loss results from the in- 

istic seattering of the electrons by the atoms of 

photographic emulsion. However, this proc- 
is the mechanism that accounts for the darkening 
film. Consequently, to a first approximation, 
ne expects a linear dependence of density upon 
suge for electrons of sufficient energy to penetrate 
film completely. 

Pele '° has made the most intensive theoretical 
study of the action of X-rays and electrons on 
photographic film. He shows that the form of the 

nsitv-dosage relationship is 


5 “rs or 
loO ergs/g). 


be 


D-DD, 


where 1) is the saturation density, and r is the dos- 


in roentgens. The quantity ¢ is a function of 
the mass of the undeveloped grains, the number of 
quanta required to affect a grain, and in the case 
of X-rays, the mass absorption coefficient of the 
silver-bromide grains. Expanding the quantity 
in the exponent, it will be observed that for small 
values of (er) one must expect a linear function for 
the density-dosage relationship. In the low-density 
region, it will be observed that the density does 
indeed exhibit this linear relationship very well. 
This independence of the response of the film to 
electron energy and the linear dependence upon 
energy loss in the film is of the type to be expected 
from Pele’s results. As shown by these curves, 
the slope of the curve is determined only by the 
tvpe of film used, that is, upon grain size, ete. For 
the du Pont 552 film only, there exists a slight amount 
of curvature for densities greater than 1.0. The 
saturation density for this film has been determined 

be roughly 3.5; thus, the curvature that is ob- 
served results from the saturation properties of 
the film. For Minimax film, this curvature is not 
‘bserved in the density range zero to six. In the 
low-energy region where the incident electrons are 
stopped in the film, the response of the film can be 
expected to be energy dependent." 

lhe lower curves of figures 3 and 4, respectively, 
show the response of Minimax and du Pont 552 film 
vhen exposed to uranium. These curves indicate a 


mire 


R. Pele, Proc. Phys. Soe. 52, 523 (1945 

F. Baker, E. G. Ramberg, and J. Hillier, J. Applied Phys. 18, 450 (1942); 
Horries, Z. Physik 199, 408 (1942); Physik. Z. 43, 190 (1047); H Nissen, 
k 192, 573 (1044 
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slope that is roughly one-half that of the correspond- 
ing curves obtained by using monoenergetic electrons 
in the previously mentioned energy interval. This 
difference in slope results from the difference in the 
energy spectrum to which the films have been ex- 
posed. The electrons from the uranium plaque have 
a continuous spectrum with a calculated most prob- 
able energy at roughly 0.1 Mev and an average in 
the neighborhood of 0.45 Mev. In the decay scheme 
of U** there are six beta emitters of interest, UX,, 
UX,, RaB, RaC, RaD, and Rak." Emitters UX, 
and RaC each have two possible modes of decay. 
For UX, the beta spectrum of one of the transitions 
has a maximum energy of 0.205 Mev and is 80 percent 
probable; the other has a maximum energy of 0.11 
Mev and is 20 percent probable. For RaC, the 
end-point energies are 1.65 and 3.16 Mev with a 
probability of 80 percent and 20 percent, respec- 
tively. The remaining end-point energies are 2.3, 
0.7, 0.029, and 1.16 Mev, respectively. From these 
data, semiempirical beta spectra’ have been con- 
structed for each of these transitions to determine a 
composite spectrum. In this manner, the above 
values of the most probable energy and average 
energy have been determined. The film itself ts 
wrapped in two layers of paper whose total thick- 
ness is of the order of 25 mg/em*®. From the curve 
of the composite U** spectrum, it can be shown that 
the fractional number of electrons absorbed in the 
paper wrapping is of the order of 22 percent, so that 
the fractional loss of electrons due to their absorption 
by the paper covering is large for electrons from ura- 
nium and small for electrons from the monoenergetic 
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beam. Consequently, the reduced density of the 
film when exposed to the uranium plaque results 
from the loss of electrons by absorption, as well as 
the modification of the electron energy spectrum 
upon passage through the paper wrapping. 
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Figure 4. 


An auxiliary experiment was performed to deter- 
mine the dependence of the Minimax dental film 
only upon the energy spectrum of the incident 


electrons. In addition to the paper wrapping nor- 
mally present, aluminum absorbers of varying thick- 
ness were wrapped around the film, and the films 
were exposed to the uranium plaque until each re- 
ceived a surface exposure equal to 41.9 ergs/g of 
tissue-equivalent material. A set of data was ob- 
tained by using monoenergetic electrons at 0.6, 0.9, 
and 1.2 Mev. The results of this experiment are 
shown in figure 5. The upper curves are the results 
for the density of the film with aluminum absorbers 
when exposed to monoenergetic electrons at the 
three energies. The lower curve is the result for 
the identical experiment using the uranium plaque. 
For zero aluminum absorber, the density of the ex- 
posed film after background had been subtracted is 
approximately 0.52 for the lower curve. This value 
is to be compared with the density of 0.45 obtained 
from the corresponding point in figure 3 when the 
background of 0.25 is subtracted. In figure 5, the 
density of the film for zero aluminum absorber at 0.6 
and 1.2 Mev is 1.10 and 0.87 respectively. Com- 
parison with 0.9 Mev is omitted as no density-dosage 





information had been taken at this energy for Mir 
max film. The densities from figure 3 at a dosage 
41.9 ergs/g are 0.90 and 0.95 for the energies 0.6 a: 
1.2 Mev, respectively. The curves of figure 5 fi 
the monoenergetic electrons give an average densit 
of 0.97 for zero aluminum absorber. This valu 
stands in good agreement with the value of 0.9 
obtained from figure 3. 
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A, Thickness of paper covering the film. Penetration by electrons; dosage 


equals 41.9 ergs/g 


It is clear from the foregoing results that a simple 
determination of the density is not sufficient in order 
to determine uniquely the dosage to which the film 
has been exposed. In addition to the density, a 
knowledge of the energy spectrum of the electron 
beam would be required. In order to circumvent 
this limitation, a step-wedge type of filter could be 
used. A plot of the filter thickness against film den- 
sities would then permit extrapolation to a filter 
thickness of 7 mg/cm. Because the thinnest filter 
would have to be quite thick in order to make the 
packet light tight and because beta rays and electrons 
produce very different types of attenuation curves, 
the proposed filter would not give a high order of 
accuracy, but in many cases it may be sufficient 
for personnel monitoring. 


The authors thank H. O. Wyckoff for the many 
helpful discussions during the course of the experi- 
ment; C. R. Horner, now with the Department of 
the Navy, who was chiefly responsible for the design 
of the extrapolation chamber and the FP-—54 elec- 
trometer housing; and J. A. Simpson, for the design 
and construction of the electron extractor. 


Wasuinoron, August 24, 1951. 
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[Thermal Converters as AC-DC Transfer Standards for 
Current and Voltage Measurements at Audio Frequencies 


Francis L. Hermach 


Thermal converters and associated e quipment that are used as ac-de transfer standards 
at the National Bureau of Standards for the precise measurement of current and voltage at 


power and audio frequencies are described. 


The standards and the equipment are primarily 
used to standardize a-c ammeters and voltmeters submitted to the 


Bureau for certification. 


The ac-de transfer may be made with these thermal converters at currents from | milliampere 
to 50 amperes, voltages of 0.2 to 750 volts, with an accuracy of 0.01 percent at frequencies 


from 25 to 20,000 cycles per second. 


The special tests to insure the required accuracy of the transfer standards are described, 


and the results are presented. 


A number of factors that limit the transfer accuracy of ther- 


mal converters have been discovered, and the results of special tests and theoretical work to 


evaluate these factors are discussed. 


The solutions, by an approximation method, of certain 


pertinent nonlinear differential equations governing the heating of a conductor by an electric 


current are given. 


1. Introduction 


lhe increasing use of electric energy for aircraft, 
induction furnaces, and induction heating, and the 
accuracy required in measurements in 
electronics, have led to increasing demands for the 
accurate standardization of ammeters and _ volt- 
meters at frequencies extending upward from power 
frequencies through the entire audio-frequency 
range. ‘To meet these demands, special instruments 
have been developed at the National Bureau of 
Standards for the measurement of current and 
voltage over rather wide ranges. They make use 
of thermal converters ' (often called thermoelements) 
like those incorporated in ordinary thermocouple 
instruments, but differ in the manner of reading 
and use. They may be used either directly to 
measure the ac-de differences of ammeters and 
voltmeters, or with a suitable potentiometer and 
to measure alternating currents and 
voltages. They were designed and are used pri- 
marily for testing electric instruments, at currents 
from 1 ma to 50 amp and voltages from 0.2 to 750 v, 
with an accuracy of 0.01 percent at frequencies 
from 25 to 20,000 e/s. 


greater 


access ries 


2. Transfer Principle 


The basic electrical units are defined in terms of | 


a concordant system of mechanical units and are 
realized by absolute electrical measurements carried 
at national standardizing laboratories to fix 
the value of groups of standard cells and resistors. 
These standards are used in conjunction with a 
potentiometer to make measurements of direct 
current, and power. It has been known 

long time, but not sufficiently realized, that 
measurement of the corresponding alternating- 


out 


voltage, 
for 
the 


posed AIEE definition 30.89.040: A thermal converter is a device that 

ts of one or more thermojunctions in thermal contact with an electric 
w integral therewith, so that the electromotive force developed at its 
terminals by thermoelectric action gives a measure of the input current 
ater. 


| 
| 
| 


higher accuracy than the first, 


current quantities depends fundamentally on certain 
standard types of electric transfer instruments that 
ideally have the same response on direct and alter- 
nating current. If the instruments are calibrated 
on direct current at the time of each use and if the 
precision of reading is suitably increased, long-time 
stability, freedom from drift, small temperature 
influence, low losses, and other normally desirable 
characteristics become of secondary importance, 
and the instruments may be primarily designed and 
constructed to have the best possible frequency 
characteristic. For difference measurements such 
instruments need not have long scales in the ordinary 
but rather may be designed for increased 
precision of reading by methods that would not be 
applicable to instruments intended for general 
service. 

Special electrodynamic transfer instruments have 
long been used at the Bureau for accurate a-c 
measurements. These transfer standards are used 
largely to standardize (test) other instruments 
submitted to the Bureau for certification. There 
are two distinct types of such tests. For the 
“straight a-c test”, the instrument under test and 
a transfer standard are connected to measure the 
same alternating electrical quantity (current, volt- 
age, or power), which is adjusted to produce the 
desired deflection of the test instrument. The 
response of the standard instrument is observed, 
then the standard is transferred to direct current. 
The direct quantity is adjusted to give the same 
response of the standard and is then measured with 
a suitable potentiometer and accessories. For the 
so-called ‘“‘ac-de difference test’’, both instruments 
are connected to measure the same quantity first 
on alternating and then on direct current, which 
is in each case adjusted to give the same deflection 
of the test instrument. From the averaged difference 
in the response of the transfer standard, the ac-de 
difference of the test instrument is computed. 
The second type of test can be made with somewhat 
and in conjunction 


sense, 
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with a test on direct current gives more information 
about the performance of the test instrument. 
Because the ac-de difference of an instrument de- 
pends upon geometrical factors that are relatively 
permanent, subsequent tests for checking the con- 
stancy of calibration need generally be made only 
on direct current. Occasionally transfer tests are 
made by using a selected low frequency (such as 
60 c/s) as the reference in place of direct current. 

An example of the accuracy with which such 
tests are made, laboratory standard instruments, 
with scales 12 in. long, are regularly standardized 
and certified at the Bureau to 0.05 of a scale division. 
For a 150-division instrument this requires an 
accuracy of measurement of 0.03 percent or better. 


3. Development and Description of Electro- 
thermic Transfer Standards 


3.1. Choice of Standards 


Several types of instruments were considered in 
choosing transfer standards to cover the full audio- 
frequency range. Electrodynamic instruments have 
been highly developed [1, 2, 3]? at the Bureau for 
measurements at the commercially important power 
frequencies and can be used with proper corrections, 
up to about 2,000 c/s. However, unavoidable in- 
ductance and stray capacitance errors limit their 
useful frequency range and present little hope for the 
desired extension. Electrostatic instruments have 


been carefully studied [4, 5] at the National Physical 


Laboratory. They are useful over wide frequency 
ranges, but because of their low-torque—weight ratio 
at low voltage are not readily adaptable for either 
the low voltage or the current ranges desired in this 
application. Electrothermic instruments, which use 
an effect produced by the heating of a conductor 
carrying a current to be measured, seemed more 
promising. Of the several kinds, thermocouple in- 
struments, in which the temperature rise of the con- 
ductor (heater) is measured by a _ thermocouple, 
seemed most feasible. Straight-wire heaters allow 
a wide frequency range, and the thermocouple 
measurements permit high precision of reading. The 
usual disadvantages of poor stability and large tem- 
perature influence are not important in their use as 
transfer instruments, and their low overload capac- 
ity is not a serious limitation in careful laboratory 
work. Fortunately, good thermal converters are 
commercially available. They have been used for 
measurements at frequencies up to about 65 Me, with 
a d-c millivoltmeter as the inliianter, and are gener- 
ally considered in the \%- to 2-percent accuracy 
classes. So far as is known, they have not been 
previously studied for use as transfer instruments at 
the frequencies and the accuracies needed in this 
application. 

For these reasons, commercial thermal converters 
were selected for study to determine their value as 
ac-de standards in a wide-range audio-frequency 


? Figures in brackets indicate the literature references at the end of this paper 





transfer voltmeter. A careful experimental a; 
theoretical study of their performance led to {| 
purchase of additional commercial thermal convert; 
of suitable ranges and the design and construction . 
equipment for using them as transfer standards f 
voltage and current. 


3.2. Description of Transfer Standards 


Each thermal converter purchased in the milli 
ampere ranges is of the vacuum type, mounted in a) 
evacuated glass bulb with supports of copper wir: 
embedded in the glass. The conductor, heated by 
the current to be measured, is a short straight wir: 
generally less than 4% em long and often less than 
0.001 in. in diameter. The heater alloy and dimen- 
sions are chosen to give a temperature rise of about 
200° C at the center where the hot junction of the 
thermocouple is fastened by a ceramic bead. This 
bead provides electrical but not thermal insulation 
between the heater and thermocouple circuits, with a 
coupling capacitance less than | yyuf and an insulation 
resistance greater than 50 megohms at 25° C. 

Because high current sensitivity is not required, 
the thermal converters in the 1- to 50-amp ranges are 
not evacuated. The hot junction of the thermo- 
couple is welded directly to the midpoint of a thin- 
walled tubular heater. These thermal converters 
are temperature compensated [6]. 

For current measurements, thermal converters in 
the series 1, 2, 5, 10, ete., with rated output electro- 
motive forces of 10 mv were purchased, with ranges 
from 1 ma to 50 amp. Thermal converters of 7.5- 
and 30-ma ratings, with appropriate series resistors, 
are used for voltage measurements. Figures 1 and 2 
show some of the thermal converters and the model 
A and model B voltmeter elements. It should be 
emphasized that the caption NBS on the nameplate 
sbown in figure 2 signifies only that the equipment 
was assembled for use at the National Bureau of 
Standards. All the thermal converters were pur- 
chased commercially, as were the components of the 
voltmeters. It should also be emphasized that the 
thermal converters shown in these figures are the 
transfer elements only. An indicator (to be dis- 
cussed in the next section) is necessary for all meas- 
urements, and a potentiometer and accessories are 
also necessary, except for transfer tests of other 
instruments. 

The model A voltmeter was developed as a proto- 
type instrument to meet an immediate need for 
voltage measurements. It consists of a 30-ma 
thermal converter connected with its heater in 
series with a commercial decade resistance box, 
modified to have two 1,000-ohm-per-step decades, 
and a 100-, a 10- and a 1l-ohm-per-step decade, all 
in series. As shown in figure 3, one end of the heater 
of the thermal converter is connected to the low-side 
terminal of the box, with the shield of the box con- 
nected to the other end of the heater and to the 
shield of the coaxial cable connecting the instrument 
to the circuit. The center lead of this cable is con- 
nected directly to the high-side terminal of the box 








Figure 1. Typical thermal converters. 


Mounted thermal converters and the model A and 
model B voltmeter elements. 


FiGuRE 2 


Like the electrodynamic transfer standards used at 
lower frequencies, the series resistor is adjusted in 
use to give the same nominal current level for each 
measured voltage. This was chosen at 20 ma, re- 
sulting in a voltmeter having a constant 50 ohms/v 
and an upper range of 400 v. The model A volt- 
meter, described in some detail in an earlier paper 
7| is now no longer used for routine measurements. 
For the model B voltmeter, a built-in 7.5-ma 
thermal converter was used, with a fixed resistor | 
having taps to give voltage ranges of 1.5, 3, 6, 7.5, | 
15, 30, 60, 75, 150, 300, 600, and 750v. The resist- | 
ince cards for this instrument were purchased com- | 
mercially, and the cards and thermal converter were 


123 


mounted in a suitable shielded box with the shield 
connected to the low-side terminal of the instrument, 
as shown in figure 3. For both instruments the 
effects of capacitance currents and of the self- 
inductance of the resistors were carefully considered ; 
the limiting factors are discussed in another section 
of this paper. 


Coaxial 
Cable 








Thermal Converter 





750 600 300 150 75 60 30 


Ficure : Circuit diagrams of the model A and model B volt- 


meter elements 


A. Model A voltmeter element, B, model B voltmeter element 


3.3. Indicator 


Because of its convenience and portability, a milli- 
voltmeter is generally used for measuring the emf in 
the usual thermocouple instrument. The precision 
of reading of the usual millivoltmeter is of course 
far too low for this application. For increased ac- 
curacy a potentiometer can be used to measure the 
thermocouple emf. However, if the thermal con- 
verter is used only as a transfer standard, the emf 
need not be measured in a “straight a-c’’ test; and in 
ac-de difference tests it is only the small change in emf 
between the a-c and d-c settings that is significant. 
The change and the full emf with which it is com- 
pared need be measured with relatively low accuracy, 
provided the change is referred to a highly stable 
base value. Thus high sensitivity and high stability 
but only moderate accuracy in the measurement of 
emf are required of the indicator. For these require- 
ments, a Lindeck potentiometer for providing the 
base value, used in conjunction with a suitable gal- 
vanometer, the deflection of which indicates the 
change, forms the ideal indicator. In a Lindeck 
potentiometer, based on Poggendorf’s second prin- 
ciple [8], the emf to be measured is balanced by an 
adjustment of the current through a fixed resistor. 
At balance the voltage drop across the resistor, 
which is the product of the potentiometer current 
and the resistance, is equal to the measured emf. 
The current is ordinarily measured with a milliam- 





meter, which generally sets the limit on the accuracy 
obtainable with this form of potentiometer 

Such a potentiometer has been incorporated in a 
panel constructed for the audio-frequency testing. It 
was designed to have adequate stability, freedom from 
changing extraneous electromotive force and the re- 
quired ranges, with sufficient precision of reading. Its 
circuit is shown in the central portion of figure 4 
The resistor, R, is made of manganin, and a special 
thermofree key [8] and copper binding posts are used 
in the electrothermally sensitive emf circuit. The 
components are mounted in the central part of the 
panel, as shown in figure 5, and are enclosed by a 
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© signifies coaxial lead with shield grounded 


grounded metal shield. No other thermal pre- 
cautions were found necessary. A _ reflecting wall- 
mounted galvanometer is used, with a voltage sensi- 
tivity (in the cireuit) of 1 mm/uy A No. 6 dry cell 
supplies the potentiometer current, which is adjust- 
able from 0.2 to 1.2 ma. At these current levels, 
repeated tests have shown that the steady drift in the 
voltage across the four terminal resistor is less than 
0.02 percent per hour, and that superposed fluctua- 
tions are much less than 0.01 percent. The four- 
terminal manganin resistor has taps at 0.1, 0.5, 2.5, 
5, 10, and 25 ohms, giving milivolt ranges, at a 
current of | ma, of these same numerical values 
The lowest range is used in transfer tests as a quick 
check of the potentiometer. With the thermocouple 
of the thermal converter connected to the “TC” 
posts of the potentiometer, but with no current 
through its heater, a galvanometer deflection of 10.0 
em to the left on the seale with the Lindeck set for 0.1 
mv serves as a valuable partial check of the required 
accuracy and correctness of the Lindeck circuit 
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Schematic diagram of circuits used for testing instruments 


3.4. Use of Instruments 


These transfer standards are used almost solely for 
standardizing other a-c instruments. The necessary 
switching and control circuits for such tests have been 
incorporated in the panel shown in figure 5. In 
addition to these circuits and the Lindeck potenti- 
ometer, the panel contains the impedance-matching 
transformers for the high-voltage circuits used in 
testing voltmeters and the high-current circuits used 
in testing ammeters. The present transformers pro- 
vide ranges up to 50 amp and 800 v, with a nominal 
input voltage of 150 v anda power level up to 100 w 
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Motor-generator sets and voltage stabilizers are 
available for tests at 60 c/s and audio-frequency 
oscillators and power amplifiers for tests at other 
frequencies. Batteries, controlled by suitable ad- 
justable resistance voltage dividers and series resis- 
tors, are used for the necessary direct-current sources. 
The potentiometer, standard cells, volt boxes, and 
resistors are those normally used for instrument 
testing and are periodically standardized to insure a 
continued accuracy of considerably better than 0.01 
percent, 

A schematic diagram of the major circuits and 
equipment on this panel is shown in figure 4. For an 
a-c test of an ammeter, the transfer thermoelement 
is connected to the AC-DC binding posts and its 
thermocouple to the TC posts on the panel. The 
test instrument is connected to the AC posts, and a 
suitable four-terminal standard resistor is connected 
to the DC posts for use in measuring the d-c current 
by means of a potentiometer. For low currents at 
the higher audio frequencies, shielded leads with the 
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lds grounded are used to minimize errors due to 
y capacitance currents. The alternating current 
\justed for the desired deflection of the test instru- 
it, and the Lindeck potentiometer is adjusted for 
‘on-scale”’ deflection of the galvanometer. The 
ndard instrument is then switched to direct cur- 
which is adjusted to give the same galvanometer 
lection and is then measured with the external 
entiometer. The direct current through the 
iter of the thermoelement is then reversed and the 
measurement repeated 


Instrument testing panel, including Lindeck 


t le ne nt 


Figure 5. 


For an a-c test of a voltmeter, a similar procedure 
s used. The transfer standard is connected to the 
\V-DV receptacle on the panel with a coaxial lead, 
vith due care to connect its thermal converter to the 
rounded shield of this lead, and the thermocouple 
connected to the TC posts on the panel. The 
est instrument is connected to the circuit with a 
pecial twin lead consisting of two coaxial lines per- 
nanently paralleled at the instrument and connecting 
0 the two AV receptacles on the panel. Similar 
eads of the same length are used to connect the high 
ule of a volt box for use with the external poten- 
ometer to measure the direct voltage, to the two 
D\ receptacles The use of these leads eliminates 
the effect of any stray magnetic field or lead resistance 
on the measurements. 
This simple procedure is reliable only if the drift 
of the response of the standard instrument is negli- 
bly small in the time taken for a series of such 


ich reversed d-c readings are common practice for highest accuracy with all 
f ac-de instruments 


readings. This is fortunately true for each of the 
thermal converters used. After an initial warm-up 
period of a minute or so, the maximum observed 
drift under those laboratory conditions of any of the 
thermal converters tested was less than 0.5 percent 
per hour, and the average less than 0.05 percent per 
hour. If it were necessary the effect of significant 
drift could be practically eliminated by using a deflec- 
tion method as outlined for transfer tests, calibrating 
the standard instrument on direct current at the 
nominally correct current or voltage before and after 
each a-c setting. 

For ac-de difference or transfer tests, which are 
more frequent and usually more important than 
straight a-c tests, the test and standard instruments 
are connected in series to the AC-DC posts ¢ in the 
case of ammeters, or in parallel with suitable coaxial 
leads to the AV-DV receptacles in the case of volt- 
meters. Suitable short leads are used to connect the 
two AV receptacles together and to connect the two 
DV receptacles together. The current or voltage is 
adjusted for the desired deflection of the test instru- 
ment, and the Lindeck potentiometer is adjusted for 
a deflection near the null position at the center of 
the galvanometer scale. The two instruments are 
then in quick succession connected to alternating, 
direct, reversed direct, and then alternating current. 
The current is adjusted to give the same deflection 
of the test instrument in each case. The resulting 
deflections of the galvanometer of the Lindeck poten- 
tiometer are observed with the potentiometer con- 
trols undisturbed. This procedure practically elim- 
inates any effect of drift in either instrument, and 
errors in repeating the settings of the instrument 
under test can be made very small by using a low- 
power microscope and setting the image of the 
instrument pointer in coincidence with a cross-hair 
in the eyepiece. The control circuits have been 
arranged to facilitate switching, and with two 
observers the procedure is rapid even with moder- 
ately fluctuating voltages as only one of the two 
instruments is set to a given deflection. 

From the average of the deflections on alternating 
and on direct current, the ac-de difference of the test 
instrument is readily computed by a simple formula.® 


4. Evaluation of Transfer Performance 


The ac-de difference of each transfer standard 
must be known to the full accuracy desired in the 
measurements to be made. As, in a sense, there are 
no absolute transfer standards, such an evaluation 
essentially consists of four steps. 

1. Choice of a type of instrument that theory in- 

dicates is suitable over the desired ranges. 

2. Study of all known effects that can cause such 

an instrument to depart from this ideal. 

3. Construction and use of an instrument in such 

a way that these effects, by computation or 
direct tests of the separate components, are 
negligible or known. 


‘One terminal of the heater is connected to the post that is grounded on 
alternating current 
’ The formulas used in computing the results of the a« 


developed in appendix | 


and transfer tests are 
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4. Comparison of actual transfer performance of 


the standard with that of some other instru 
ment, preferably of a quite different type, 
whose performance has also been evaluated 
The evaluation of electrothermic transfet 
standards was considered of primary importance 
because of the high order of accuracy demanded in 
this application. The unique facilities of the Bureau 
made this quite feasible, particularly for step 4, 
which is necessary to guard against errors from un 
suspected causes. Prototype instruments such as 
these, in which previously available ranges are 
greatly extended, cannot be directly compared with 
other types of standards over their full range. The 
evaluation therefore consisted of the study, compu 
tations, and special tests of steps 2 and 3 during the 
design and construction of these instruments and 
their panel, the intercomparisons of thermal con- 
verters of adjacent current ranges in the series, the 
comparison at low frequencies of selected thermal 
converters with the electrodynamic transfer stand 
ards at the Bureau, and the comparison of these 
thermal converters with certain portable hot-wire 
and electrostatic instruments also available 
A thermal converter can used as a transfer 
standard for alternating-current measurements over 
the range of frequencies for which the response char 
acteristic, =f(/),is the same as with direct current, 
where / is the output emf, and / the current through 
the heater of the thermal converter. For a volt 
meter element consisting of a resistor in series W ith a 
thermal converter, the magnitude of the effective im- 
pedance Z=\// must also be equal to the effective 
resistance on direct current, where V is the applied 
voltage In addition to these two requirements foi 
instrument testing, the fundamental principle that the 
standard and the instrument under test “see’’ exactly 
the same quantity must be carefully observed. 
The response characteristic of the types of thermal 
converters used were believed to be suitable at fre- 
quencies far higher than those used in this applica- 
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cation, and their direct-current and low-frequen 
responses were verified by special tests described i 
other sections. The effective impedance of eac! 
voltmeter element was approximately compute: 
from inductance and capacitance data supplied b 
the manufacturer of the resistance cards, and wa 
verified by special tests to be described later. Thy 
capacitance and resistance of the Lindeck potenti 
ometer to ground were measured and found to by 
1,600 wuf and 10,000 megohms at an ambient relatiy: 
humidity of 40 percent. Thus the computed stray 
currents are sufficiently small so that no appreciabl: 
fraction of the currents measured will be diverted 
through the thermocouple at audio frequencies wher 
one end of the heater is grounded. The 
shielded leads where necessary, the avoidance of » 
ground between the test instrument and the standard 
in the current circuits, and the use of the special co 
axial lines in the voltage circuits assure that the 
fundamental principle observed even at the 
smallest current and voltage in the present ranges 
The comparison tests of step 4 were planned to 
give maximum possible assurance against errors due 
to unsuspected causes. Each comparison consisted 
of repeated careful transfer (ac-de difference) tests 
as described in section 3.4. The 1-, 2- and 5-amp 
thermal converters were compared with the NBS 
electrodynamic transfer ammeter at selected fre 
quencies from 25 to 200 c/s. The model A voltmeter 
was compared with the NBS suppressed-zero electro 
dynamic voltmeter at frequencies of 30 and 60 c/s and 
voltages from 20 to 240 v. Both of these electro- 
dynamic instruments have been carefully studied 
{1, 2}. They can be read with a precision of better 
than 0.005 percent. The results of these compari- 
sons, corrected for the small known transfer errors of 
the electrodynamic instruments, are shown in table 1. 
They are estimated to be accurate to 0.005 percent or 
better. The results of these tests disclosed that the 
model A voltmeter and the electrodynamic voltmeter 
were in excellent agreement at these frequencies 
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y disclosed significant disagreements between the 
pere-range thermal converters and the electro- 
amic ammeter, which were considered to be due 
low-frequency error in the 1- and 2-amp thermal 
verters and a small error independent of frequency 
he 5-amp thermal converter. 

he model A voltmeter was also compared with 

» new high-grade portable pivoted electrostatic 

imeters, having scales over 5 in. long, at selected 

quencies from 30 to 20,000 c/s. For increased 
cision a low-power microscope was focused on 
pointer of the electrostatic instrument under 

‘. This instrument was always set to the same 

(lection and was tapped lightly before each read- 

Repeated sets of careful readings were taken 
each frequency, and all tests were at the full-scale 
flection of the electrostatic instrument. Because 
these repeated readings and the favorable condi- 
ns, the averaged results are estimated to be accu- 
te to 0.01 percent. The results, as shown in table 
verify the absence of capacitance errors in the 

thermal voltmeter. Direct comparisons at the 
higher frequencies could not be made at lower 
voltages, but calculations indicate that skin effect 
and the effect of inductance should be considerably 

ss than 0.01 percent at 20,000 c/s. 

Similarly, the 5-amp thermal converter was com- 
pared with a high-quality portable 5-amp hot-wire 
ammeter at frequencies from 25 to 20,000 c/s. The 
results, as shown in table 1, indicate excellent agree- 
ment at frequencies of 1,000 ¢/s and above. 
at lower frequencies showed very puzzling discrep- 


Tests 


ancies, Which, in a separate investigation, were 


traced to vibration errors © in the hot-wire instru- 
ment. For that reason the results at frequencies 
less than 1,000 c/s are not shown in this table. 

\ plus sign in table 1 indicates that with equal 
response of the thermal converter on direct and 
alternating currents the response of the other instru- 
ment was greater on alternating than on direct 
current, 

Each thermal converter was compared with the 
next higher- and lower-range thermal converter in 
the series at two currents and at selected frequencies 
from 25 to 20,000 c/s. The results of each of these 
comparisons at the higher of the two currents at 
which tests were made are shown in table 2. They 
are estimated to be accurate to somewhat better 
than 0.01 percent, with very little possibility of 
systematic error. A plus sign in the table indicates 
that with equal response of the lower-range thermal 
converter on alternating and on direct currents, the 
response of the higher-range thermal converter was 
greater on alternating than on direct currents. This 
series of tests was made to extend stepwise the 
results of low-frequency comparisons with the elec- 
trodynamic ammeter, and to indicate any high- 
frequency effects that might well be expected to 
differ in different thermal converters.’ The group 

Electromagnetic forces between the hot wire and its closely adjacent iron 

inting plate resulted in vibration of the wire because the plate was magnetized 

the damping magnet of the instrument This vibration caused additional 
ng and thus an error that was as large as 0.2 percent at certain frequencies 


The percentage errors due to any high-frequency effects would be expected 
independent of the current level and could thus be evaluated by these tests 


TABLE 2. Results of intercomparisons of thermal converters 


Comparative ac-de difference 
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* This 30-ma thermal converter was the one used in the model A voltmeter 

of tests disclosed small discrepancies independent 
of frequency in a few thermal converters, and ver- 
ified the absence of any high-frequency errors in 
the thermal converters tested. However, they could 
not be used to assign numerical values to the errors 
of individual thermal converters in the milliampere 
range as the differences were found to depend upon 
the current. The concluding series of tests to eval- 
uate these errors were comparisons of a shunted 
thermal converter with each of those pairs of thermal 
converters that showed discrepancies in the previous 
intercomparisons. The model A_ voltmeter was 
shunted with a bifilar resistor of computed skin effect 
and time constant for tests at currents from 1 to 40 
amps, and with a suitable high-quality audio-fre- 
quency decade resistance box for tests at lower 
currents. The results are shown in table 3 and are 
estimated to be accurate to better than 0.01 percent. 
A positive sign in this table indicates that with equal 


Results of comparisons of selected thermal converters 
with shunted model A thermal voltmeter 


TABLE 3. 


‘omparative ac-de difference 
Test 


current 


Thermal converter 
range 


100 c/s woes 
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response of the thermal converter on alternating and 
on direct current the response of the shunted model 
A voltmeter was greater on alternating than on 
direct current. The results of tests of the 1-, 2- and 
5-amp thermal converters, duplicating those with 
the electrodynamic ammeter, are not shown in the 
table. The average difference between eight such 
duplicate results obtained with this ammeter and 
with the shunted model A voltmeter was 0.004 per- 
cent, and the largest difference was 0.013 percent. 
A few additional check tests of these thermal con- 
verters with the shunted voltmeter were made at 
other frequencies but are not tabulated. 

The results of all these tests were combined to 
assign reliable values of ac-de difference, 4, to each 
thermal converter in the series from 1 ma to 50 amp 
and to the model A voltmeter. These values were 
less than 0.01 percent for the model A voltmeter and 
for all the thermal converters of ranges less than 100 
ma. The values for the other converters (rounded 
off to the nearest 0.005%) are shown in table 4. 


Observed percentage ac-de differences of thermal 
converters (to the nearest 0.005 percent) 


TaBie 4. 


(a) Differences independent of frequency 


Percentage of rated curren 


ma 
1 through 5% 
mo 
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(b) Differences dependent upon frequency 


Frequency (cycles per second) 


Test 
current 


Converter 
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® Leas than 0.01 percent throughout, 


The model B voltmeter, which was developed to 
give a wider voltage range and greater convenience 
of use than the model A, was compared with the 
model A voltmeter at rated voltage on each range 
up through the 300-v range and at 400 v on the higher 
ranges, at selected frequencies from 30 to 20,000 e/s. 
The results are shown in table 5 and represent 
directly the ac-de difference, 6, of the model B volt- 
meter. They disclosed significant but not unex- 
pected errors for the higher voltage ranges at the 
higher frequencies. 


TaBLe 5 


Range of 
model B 


Volts 


3 through 75 


1”) 
300... 


thermal voltmeters 


Results of comparisons of model A and model j 


ac-dc difference of model B voltmeter 


Voltage a ee ee 
25 1,000 5,000 
c/s c/s e/s 


, 
Rated 
do 


do 


15,000 | 20,000 
c/s c/s 





+0. 014 
+. OZ 


+0. 044 
+. 004 


ao wo 


7™) +).003 +0. 011 


*® Less than 0.01 percent, 25 to 20,000 c/s each range, 


As a result of these extensive series of tests, the 
transfer performance of each of this first set of ther- 
mal converters and of each of the two voltmeters 
was uniquely determined. Small but significant ac- 
de differences of both known and unsuspected causes 
were found. Corrections for these differences can 
readily be applied when these thermal converters are 
used as transfer standards. However, an investiga- 
tion of such differences is an especially important 
objective of a development such as this. As a result 
of such an investigation, a number of factors that 
limit the transfer accuracy of an electrothermic in- 
strument were discovered and evaluated, with the 
results given in the following sections. 


5. D-C Response 


In all measurements with a transfer standard the 
average of the d-c response of the standard for the 
two directions of current is taken as the reference upon 
which the a-c measurements are based. For all 
three types of such standards, the reversed d-c aver- 
age gives a better basis than does the response for one 
direetion alone as effects that lead to small differ- 
ences in response for the two directions of direct cur- 
rent are reduced to the second order of smallness 
when the average is taken. 

Reversal differences may be expected in a contact 
thermal converter because of the flow of heater cur- 
rent through the hot junction of the thermocouple. 
This can cause a voltage drop in the thermocouple 
circuit, which may lead to second-order errors be- 
cause of the nonlinearities of the thermal and electric 
circuits. This cause of reversal difference is mini- 
mized by the special constructions used in most 
high-range thermal converters, and is eliminated by 
the use of the small insulating bead between the 
heater and the thermocouple in the low-range ther- 
mal converters. Surprisingly, however, early tests of 
the first thermal converters purchased for the model 
A voltmeter showed appreciable reversal differences 
(up to 0.2%) even in these insulated-heater thermal! 
converters. It was soon realized that this was prob- 
ably due to thermoelectric effects in the heater cir- 
cuit. Peltier and Thomson heating of the heater can 
cause dissymmetry in the temperature rise of the 
heater, and thus cause differences in the emf of the 
thermocouple unless the hot junction is exactly at the 





point of the heater.* This supposition was veri- 
by an analysis of the data, which showed that 
differences were larger in those thermal converters 
) heaters having larger thermoelectric effects. 
ther verification was secured by tests of special 
rmal converters, constructed by one of the manu- 
urers * at the Bureau’s request, with the hot 
ection of each deliberately and considerably off- 
iter. The reversal differences were much larger, 
| the signs and magnitudes were found to be 
whly in agreement with calculations based on 
rmoelectric data for the alloys used. 

Early intercomparisons of these few thermal 

iverters indicated that no significant ac-de dif- 
ferences should be expected, even with rather large 
reversal differences from this cause, provided that 
the average for the two directions of current is 
taken as the d-c response. However, intercom- 
parisons and tests of the much larger number of 
thermal converters later purchased for current 
measurements showed that a good proportion had 
measurable ac-de differences, all of the same sign 
and independent of frequency. This suggested 
a d-c error, for which these thermoelectric efforts 
seemed the most likely cause. This was inves- 
tigated theoretically (see appendix 2) by consid- 
ering the temperature rise of a homogeneous con- 
ductor heated by an electric current and cooled 
solely by conduction to two relatively massive 
terminals; i. e. by neglecting, for simplification, 
radiation and convection losses and the heat ab- 
stracted by the attached thermocouple. 

With only the normal joule, or resistance, heating, 
and with the terminals of the conductor at equal 
and fixed temperature, the temperature distri- 
bution along the conductor is parabolic, and the 
temperature rise at the midpoint is, as Goodwin 
showed," equal to 

Vy? 


co (1) 
Spk 


where V, is the voltage drop across the conductor, 
and p and k are the electrical resistivity and ther- 
mal conductivity, respectively. In addition, how- 
ever, there is Peltier heating or cooling at the junc- 
tions of dissimilar metals (the two junctions between 
the heater and its two supports) and Thomson 
heating along each half of the wire. Unlike joule 
heating, these are dependent upon the direction 
of the current flow and can thus cause a dissym- 
metry of the temperature distribution along the 
conductor, which reverses when the direct current 
through the conductor is reversed. On alternating 
current, even at the low frequencies in which we 
are interested, the reversal of current occurs so 
rapidly that the thermal inertia of the wire prevents 
any such dissymmetry, and the temperature dis- 


* This dissymmetry has actually been utilized to measure the Thomson effect. 


* M. Rosenfield, of the Field Electrical Instrument Co., who first showed that 
the reversal differences were dependent upon the type of heater material. 

* See [6]. In addition to this classic paper, Goodwin has written an excellent 
eries of articles on ‘Thermal problems relgting to measuring and control serv- 

es"’, which appeared in Weston Eng. Note’, 1948-50 (Weston Electrical Instru- 
ment Corp., Newark, N. J.). 





tribution is unchanged by these thermoelectric 
effects. Calculations show (see appendix 2) that 
with direct current the temperature rise at the 
midpoint of the conductor, where the hot junction 
of the thermocouple is attached, is unaffected by 
the Peltier heating, for completely symmetrical 
construction of a thermoelement. They show, 
however, that it is changed by Thomson heating, 
and becomes, approximately," 


6, of 1 -; (7 )} (2) 


where 6, is the midpoint temperature rise in the 
presence of Thomson heating, in ° C; 6, the rise 
without Thomson heating; ¢, the Thomson voltage 
coefficient in volts-degree ~', V,, the voltage across 
the conductor in volts. 

This was derived by considering ¢ a constant, 
whereas both experimental results and the electron 
theory of metallic conduction indicate that o— BT, 
where B is a constant and 7 is the absolute tempera- 
ture. However, a solution by a_ perturbation 
method,” of the nonlinear equation that results when 
o=BT, shows that the above results give a suffi- 
ciently close approximation for the temperature rise 
of about 200 deg C encountered in typical thermal 
converters. Thus eq 2 shows that an ac-de difference, 
or transfer error, should be expected in a thermal 
converter, and enables the approximate calculation 
of this error. Such calculations have been made for 
some common metals and some of the alloys used as 
heaters in thermal converters, on the assumption of 
a 200 deg C rise at the midpoint, and a voltage drop 
of 0.2 v at rated currents." The results, converted to 
ac-de difference for the same temperature rise, rather 
than temperature difference for the same current, are 
shown in table 6. For the computation of these 
results, the Thomson coefficient ¢ was considered to 
be equal to B(@)+-2/30,+-273), where B is the second 
derivative of the characteristic curve of emf versus 
temperature for a thermocouple composed of the 
indicated metal and lead."* 4 is the ambient temper- 
ature indeg C. The results show approximately the 
errors to be expected for these materials. 


TaB.e 6. Calculated ac-de difference at rated current due to 


Thomson effects in thermal converter heaters 


ac-dc 


Metal difference 


Cc Percent 
+0.8xK10-% <0. 005 
Platinum 3.2 013 
Manganin —0. 08 <. 005 
Constantan (60 Cu-40 Ni) 7.3 O64 
9 Ni-10 Cr 3.2 013 
85 Pt-15 Ir 2.1 005 


Copper 


'! At other points along the conductor the average temperature rise for the two 
directions of direct current is also changed by Thomson heating. 

” This method was suggestéd by Chester Snow. 

% From eq 1, a series of thermal converters having similar thermocouples and 
the same rated output emf at rated current, and having heaters that obey the 
Weidemann-Franz law will all have the same temperature rise and the same 


voltage drop across the heater. The temperature rise is about 200 deg C, and the 
voltage drop about 0.2 v at rated current for most thermal converters. 

“ This weighted average absolute temperature gives results in closer agreement 
with those computed by the more exact method than does the unweighted average, 
60+ 1/204+-273. 





The results may be compared with the observed 
ac-de differences shown in table 4. From information 
supplied by the manufacturers, the 1-, 2-, 5-, 10-, 20-, 
50-, and 200-ma thermal converters have heaters of 
either carbon, nickel-chromium, or _platinum- 
rhodium-ruthenium alloys. The observed ac-de dif- 
ferences of these converters were less than 0.01 per- 
cent. The 100-, 500- and 1,000-ma thermal convert- 
ers had heaters of constantan (Advance). The 
observed differences for these converters (at rated 
current) were close to the calculated value of —0.06 
percent. The 1- and 2-amp converters, of platinum- 
iridium, had observed differences (from table 4,b) 
less than 0.005 percent, except at low frequencies, 
also close to the calculated values. The 5-, 10-, 20- 
and 50-amp converters, of a special platinum—3- 
percent-copper alloy, had differences of —0.010 to 

0.015 percent. Thus the agreement between the 
computed and observed values serves as excellent 
verification of this source of error in transfer thermal 
converters. Very fortunately, the error is signifi- 
cantly large only for materials having unusually 
large thermoelectric effects, such as constantan. 
However, this does lead to the interesting conclusion 
that for operation at a low temperature a thermal 
converter should have large thermoelectric effects 
in the thermocouple, but that for the highest transfer 
accuracy it should have very small thermoelectric 
effects in the heater. 

Additional experimental verification of this cause 
of transfer error was secured by constructing” at the 
Bureau two 5-amp thermal converters differing only 
in that one had a heater of manganin and the other of 
constantan. These were 5-amp thermal converters 
with a small Borax bead between the hot junction of 
a copper-constantan thermocouple and the center of 
the heater wire, to eliminate contact effects. Careful 
ac-de difference tests with the standard electro- 
dynamic ammeter as the reference established agree- 
ment within 0.005 percent for the manganin element 
and indicated an ae-de difference of —0.06 percent 
for the constantan element at an emf of about 5 mv. 

As a result of this work, a second set of thermal 
converters in the 1- to 1,000-ma range was purchased 
under specifications prohibiting heaters of constantan 
or similar alloys. All of these converters were found 
to have ac-de differences less than 0.01 percent at 
audio frequencies. 


6. Low-frequency Response 


It has long [10] been known that there is a low- 
frequency limit to the correct rms response of a ther- 
al converter, for reasons analogous to those causing 
a similar low-frequency limit in all other types of 
transfer instruments. Below this limit the response 
characteristic E= (1) is not independent of frequency 
(except for symmetrical square waves), and may be 
expected to be in error by an amount that is inversely 
proportional to the square of the frequency. As 
table 4, b, shows, such an error was observed for the 
1- and 2-amp thermal converters. 


’ These thermal converters were constructed by E. 8. Williams. 








When a constant direct current is applied, the 
heater of a thermal converter reaches a constant 
temperature. When an alternating current is ap- 
plied (after the transient state is over) the tempera- 
ture of the heater varies cyclically about a mean 
value, with an amplitude that decreases with  in- 
creasing frequency because of the thermal capacity of 
the element. At frequencies high enough so that the 
cyclic variations are negligible, the temperature rise 
of the heater is independent of frequency (provided 
that at very high frequencies the skin effect and 
electromagnetic radiation are negligible and that 
standing waves do not alter the temperature distribu- 
tion). However, at lower frequencies the mean 
value, to which the d-c instrument or potentiometer 
connected to the thermocouple responds, may not be 
the same because of nonlinearities in the thermal and 
electric circuits. These nonlinearities are caused 
principally by radiation losses, the dependence upon 
temperature of the electrical resistivity and the ther- 
mal conductivity of the heater, and the curvature of 
the emf versus temperature characteristic of the 
thermocouple. A theoretical evaluation of these ef- 
fects for the conditions in a thermal converter, in 
which most of the cooling is by conduction to the 
heater, is difficult, but does serve as a guide in deter- 
mining the approximate limiting frequencies for an 
arbitrarily chosen error and for determining which 
factors influence that error. Such an approximate 
analysis has been carried out by a perturbation 
method (see appendix 3). The results show that the 
average (time average) temperature rise at the mid- 
point of a conductor heated by a sinusoidal electric 
current and cooled solely by conduction to massive 
terminals is approximately 


ee 
Bac = ge [1 +300" | (3) 


and that the percentage ac-de difference of a thermal 
converter having the hot junction of a thermocouple 
fastened to the midpoint of such a heater is approxi- 
mately 
6= —25H6,.° (4) 
where 
6.-—the midpoint temperature rise in de- 
grees centigrade on alternating cur- 
rent of rms value J amperes 
64-—=the midpoint temperature rise in de- 
grees centigrade for the same value 
of direct current. 
h=a—8—2NTZ’, where a is temperature 
coefficient of electrical resistivity. 
8=the temperature coefficient of thermal 
conductivity of the conductor, 

N= ptK/ak. 

p and a=the perimeter and area of the cross 
section of the conductor, respec- 
tively. 

¢=the emissivity. 

K=the Stefan-Boltzmann constant in 
watts-cm~*-(deg K)~* 

k=the thermal conductivity in watts- 
em~'-deg™! 





the ambient temperature in degrees 
Kelvin 

d/l? and is (approximately) the ratio 
of the crest value of the cyclic 
fluctuation of temperature at the 
midpoint to the average temperature 
rise at the midpoint 

2xf, where f is the frequency of the 
applied current in cycles per second. 

the thermal diffusivity of the heater 
material in em?-sec™'. 

one-half the length of the conductor in 
centimeters. 

h+B/2A, where A and B are the con- 
stants in the emf-temperature char- 
acteristic for the thermocouple it- 
self, that in the expression, 
E=A6+-(B/2)8 


ese equations indicate the factors upon which 
the low-frequency error depends. They show that 
the temperature rise can correctly be used as a 
measure of the rms alternating current, at fre- 
quencies high enough so that the second term in 
the bracket of eq 3 is negligible for the accuracy 
desired; in other words, at frequencies for which 
the heater integrates reasonably well. This is true 
even if the thermal and electrical conductivities 
of the heater are dependent upon temperature. 
\s a corollary, these equations show that a thermal 
converter correctly indicates rms response, except 
it low frequencies, even though its response char- 
acteristic E=f(/) is not quadratic. 

The above equations are valid only at frequen- 

es for which the second term in the brackets of 
q 3 is small compared to 1. 

Equation 4 may be used to calculate approxi- 
mately the frequencies for which ac-de differences 
may be expected to cause an error of an arbitrarily 
chosen value. This has been done for some typical 
materials for an error of 0.1 percent, a temperature 
rise of 200° C, and a typical heater length of 0.4 em. 
lhe results are shown in table 7."* 
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Frequencies for which the computed ac-de differ- 
ence is 0.1 percent 


gnifies a thermocouple of Chromel P and Alumel, and TC-B a thermo- 





{ copper and constantan, for which B/2A=—0.1 x10 and + 0.9xi0—, 
vely} 


Frequency for 


=(.1* 


Approximate temperature 
coefficient 
_| Thermal 
diffusivity 
Thermal 
conductivity 


ter material 


Electrical 
resistivity 


%/°C cem?-sec 
—0. 2x10% 1.1 
on 
075 
. 04 
058 
> 


heater 0.4 cm long with a temperature rise of 200° C. 


nated from Lorenz’ Law 


Tect of radiation losses has been neglected in pe this table. 
this effect is significant only in the very low and very high range 


mverters 


485694—52 


It is evident that the most significant factor in 
the error term of eq 4 is the length of the heater 
as the error is dependent upon the fourth power 
of this quantity.” It is of interest to note that 
the lengths, 2/, of the heater of the 1- and 2-amp 
thermal converters are only 0.28 and 0.19 em, 
respectively, and that the length of each of the 
higher-range thermal converters is about 0.5 cm. 
The value of H can be determined from the depar- 
ture from “square law” of the characteristic equa- 
tion H=f(/) of a thermal converter. For these 
thermal converters it was considerably less than 
the factor given for platinum in table 7. The 
estimated frequencies for a 0.1-percent error, on 
the assumption that the diffusivity of the plat- 
inum-iridium alloy used in these thermal converters 
is the same as that of platinum, were about one-half 
the observed values. 


7. High-Frequency Response 


As these thermal converters are inherently suit- 
able at frequencies up to 100 Me or more, the upper- 
frequency limit is set by the circuits in which they 
are used. The precautions indicated in section 4 
are sufficient to control capacitance and leakage 
currents in tests of ammeters at frequencies con- 
siderably higher than the 20,000 c/s established 
as the present upper limit by the demand for this 
testing service. The actual upper-frequency limit 
has not been established. 

The magnitude of the effective impedance sets the 
upper frequency limit for the model A and the model 
B voltmeters. For an accuracy of 0.01 percent at 
20,000 c/s, a time constant of less than 11077 and a 
skin effect less than 0.01 percent are necessary for 
the impedance of each step of the model A and each 
range of the model B instruments. The time con- 
stant of each decade of the model A voltmeter was 
computed from the inductance and direct-capaci- 
tance values furnished by the manufacturer of the 
resistance box used. This was less than 5 107° for 
each decade. The skin effect was also negligible for 
each decade. The resistance box is connected as a 
three-terminal resistor, with the shield connected to 
the grounded side of the line. For this connection 
the errors caused by the capacitance currents to the 
shield are much less than those that would occur if 
the shield were connected to a resistor terminal. 
They are significant only for the two 1,000-ohms- 
per-step decades. The effects of the capacitance of 
the resistance cards and switch studs of these decades 
to the shield were computed from the manufac- 
turer’s capacitance figures of about 10 uuf per decade 
by considering the capacitance as distributed along a 
uniform transmission line, composed of resistance 
elements and short-circuited at the end. From the 
usual transmission line formulas, the magnitude of 


" For a given heater alloy, the low-frequency range can be extended by using 
longer heaters, with proportionately greater cross-sectional area to give the same 
temperature rise for the same current. If solid heaters are used, the improved 
low-frequency performance will be accompanied by increased error due to skin 
effect at high frequencies. 
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the current at the short-circuited end was computed 
as approximately 
v 1—0.0056 (wR) 
Re 20 lw ‘ 
where 
\'=the applied voltage 
] =the current at the short-circuited end 
?=the resistance in the circuit 
(=the total capacitance to 
windings in use 
w=2ef. 
At 20,000 c/s for a capacitance of 30 yuf this is 
changed by less than 0.005 percent at 20,000 ohms, 
the maximum setting. The effect of the capacitance 
(about 10 yuf) of each rather large switch blade to 
the shield was separately computed, and found to be 
less than 0.005 percent at 20,000 c/s. Thus it was 
concluded that the errors of this voltmeter should 
be a maximum of 0.01 percent at the highest settings 
and frequencies used, and less at lower settings and 
frequencies. The performance at the top settings 
was verified by the tests outlined in section 4. 

The impedance of the model B voltmeter was simi- 
larly computed, and only the capacitance to the shield 
was found to be significant. As a rough but ade- 
quate approximation this may be considered the 
same for each step, which leads to the approxima- 
tion of considering the capacitance as distributed 
along a nonuniform transmission line composed of 
resistance elements and short-circuited at the end. 
The effective impedance was computed by a pertur- 
bation method,” on the rough assumption that the 
resistance per unit length was proportional to the 
square of the electrical distance from the output end. 
The magnitude of the current at the short-circuited 
end was found to be approximately 


shield of the 


Flt — 0.0022 (wCR)’]. 
‘ 


(6) 


From the measured capacitance of 50 yyf and the 
total resistance of 100,000 ohms for the 750-v range, 
the change in current was computed as —0.09 per- 
cent at 20,000 c/s, which is in very good agreement 
with the observed results of table 5. Moderately 
good agreement for the other ranges was also ob- 
tained, but the change in current was not observed 
to increase in proportion to the square of the fre- 
quency as indicated by eq 6. However, the agree- 
ment 1s good enough to show that the ac-de difference 
of this instrument is due to capacitance currents to 
theshield. These currents cannot easily be reduced in 
a shielded multirange instrument of reasonable size. 
(A shield is desirable to control the currents.) They 
could be compensated by connecting appropriate 
capacitors across sections of the resistors, but this 
seems less desirable than the procedure of applying 
small corrections and has not been done in the model 
B voltmeter. Much of the capacitance is _ that 
between the binding posts of each range (see fig. 2) 


1 See appendix 4 


| 





and between the high-side binding post and th i 
shield. This shield is brought as close as 2.5 er t, 
the row of high-side posts. 


8. Other Effects 


The intercomparisons of insulated-heater ther nq] 
converters disclosed an unsuspected ac-de differ:-nee 
that was dependent upon the voltage between the 
heater and ground and upon the heater current of 
each of these elements. This was traced to a very 
small leakage current flowing through the ceramic 
bead between the heater and thermocouple and thep 
through the stray capacitance of the Lindeck poten- 
tiometer to ground. This current caused additiona| 
heating of the bead and thus an error as on direct 
current the insulation resistance of the Lindeck was 
much higher than the capacitive reactance, even at 
20 c/s. In addition, the ceramic bead of each of 
these thermal converters exhibited marked dielec- 
tric absorption, so that its a-c resistance was much 
less than its d-c resistance. For a typical thermal 
converter at rated current a leakage current of 0.7 ua 
resulted from a difference of 50 v between the heater 
and ground, resulting in a power dissipation of 35 yw 
in the bead. This was sufficient to cause a 0.2-per- 
cent error, Which was proportional to the square of 
the voltage difference and was very markedly de- 
pendent upon the heater current.” The reactance 
of a capacitance as small as 100 yyf is less than the 
computed bead resistance, even at a frequency as 
low as 50 c/s, so that in almost any application this 
source of error makes it necessary to maintain the 
heater near ground potential. 


9. Other Applications 


This equipment was designed for the testing of 


other instruments. Portability and the ease of mak- 
ing other a-c measurements were not considered to 
be important. However, the thermal converters and 
voltmeter elements are readily portable and have 
been used, with a small portable millivolt potenti- 
ometer as the indicator, for measurements at other 
locations. For measurements of current or voltage 
to 0.01 percent, it is necessary to use a transfer 
method, calibrating the thermal converter or volt- 
meter on direct current at the time of use. For a 
more moderate accuracy of about 0.1 percent, it 
should be feasible to use a low-range thermal con- 
verter, shunting it for current measurements and 
adding a series resistor for voltage measurements 
Only the thermal converter itself need be transferred 
to direct current if the resistances of the shunts and 
the series resistor are properly adjusted. Thus d-c 
sources of extended range, volt boxes, etc., should be 
unnecessary. The whole apparatus could be built 
in a convenient case, contaiming a small potenti- 
ometer for the measurements of the current through. 
or the voltage across the thermal converter heater. 
1? This was due to the marked decrease in the resistance of the ceramic bead 


with an increase in temperature. With a fixed voltage across the bead, the error 
was in some cases proportional to the tenth power of the heater current 
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. suitable low-range thermal converter should also 

-e feasible an a-c potentiometer with an inherent 
racy of 0.01 percent. The basic circuit for a 

ur form of such a potentiometer (to measure the 
rnitude but not the components or phase angle of 
alternating voltage) is shown in figure 6. With 
switch thrown to direct-current the potenti- 
eter is standardized in the usual manner, and the 
flection of the indicator of the thermal converter is 
served. The switch is then thrown to alternating 
rrent, which is adjusted for the same deflection 
this indicator. The circuit is surely not new, but 
proved thermal converter makes possible a con- 
nient transfer with accuracies comparable to those 

a d-e potentiometer. The other basic limitation 
mains that any a-c potentiometer measures only 
the fundamental component of the unknown voltage. 
However for an alternating quantity having a har- 
monic content as large as 2 percent, the rms value 
differs from the fundamental by only 0.02 percent. 
\udio-frequency oscillators and amplifiers with am- 
plitude distortion much less than this are now avail- 
able. Such an oscillator could be used to supply the 
potentiometer directly and to supply the test circuit 


through an amplifier preceded by a suitable phase- | 


shift network. 
10. Summary 


The electrothermic transfer standards used at the 
Bureau for the precise measurement of voltage and 
current at audio frequencies, and designed primarily 
for the standardization of electric instruments, have 
been described. The extensive tests to establish the 
transfer performance of the standards have been out- 
lined, and the results have been given. The causes 
of the observed ac-de differences have been outlined 
and the errors evaluated. 
the temperature rise of a conductor heated by an 
electric current and cooled by conduction to relatively 
massive terminals have been solved to establish the 
steady-state midpoint temperature rise on direct 
current with Peltier and Thomson heating, as well as 
the ordinary resistance heating, in order to obtain 
the d-c error of a thermal converter. The equations 
have been solved to establish the average midpoint 
temperature rise on alternating current, in order to 
obtain the low-frequency error of a thermal converter 


Che equations governing | 


in which the electrical and thermal conductivities of | 


the heater are temperature dependent, in which some 
loss of heat occurs by radiation, and in which the 
characteristic of the thermocouple is nonlinear. 
Possible applications and modifications of the equip- 
ment for measurements under less stringent con- 
ditions have been suggested, and the application to 
an a-c potentiometer of excellent inherent accuracy 
has been diagrammed. 

The results of this work establish the excellence of 
thermal converters as transfer standards of highest 
accuracy, comparable over wide ranges with that 
obtainable with other types of standards at much 
narrower ranges. The requirements for such stand- 

rds are feasible and the manner of use to insure such 
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accuracy not unreasonable. The cost of the required 
number of standards and the equipment is not pro- 
hibitive. Thus electrothermic instruments may be 
ranked on a par with electrodynamic and electro- 
static instruments for a-c measurements of the high- 
est accuracy. 

As a result of this work, thermal converters of 
negligible or small known transfer error are now 
available as transfer standards with which the trans- 
fer performance of other thermal converters may be 
evaluated at the NBS to 0.01 percent at audio 
frequencies. 


The autbor acknowledges the help of Murray Blitz, 
who made most of the intercomparisons described in 
section 4 of this report and checked all the tables and 
equations, and of Earl Williams, who constructed all 
the equipment and assisted with the measurements. 
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Appendix 1. Development of Working 
Equations 


By definition, the percentage ac-de difference of an instru- 
ment is the percentage difference in the quantity required 
to give the same response on alternating and direct current. 
It is 

Qac—Qae 
Q 


de 


6= 100 


where Q is the quantity (current, voltage, or power) that 
the instrument measures. 

In a “straight a-c test’? the direct current or voltage is 
adjusted to give the same response of the transfer standard 
observed on alternating current. The direct current or 
voltage is measured with the potentiometer and its accessories. 
Therefore, from the above definition, 


Qec=Qae (1 + as, (8) 





is the quantity measured with the potentiometer 
difference of the transfer standard, and Q, 
required for the observed deflection of the 


where Oy, 
4, is the acl 
is the a-c value 
teat instrument 
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Ficvure 6. Elementary form of Polar a-« potentiometer 


When a thermal converter is used for transfer tests the 
small percentage difference, 6,, between the values required 
for the same indication of the test instrument on alternating 
and direct current is determined from the resulting deflection 
of the Lindeck galvanometer 
thermal converter, the small percentage difference in output 
emf resulting from a small percentage difference in heater 
current must be evaluated, This relationship may be 
expressed as 

AE Al 
nm ’ 


BE / 


where E is the output voltage, 7 the heater current, and n 
a number that is generally to 2, but may range from 
from 1.5 to 2.2. The characteristic, n, may be evaluated on 
direct current by measuring with a potentiometer the small 
percentage differences in E resulting from small percentage 
differences in heater current successively with a po- 
tentiometer at selected current levels, then computing n 
by the above expression and plotting n versus FE. The 
characteristic need not determined with an accuracy 
better than | percent and is relatively permanent 

If the Lindeck potentiometer shown in figure 3 is used to 
measure FE and AEF in millivolts in a transfer test, then at 
balance E~ V,.% the voltage drop across the four-terminal 
resistor If the galvanometer deflection D, is directly pro- 
portional to the current through the galvanometer™ and if 
the potentiometer voltage is unchanged during a determina- 
tion, 


close 


set 


, ADR, 
AE -° 
where AD is the difference in galvanometer deflection corre- 
sponding to a difference AZ in the applied emf, R, is the 
resistance of the galvanometer circuit, and S is the current 
sensitivity of the galvanometer in centimeters per milli- 
ampere on the scale. In the Lindeck potentiometer used in 
this setup, the galvanometer sensitivity is adjusted so that 
numerically 

100K, _ 
8 


These relationships may be combined with eq 9 to give 


100 (10) 


Al 4D. 
I av, 


— 
* An accuracy of | percent is sufficient for this measurement 
® The use of a straight instead of cylindrical scale introduces an error that can 
be shown to be less than | percent if the galvanometer-to-scale distance is at 
least three times the scale length 


Therefore, for each transfer 


Therefore, for a transfer test of an ammeter with a ther 
converter that has no transfer error, the required a 
difference of the test instrument is by eq 7 and 10 


Da. 


te nV, 


by 100 (f*, =o) Da 


where the alternating and direct currents are adjusted for | 
same response of the test instrument, and the result 
galvanometer deflections are D,, and D4,., respectively 

If the transfer thermal converter has a percentage a 
difference, 5,, at the current and frequency at which the tes 
made, then the ac-de difference of the instrument under ¢ 
is, to the desired order of approximation, 


Da 


by 160 ( foe! la ) D.. F 3. 


i nv 
Similarly, for a voltmeter test, 


br - 100( Ve Ve ) Dec— Da 
V ae nV’, 

in which 6, is the percentage ac-de difference of the range 
setting) of the standard voltmeter used 

In this deflection method of using the Lindeck potentiom 
eter for rapid transfer tests, the resistance of the circuit seen 
by the galvanometer must be constant to better than | per 
cent. The resistance of the thermocouple of each therma 
converter has been adjusted to the same value so that th 
simplified formulas may always be used. The temperature 
coefficient of the resistance of the Lindeck circuit has beer 
computed as less than 0.1 percent per deg C, and computations 
show that the resistance seen by the galvanometer is changed 
by less than 0.2 percent by the adjustable battery circuit, for 
all ranges of the potentiometer except the topmost 


Appendix 2. Effect of Peltier and Thomson 
Heating 


The problem is to find the steady-state temperature rise of a 
uniform conductor of length 2/ and uniform cross-sectional 
area, a, carrying a constant current, 7, and cooled solely by 
conduction to its terminals. The joule heating of a length d 
occurs at the time rate /*pdr/a watts, where p is the electric 
resistivity in ohm-cem. The difference in the heat conducted 
per second across the two ends of this differential element is 
ak(d*6/dz*)dzx watts where k& is the thermal conductivity in 
watts-cm™'-degrees™', and @ is the temperature of the 
element in deg C. Equating these and dividing by dr gives 
the differential equation governing the temperature of the 
conductor 

ak ae (14 
dx? 
which is of the form 
a’ b, (l4a 
where b= J%p/a*k. 

If the terminals remain at equal temperatures 6, and if the 
origin is at the midpoint of the conductor, the boundary con- 
ditions are @=% at r=+l. The solution, by direct integra- 


tion and application of these conditions, is 


1 a=b (na 
and at the midpoint 
(15 





If the terminals are of different metal than the conductor 
Peltier heating and cooling (proportional to the first power o 
the current) may occur at the junctions of the conductor wit! 
the terminals. These will be equal and opposite in sign a‘ 
the two terminals, and, if we can assume complete therma 
| symmetry, their effect will be to raise the tempertaure o 


f 
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rminal by a small amount Aé@ and to lower the other by 

al amount The differential equation then has as its 

larv conditions, @=6@)+ A@ at a l/ and @—8@ Aé at 
The solution is 


b.. Aéz 
A @ I + 16) 


therefore unchanged at the midpoint.“ Of course, 
te therma! svmmetry is never attained, so that Peltier 
may cause a change in the midpoint temperature, 
ver, for massive terminals the temperature changes, 
very small, so that Peltier heating need not be con- 

d significant 
homson heating occurs in homogeneous conductors having 
perature gradient in the direction of the current flow 
he one-dimensicnal case here considered, the Thomson 
n an element of length dr having a temperature differ- 
is o/ d@ watts where a is the Thomson voltage coefficient 
eriment and theory indicate that o is proportional to the 
ite temperature, but, as the Thomson heating is small 
pared to the joule heating, we might, as a first approxi- 
on, consider o constant } quation 14 then becomes 


0 . Ip de 
eb + —" +o) 0. 17 
di a da 


is of the form 


the boundary conditions @= 6) at 2 l, where « aliak 
| steady-state solution of this eq is 


»/ é . ur 
a4 ( : ) br 
e\tanhq = sinh gq ‘ 


q l At the midpoint, where r=0 
bi? /cosh ¢ l 
6, 0 ( d )- 
q sinh q 
Expanding the hyperbolic functions gives, by neglecting 
er-order terns, 
b/ ri g 
On — Oy = 0 (1 z ) 64(1 : ) IS 
2 12 12 
sy substituting the values for 6 and ¢ and noting that 


2/pl a= V, the vo:tage drop across the conductor, eq 1 and 2 
he text can be obtained. They show that Thomson heat- 
causes a small (second order of smal!ness) decrease in the 

nperature rise at the midpoint of the conductor, and can 
ead to a transfer error in a thermal converter However, 
approximation of considering ¢ a constant is a drastic one, 
ud a more accurate approach will be to consider o= BT, 
where B is a constant,” and 7’ is the absolute temperature 

Equation 17 then becomes 


h is of the form 


T’ + wTT’ +b=0, (19 
“ the boundary conditions T=T7T) at z +/, where 
BI /ak 
This is a nonlinear equation and very difficult to solve 
Chester Snow showed that it could be reduced to a Ricatti and 


to a Bessel equation, but suggested an approximation or 
irbation method to give the desired results more simply 
| directly. The answer is not significantly different from 
nee Peltier heating depends upon the junction temperature, a second 
hange in midpoint temperature may occur 
&E/dT*, where F is the emf of a thermocouple composed of the metal in 
1 and a metal such as lead, which has very small Thomson heating 


that of eq 18. The method is powerful, however, and was 
used later in evaluating other perturbing effects, so that it 
will be outlined here. The Thomson heating, which is much 
smaller than the Joule heating, can be considered as perturbing 
the temperature distribution of the conductor. — Its effect can 
be approximated by solving the differential equation without 
the perturbing term, then substituting the result, and its 
derivatives as necessary, in the perturbing term only, and 
solving the resulting equation 
Without the Thomson heat term, eq 19 becomes 


r b 


with the boundary condition given, the solution of which is 


where the subscript 1 indicates the first approximation. 
Placing this value of 7, and 7,’ in the perturbing term of 19, 
reduces this to a linear equation 


an wher , we 
i's b 9 + wha T 4 T'»)> 


where 7',— 62/2, the solution of which is 


4 wher wh 
Pi— Teg (Pa 0 @ (Pat Toda 


At «=0 this reduces to 7T,— Ty as only first-order terms are 
given, so that the process must be repeated to find the second- 
order term indicated by eq 18. The algebra rapidly becomes 
very lengthy and tedious, but the procedure is straightforward, 
and the resulting temperature rise at the midpoint is found 
to be approximately 


T To=Tall—3.4X 10-22 10.5 X 10> *z9z4g— 8.3 X 10-225], 


where z= wT)l, and zg =wT4l 


Appendix 3. Low-Frequency Error of a 
Thermoelement 


The problem is to find the steady-state average (over an 
integral number of cycles) temperature rise at the center of a 
uniform conductor of length, 2/, and cross-sectional area, a, 
carrying a sinusoidal current and cooled chiefly by conduction 
to massive terminals that are at fixed and equal temperatures 
Assume that the current remains sinusoidal and that only 
Joule heating of the conductor occurs. If all the cooling is by 
conduction to the terminals, the differential equation govern- 
ing the temperature of the conductor is, under these assump- 
tions, 


078 . wp oO 
ak t : ams : (21) 
or a ot 
with the boundary condition 6(+ Lf 0, where 


é=temperature in deg C 

i? =the square of the instantaneous value of the current = 
P?(1—cos 2) 

rms value of the current 

w= 2xf, and f is the frequency of the current in cycles 
per second 

i=time in seconds 

m==density of the conductor in grams-cm 

s=thermal capacity of the conductor in joules-grams~! 
-deg C-! 


135 





The procedure will be to consider p and k constant and to 
solve for the periodic steady-state temperature rise, then to 
approximate p and k as p= po (1 + a@) and k= ke (1 + 88), to tm 
troduce radiation loss, and to solve the resulting equation by 
the perturbation method just outlined For the 
sidered a@ and 8@ are considerably less than 1, and radiation 
than the Joule heating of the conductor 


cases con 


With a and 8-0, eq 21 is the form 


a ” 
q h(i 1—cos 2wl 
) i 
where g = ma/h ancl / T°? po! ath 
let @~ ulrd viz Then 
“a yeu d*y 
dy da du 
Let yy Or h Then ¢ (ha 2 Cy 
Applying boundary conditions gives 
OC +L (+1049) <0 
alt r LO +9 i 0 


We may take ¢(/ ¥ l 0 
Then ¥ b/20P i and eq 22 becomes 


Oru Ou 
; gq b cos Qwt 
dy ~ Of 
with u it 0 
From complex number theory if u is 
can be represented by the real part of Ue 


harmonic in time, it 
*' where U isa 
complex function of the real variable z, and j«y-—! Simi- 
larly, 6 cos 2wt can be represented by the real part of be’ 
where b is a real number. Substituitng these quantities and 
the appropriate derivatives in eq 23 and dividing by ¢ 


results in 


wt 


jlwt 


t | 
‘ j2wgl b 24 
ds 
with U«O when « +/. The solution of this is 
1 hb [ , cosh =| 25) 
rm cosh pl 
where w= yj2eg At the midpoint 


h h 
rr , , 1 ] 7 ) 
“ cosh pl . 2wg 


as at the frequencies we are interested in cosh yl -l. There- 
fore, u-real part of ['¢"*'= real part of 

| jb . j b 

= (cos Jw sin 2Qwt ) sin 2wl 

| 2wg / | 2wg 

So that, to this approximation, at the midpoint 
br / l , , 
Omni 9 ( l ol sin 2wl ) 04(1—q sin 2wl (26) 

where 0467/2 and q=1/wgl*. Note that q is the crest value 


of the cyclic temperature fluctuation divided by the average 
temperature rise at midpoint 


In general, (<= 1),+ 7!'5, where /’, and 7’, are functions 


of r, so that ul’, cos 2wt—T', sin Qwt 
and 


b 
, 


6, (i r U. cos Qwl { » Sin Swi 


te 
~~ 
~~ 








value over an integral number of evek 


Thus the average 
the temperature rise is the same as that found in append 
for the same numerical value of direct current, so that uw 
temperature rise of the conductor 
current without ac-de er 


these conditions the 

be used as a measure of the 
We shall now consider that p= pe(l+a@) and k=—ky(l 

and that some small fraction of the total heat is lost 

radiation If pis the perimeter of the heater in centimet 

then the heat radiated per from the surface \ 

differential length, d,, of the heater is 


second 


dil pt K(T*— T))dz, 

where ¢ is the emissivity of the heater material, a numer 
constant, and A is the Stefan-Boltzmann constant in wat 
em™~*(deg K)~*. TJ’ is the absolute temperature of the 
ment; 7) is the absolute temperature of the surroundi 


Introducing these terms in eq 22 gives, as fé<1, 
) 0 Kh, 
.—g( 1— 0) b( 1 —cos Qut)(1+90) +78" (Ts 7 
a yt ah 
where »=a@—8. The perturbing terms are 
Ni 
1) 3g0 
7 ry; 
2 bol cos 2wi nf 
nitK 
3) ee (T—T? 
ak 


As these perturbing terms are relatively small, the solution 
will differ by only a small amount from that given by eq 27 
Hence a sufficient expression for these small terms will be 
substituting the known values of the dependent 
variables y from eq 27, and its derivatives as necessary 
in these perturbing terms. When this is done, it is found 
that the first of these terms reduces to functions of z, each 
multiplied by sin nwt or cos nol, where n=2, 4. The second 
reduces to similar harmonic terms plus the terms —b[,/?’ 
nl’,/2)|, where P<(P—2*)b/2. The third term becomes 
unnecessarily complex It can be simplified by making the 
further approximation that for this term 


given by 


T— Ty 0,— Pq sin 2wt 
and 

T*—T $= (To+0,)*—T§ = 40,73 +-607T), 0-<T 
and also that k=ko. With these approximations, this term 
reduces to harmonic terms plus the terms 


INTIP +6NTIP( 14%), 


where N = pt K/ako. 


Thus eq 28 becomes 


ot — 9 = —0( 1 4+-9P vs) { 2NTiP| 27, 3P 14 ($) | 


oY ad ol 


ne 
‘ . 
ot F,, cos nwt + 


G,, Sin net ), 29 


where F,, and G, are known functions of x alone. We again let 
Or, Dm=ulz, O+¥(2) and set O%y/dz* equal to the terms 


independent of ¢ in eq 29, so that 


" ne 
Oru Ou ’ , 
- STF n COS nol +G, sin not) 


or ot 


(30) 


ne? 
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ach frequency of the harmonic driving terms in the 
ind side of this equation, the solution of the equation 
unishes at the end points will be of the form Ul, 
sin nowt, where U/, and Ul’, are functions of z. The 
value, over an integral number of cycles of the 
frequency, will therefore vanish Since this is a 
juation, the solution is the sum of such solutions 
ill therefore also vanish Thus the average value 
ich we are interested in, is equal to ¥ (2 
iluation ~ (2) in the same manner as before gives for 


ilue of ¥ at the midpoint where z— 0, 
q 
») | 
31 


cos 


the assumptions that [ and Uo(0) <0. 
required expression for ’, can be evaluated from eq 25 
| 27 by noting that 
ut ‘ j2wg I yl jy and “M ml JY a 
’ yw, so that cosh pz cosh yx cos yx 
ind cosh pl j sinh yi sin yl. 
cies we are interested in pl 5, so that 


+j sinh yz sin 


At the fre- 


cosh yl cosyl 


cosh ul 5 (cos yl rj sin yl). 
Sines real part of 
t cosh ur , ‘ 
Vi re 2 1 - £ e -1| U, cos 2wi—U, sin et, 
|2wg cosh yl | 
will be found that 
b 
U. {sinh yz sin yr cos yl—cosh yr cos yz sin y/]. 
wage’ 


Performing the desired integration results in 


b 
U. cos yl cosh yz cos yz 


z sin yl sinh yz sin yr). 
2way'e7 


From which 


lal —t) and U,(0)< <U.(). 
g 


rherefore, eq 31 becomes 
(0) =0..=0.1 14¢4+(7%—1 NTw# )o g 
dis pete’ @ 
h 
aa[ | + M4 eaa? | (32) 


where M is equal to the remaining terms in eq 31 and h=y 
2N TRE. 

By the same procedure, eq 28 can be solved for the d-c case 
by considering 00/0t=0 and iJ, a constant of the same 
numerical value as the rms value of the alternating current. 
rhe temperature rise at the midpoint will be found to be 


Bac O41 4+- M]. (33) 
rhus the temperature rise on direct current is the same as on 
alternating current when the frequency is high enough so that 
the third term in the brackets of eq 32 is sufficiently small 

pared to unity. Therefore the temperature rise at the 
lpoint on alternating current is 


hO aq? a ‘ 
B., eae 1 tar =o ] =e [1 n eat | (34) 


s small compared to unity, 


If the hot junction of the thermocouple is fastened to the 
midpoint of the heater, the emf of the couple may be repre- 
sented by 


_B 


where 6,, is the temperature of the hot junction in deg C, and 
A and B are constants to be evaluated at the temperature of 
the cold junction. On direct current 6,64. and 


y 


E = Ae, [1 j = 040} 


On alternating current of the same rms value, 6, mav be 


m 


considered to be closely enough for these purposes 


h 
Om =O4 [ rs 0a | 1—q sin 2wt 


The average value of emf over an integral number of cycles of 
alternating current is then (to the desired accuracy, and if 
(B/2.A) Og. is small compared to unity 


where HT =h+B/2A. 

For table 7, it is desired to compute the frequency at which 
the ac-de difference of a thermal converter due to this effect 
is equal to 0.1 percent. From eq 7 and 9 and for n=2 the 
difference in current for the same emf is one-half the difference 
in emf for the same current, and is opposite in sign. There- 
fore the percentage ac-de difference, 5, is approximately 


25110 4d? 
é 25Hé,.@ : , 
i (Qaf) 2" 
where in addition to the symbols already defined, d= 1/g= 
ko/ms and is the thermal diffusivity in appropriate units. 


Appendix 4. Effect of Shield Capacitance 
on the Model B Voltmeter 


We may approximate the effects’ of capacitance to the 
shield of this voltmeter by considering the voltmeter as 
consisting of a transmission line of r,=—rer? ohms per unit 
length, having a constant capacitance cy farads per unit 
length to the shield, witha sinusoidal voltage of rms value V 
volts applied to the high-resistance end. The other end may 
be considered short-circuited, since the thermal converter, 
which is connected between the low-resistance end and the 
ground post, is of low resistance. The shield is also connected 
to this ground post. We wish to find the magnitude of 
current, Jo, at the short-circuited end For an element of 
length dr at a distance z from this end, we have 


dV 


ete (35 
dz J ' ’ ) 


Ir, = Trex? and al 


iz 


where d\ is the rms voltage drop across the element, and 
dl is the rms current from the element to the shield. Differ- 
entiating the first of these and substituting the values of J 
and d//dz in the resulting expression gives 


@V 2dV : 
Mr? 0 (36) 
dx* x dz wart , . 


subject to the boundary conditions V=0 at r=0 and V=V, 
at z=/ and where \? = jwcoro. 
When A=0 this becomes, with the 


conditions, 


same boundary 


fy 2dV 
0 
dx* x dr , 
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the 


solution of 


whi 


We now put this value of 


Me i eq 36, to get 


in the perturbing tern 


contaming 


a\ 2» dV VV 
dz r dr l 
The solution of this, with the same boundary conditions, is 
V2’ el bg Nea! 
J (1 = )4 a 
\s A 


jut we 
Putting the above 


value of |} 
solving as before 


“ives 





will need to re peat this proces 
terms in A‘, for we are interested in the 


to secure 


and / 


last term of eq 36 and 


magnitudes of | 


in the 
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From eq 35 


F 1 dl 
ror? dz 
Performing the indicated differentiation of eq 37 and 
stituting in 38, and noting that the total resista: 
R <= rel*/3 leads to 


Va n2/4 NArt x24 Mfrs 
/ ! +88 10°>*A24 (1 )+ 4 

R 2s : 28 672 
Solving for I, at 2 0 


is / 


, and noting that 
Ce ZIVes 


the total capacit 


J eL JOLOTWC R —7.9 10 Rr) | 


Since the terms in the bracket are small compared with 


£ 
99 » 
I AL 0.0022 (wo K ] 


Wasuineton, August 16, 1951. 

















»me Nonnegative Trigonometric Polynomials Connected 


With a Problem in Probability 


Eugene Lukacs! and Otto Szasz’ 


Let 0b <b b, be integers and let g(@) be the Vandermonde determinant 
formed from 6}, t , 6 with the first row replaced by sin? (b,0/2 1, n The 
function g(@) is then a cosine potynomial. In connection with a problem in probability, the 
question arose as to when g(@) is a nonnegative trigonometric polynomial. This question is 


answered for six classes of such trigonometric polynomials 


a previous paper * we gave a necessary condi- | the generalized Vandermonde determinant formed 
that a polynomial without multiple roots must | from the numbers 47, 63, , & with the exponents 
<fy in order that its reciprocal be the Fourier | k, 1, 2, . , (n—1). Clearly V, ' ne 
sform of a distribution function. Imposing a | V,_,—0, while V, is the ordinary Vandermonde 
ther restriction on the polynomials it is easy to | determinant of the numbers 43, . _ &. Denote 
ive the following condition further by A,, the minor of the element in the first 
lhe reciprocal of a polynomial whose roots are ;, , 
. : . : “OW : m . ls at | > “4 ym 
single and have the same imaginary part is the | °™ and m-th column of |) so that 1; a A. 
Fourler transform of a distribution function if and | [pn this paper we study the function 
if 
The polynomial has one purely imaginary root be b.0 b.0 
a0) and n pairs of complex roots + &-+-ai sin? oy sin? ——» +++ om* — 
h b . b,., k e are = i +0 
2) The determinant he Mary be 
5,0 b.0 b,0 q(0)—2 
sin“ ——» a > * ’ sin* —5 ‘ ” 
b b2, o9 b? 
> for all @. be" a EE b? 
Vo+ >>| 1)"A,, cos 5,8 (2) 
> > m=1 
h 4 b3 le a b? a. 


This condition follows easily from formula (4.4) of | For the discussion of g(@) we need the following 
the reference cited in footnote 3. It is therefore | lemmas: 

of some interest to study this determinant and to Lemma l. If b, dy, , 6, are integers, the 
investigate for what values of b,. by, er * b,, it is determinant (2) can be factored so that 

a nonnegative function of @. 

In this paper we consider this determinant only 
for integer values of the 6; and show that it repre- 
sents for certain configurations of the bya nonnegative | where A(x) is a polynomial in « of degree b,—n 
trigonometric polynomial. Certain relations for | Jemma 2. If b, bs... .. b, are odd integers we 
generalized Vandermonde determinants of odd in- | paye 
tegers are also obtained. 

We introduce first some _ notations. Let g’ (8) =sin @ (1 —cos*6)"~! B(cos 0) = (sin 6)?"~'B(cos 8), 


q(@) ] cos 6)" Al(ecos @), (:3) 


Ob <b . . . <b, be n integers and 4) 
| {d 
bP, OP} sass b | 
é ; | where B(x) is a polynomial in « of degree 6,—2n-+1. 
b?, ae b? 
F | 
| (1) | To prove lemma | we differentiate (2) with respect 
to @and then set @=0. This shows that 
| g@-'(0)=0 for j=1, 2, . . . ad inf 
gae-) ite 62 (*-» | 
aS st | 94 9 { ) x 
vtional Bureau of Standards and Our Lady of Cincinnati College | q (0) 0 for) 0, I, =? n I (2) 


versity of Cincinnati and National Bureau of Standards 
Lukaecs and O. Szfasz, Certain Fourier transforms of distributions, Cana- | k k , . 
Math. 3, 140 to 144 (1951). | g™(0)=(—1)'V, fork=—n, (n+1),...adinf 
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From (5) we obtain the expansion of g(@) into a series 


qg(@) 6) 


On the other hand, (2) indicates that g(@) is a cosine 
wlynomial of degree 6, which vanishes for @—0 
Therefore it has the form 

q(0) cos 8) cos 6 


pI l B,| l 


B, (1—cos 6)’ (7 
n 

Comparing (6) and (7) it is seen that B,=— £, 

8,,=0 so that (1—cos 6)" is a factor of g(@); 

establishes lemma 1. 
To prove lemma 2 

2, ...,n) are odd numbers. 

(2) that 


this 


we assume that all the 4,(i=1, 
It follows then from 


q(x—0)=2V,—q(0), 8) 
therefore 
q' (w—0)—q' (0). 9) 
From (3) we have 
q’ (0) (1 cos 6)"~' sin O(n. A(cos 0) 
(1 COs 8)A’ (cos 6)] 


Substituting this into (9) we obtain 


(] cos 6)" sin A\n. A(cos A) (] cos 80)A’ (cos @)] 
(| cos @)""' sin An Al cos @) 
(1+- cos @).A’(—cos @)}. 

This shows that nA(cos @)—(1— cos @)A’(cos @) and 


therefore also g’(@) has the factor (1+ cos 6)" This 
completes the proof of lemma 2, 
If we introduce as a new variable 
zr cos 6 (10) 
and write 
P(x) gare cos 2) (11) 
we have from lemma 1 
P(2)=(1—r)"A(z), (12) 
and similarly from lemma 2 
P’ (x) (1—2*)*"'B(z). (13) 


We substitute next the expansion (6) into (8) and 
see that 


’ = Vv 2 V 
9» _ a) ya 1 a at. -. x A j i : 2k 
2V, ee, | l aR? | l jap * 6) 
(14) 














0 we obtain fro 


Differentiating (14) and setting @ 
this equation a number of relations for generaliz: 
Vandermonde determinants formed of odd integers 





9V a 1)*-! Vv, rt } 
- 0 
4 (2k)! 
>* (—1)"! Ve , r*=(0 
= (2k—m)! 
for m=1, 2,...(2n—1) 
a ry 
Se i 1)*-! Vir ‘ 
» St (2k—2n—2m)! 
{ Ltt g2ie tm mass 
for m=0, 1, 2,... ad inf. 
. — 
we? 
Ss. (—1)" Ul 0 
hottm+i (2k—2n—2m—1)! 
for m=0, 1, 2,... ad inf. J 


In the following we discuss several configurations 


of the integers 6, b), . . ., 6,, which lead to non 
negative trigonometric polynomials g(@), the results 
are given in statements, labelled (A) (B) .. . (F 


(A) If the 4, are the first n consecutive integers, that 
is if b,=2 for 1=1, 2, ., n, then the trigonometric 
polynomial g(@) is nonnegative. 

Proof: From lemma 1 we see q(@) 

where A is a constant. Hence A=2 
n 

Vo+S5A,,. Therefore A>0 
m= 1 


and consequently g(@) nonnegative for all values of ¢ 


A(l 


"g(a); from (2 


COs A)" 


it is seen that (1) 


(B) If the 6, are the first n consecutive odd integers, 
that is, if 6,27 —1 for i=1, 2, _ n then the 
trigonometric polynomial g(@) is nonnegative. 
Proof: Since b,—2n—1 the polynomial B(x) of 
lemma 2 reduces to aconstant B. From (13) we have 
P’ (2x) B(i —z*)*"' or, in view of P(i)=0, 
1 


P(xz)=B ' dt. 


Jt 


(1 —f#°)* 


Formulae (11) and (2) show then that ?(0) 
Vo. 0, hence 


g\r 2) 


Therefore P(x) >0 for |z|<1. This completes the 
proof of statement (B) since the inequalities 
P(x) >0 for |x) <1 and g(@) >0 for all @ are equivalent 


(C) If 6,—n+1, that is if the numbers 4, . . ., 4, 
are obtained from the first consecutive (n+ 1) integers 
by omitting the integer k (1 <k<n) then the trigo- 
nometric polynomial g(@) is nonnegative if and only 
if 2k? >(n+1). 

Proof: In this case P(x) has degree (n+ 1) so that 
according to lemma 1 A(z) is a linear function. 


(16) 


A(x) =a-+ br. 
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substitution (10) transforms cos ké into 7;(x) 
T.(r) is the k-th Tchebicheff polynomial of 


rst kind Formulae (2), (3) and (16) show that 


] r)"(a+ br) V >» 1)"4,,2'9. (z). (17) 
coethicients a and b of A(x) can be determined 


, after some elementary computations it is 


] | 
1)*—k? J 


A(x) and P(a 


17 
that 


_ wih 
Ala 2°A4.1 2-4 
4 


the 


1S 


he interval |a functions 


r)" A(x) have the same sign so that it is sufficient 

etermine when A(z) is nonnegative in |z|<1. 

| m (18) itis seen that A(x) >0 in |r!<1 if and only 
| k?>(n+-1)* —F&, that is if 2k? >n+1. 

D) If the 6, (@=—1, 2, . .n) are all odd numbers 

if b, —=2n+1 so that b,, bo, . 6, are obtained 

n the first (n+-1) odd integers by omitting the 


l odd integer 2k —1 (1<k<n), than the trigono- 
ric polynomial g(@) is nonnegative if either one 


the other of the following conditions is satisfied 


| order 


] 

k>= (1+ y142) i 
l - 

k< 5 (1+ y1+2n) il) 

Mt Qkhik l 
2'n(n+1)—k(k—1 
here is the root of the equation 
1 —#*)*-"(¢ dt=0 19) 


which falls into the interval (0, 1). 

If neither (i) nor (ii) is satisfied, then the function 
g(@) assumes also negative values. Moreover it is 
possib le to simplify (ii) for large n by proving the 
Corollary to statement (D). 


If for large n and k< [1+ 11+ 2n]/2 also 
n+; )[n 2k(k—1)] 
— — <p=0.5939157... 
2[ni(n + ]) k i 1)| . ; , 
then g(@) is nonnegative; if on the other hand 


(n+5)[n—2kk 


1) v 


2[n(n 1 


en there are integers n and k for which g(@) assumes 
itive and negative values. 


To prove statement (D) we need the following 
lemma 
Lemma 8. Let 
P&a 1 —t?)"-"'t*dt for s=0, 1 (20 
and 
Q(z, =P, (a P(x) (21) 


The polynomial ((z, ¢) is nonnegative in the interval 
1<x<11f and only 


§ 22 
Here 2, is the root of the equation 
R(z)=Q 1—1)"-"2— 2 dt=0, (23) 


which falls into the interval (0, 1 


n- 
v"2 


p exists and is the root of the equation 


Moreover let Wn 


Then lim w, 


: : , l 
which is located in the interval ( On = ) 


Some values of 2, as well p were computed. 
The proof of this lemma rather lengthy. In 
to avoid interrupting the discussion of the 
various nonnegative trigonometric polynomials, this 
proof will be given in the last section of this paper. 

We proceed now to the proof of statement (D). 

According to lemma 2 2(z) is a polynomial of the 
second degree. Since the 6; are all odd numbers 
P’ (x) and B(x) are even functions of x so that one 
obtains from (13) and (2) 


Is 


P"(z 1 —z*)"*—"'(a+ cz’ y DH Ya * T, (2). 
m= 
This relation permits us to determine a and c, and it 


is seen that 


P'(x 
20/0 DA. \ } 7] Qkik 1) 
o\en es d” > 2QIn(n+1)—k(k—1)] 
If we use for brevity the notation 
¥=2**(2n+1)A, (25) 
and 
+ 26 
ae 2in(n+ 1 k k 1)] ( ) 
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we have 


(KE 


This completes the proof of statement 
follows immediately from (30). 


wh 


ae y¥(1— -—_ Certain particular cases can be discussed eas 
si / a (a) Ifk=p—1l=—nthen N(n,n,n+1)—Lin,n,n 
Pia . | / / ; dt os 0 so that (U(r) tr? 4+-4ner+2n n—?2 The d 
criminant of G(r) is then l6(n—1)(n+2) and 
negative for n->1. Therefore, Q(z) >0 for all z 
If (1) is satishied we see that ¢2 <0, so that from (27 that g(@) must be nonnegative 
P(4)<0 for (a j In this case P(r) is non (b) If k<cn while p=(n+1) then Lin, k, n+1 
increasing in the interval |< z< 1, since PCI y | and 
this means (7) >0 for 2 1 so that also g(@)>0 V j | (n-+-2)(2n+ 1) (n k*) 
, = j (7 . 7 } 
We assume next that (it) is satisfied, then ¢ 0) ine n-+ 2) he? 
Since y -O and Ps y¥U(r, Sax), We see that state 
. (Dp ’ ; j i 
mane : it ae a ee eee (ir) and its discriminant can be computed eusil 
(FE) the numbers by, y wre obtained TPM | ond it can be shown that g(@) is nonnegative if ar 
the first consecutive (n+-2) integers by omitting the only if 
two integers # and p where 1 <k- p-—Isn then ; ° 
> . (7 AG & 
g(0) is nonnegative if and only if a certain polynomial i : : 
. » 
Mir) of the second degree is nonnegative for 4 | oll 
Here 
(ec) If fh n land pon then Ninn lin OQand 
Vir) hy itn — Lin, h ple 2n n—2Z 
] | nmin 2 
{ } , ann , nh) 
tnLin,k, p Vin, k, p), 7 { Qin +] 
where 
(n+ 2) (n+ 1)? —k] (n+ 1 p"| The discriminant of Q(z) is then 
Lin, k, p) 
, ” 
l(n & h |} {(n & P| 
, ;, 29 16 . = a 
(n 2)(2n 3)(n h(n p . on' v7 Th Tn 18), 
Vin, k, p) O(n-+- 1)? 
I(n 2? h | {in 2 PP} 
which is negative ifn >8 so that in this case g(@) >0 
Proof: In this case b,—n so that according to | If n=2, Q(z) has two roots in the taterval 
lemma } A(r) is a polynomial of the second degree, I<r< 1 so that g(@) assumes also negative 
A(r)-ar+br+e. From (12), (2) and(3) we have | values 
We proceed to discuss a more complicated cas 
: by assuming that the 6,, bd», b, are obtained 
Pir) r)"(aa bate e* : ~ . 
A from the first (n+ 2) odd integers by omitting two 
Vat S0(—1)"4,, 7's (2 (30) | odd integers (24 —1) and (2p —1) where 


From this relation we can determine the coeflicients 
a, b, and ec. A’ somewhat tedious 
computation shows that 


elementary 


so that b, 2n+3 
We see then from lemma 


nomial of degree four. 


that B(x) is a poly 
By the procedure employed 


A(z) = (ar? + br+-e) = 2"-'A, Q(z) in case (C) it is possible to determine the coeffi 
cients of this polynomial 
We obtain 
P(x) y(l—a Bix) y(1 —2*)"~"(42'+ ar?—b) (32 
Here 
y= 2°*(2n+3)4, >0, 
P 2 2(n + 1)7(38n4+-5)+(n 1)(2n+- 5)[k(k—1)4 p(p—1)) 4kik l)p( p—1) 
(n+ 2)(n 1)—kik 1)][(n+-2)(n+ 1) p(p 1)] (33 
h n+ 1)(3n +4) + 2(n + 1) [A(R — 1) + p( p—1)] —4h&(k — I) p(p—1) 
[(n-+ 2)(n + 1) — kA — 1)][(n + 2)(n+ 1) — p( p—1)] ) 


l 
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trigonometric polynomial g(@) is nonnegative if | We next derive the condition for the validity of 


only if (37) when a0. Clearly h’(r) =8r +a, h’’(1r) =8, so 
that A(r) has a minimum at +r a/s. 
*) We consider first the case {/a!/S8 a/S<1, we 
7 l {*)" (4¢'+ at’ b)dt>0 for all x <1. have then 
34) | min A (r) h( =) Be ase. 
0<r<1 Ss 16 


(he nonnegativity of the polynomial ?(2) depends 
the nature of the roots of the polynomial 

tr'+ar*—b. This introduces the need to 
tinguish a number of cases. For this purpose | ™ 


Relation (37) is therefore satisfied if —S<a<0 and 
a*4+-166<0. If !a|/8>1 we see that h’(r) <8(r —1) <0 


O<r<1, so that min A(r)—h(1)—44+a~—6. 
0o<T<l 


The function A(r) is nonnegative if a<—8 and 
t{+a—b>0. Buta®?+16b—(a+-8)? —16(4+a~—6) for 
any a, hence if a?+- 166 <0 we always have 4 +a—b>0., 
7 The alternatives discussed for case a0 can there- 


write 


tr°+ar b t(7 Tir ale 








a—ya*+16b | a+ ya*-4 ad *’’ | fore be jomed so that condition (8) is established. 
8 ’ 8 | To complete the discussion we must consider the 
remaining cases, namely 
We are now in a position to formulate statement (F) ’ 
| If the numbers 6,, b:, . . ., 6, are obtained | (VY) a>0 and 6>0, 
om the first consecutive (n+-2) odd integers by | and 
tting the two integers 2k —1 and 2p—1 (where | (§’) a<0 and a?+166>0. 
k<p—1i<n) then the trigonometric polyno- 
mial g(@) is nonnegative if and only if one of the | In case (y’) we see from (35) that 7,<0 while 7.>0. 
Jlowing four, mutually exclusive, conditions is | _. “1 
tisfied Since P (4) =4y (1 —#?)"-'(@ — 7) (P — r2) dt, it fol- 
LLISTICS r 
) a20 and b: 0, lows that P(x) <0 if r.>1, that is if ya*+ 166 >a-+8, 
8) a<0 and a?+ 166 <0, but this occurs if and only if b>4-+-a. 
a>O and 0<b<4+a and also Thus P(x) <0 for l<r<1 if a>O0 and 6>4-+-a, 
, RV 72) —rRy(y rs) > 0. therefore (y’) reduces to the case 
a’ 
><a ft and ——-<b<a+4 and also (y’’) a>0 and 0<b< 4a. 
16 “ 
: tr.) Rly rt , if r,> : 
Rady 2 By 2) apy if 0 | In this case both roots 7, and 7. are real and a 
d Ky —y 11) — teh (—y 11) 20 | simple computation shows that min P(r) =P( —y7r.), 
~I<7sl 
” Her hence P(x) >0 in |z| <1 if and only if 
*1 
R (1 —t*)"*~"t*™ t?@—2*)dtfori=1,2, (36) *1 J 
‘ P(—wvVr2)=47 (1—t?)*"'@—7) (@#—1,) dt>0. 
. VT2 
so that Ri(2)=R(z) as given by (23). That is if 
| : ’ - 
d . ’) wor If A(r) Ar + ar 4 al for O< ‘> ® then (] £2)" -1 (92 r.) dt 
re y(1 —2*)"“"h(2*) <0 for l<z<l. Since | J-v7, 
P(1)=0 we see that P(x)>0 for l<r<+l] *) 
sY ¢ | rT (1 2yn—1 (¢2 x “ 
Clearly A(+) (27+ ) ———}, therefore , -- t?)""'(@— 7.) dt>0 
; 4 16 J-—vr2 
, h(r)>0 in 0<+r<1 if and only if (37) | OT, Written in a more concise notation 
a\?_ a?+-166 ; “ RAy ra) — Ry ro) >0, 
(27+ ) > = nmO<r<l. (38) 
4 16 . } 


so that (y) is established. 

We proceed with the discussion of the last case 
(6’)a<0 and a’?+-166>0 and show first that always 
a+4>0. From (33) we obtain easily 


If a>0, (38) holds if and only if 6<0, but this is 
exactly condition (a). 


(a+4) 2(n+1)*(n?+n—1)+(n-4 1) [A(kK—1) + p(p—1)]—2k(k— 1) p(p—1) 
2 [(n + 2) (n+ 1)—k(kK—1)] [(n+ 2) (n+ 1)— p(p—1)] 
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On account of (51) this fraction has always a pos 
itive denominator so that its sign is determined by 


the numerator which we denote by 


We have then 


of 


’ r 
the sign 


Min,kyp) 


M 7] kp Pp 1) | | i 1) Dhih l | 


hh 1) 


Ll) itn ? 


1)*(n?-+-n—1) 


(n l)p(p—1)+ 2(n 


If 


1) 2h(k-1)>0 then M(n,k,p)>0 on 
the other hand (n+1)—2k(k—1)<0 then = also 
(n-+1) —2p(p—1)<0. Therefore M(n,k,p) decreases 
for fixed k and increasing p as well as for fixed p and 
increasing k. Hence M(n,k,p) attains its minimum 


value if k and p are as large as possible, that is 


If ii 


Vin,k,p) > Min. nn+1) 


2(n+1)?(n nl) 2n*(n 


1)? —2n?(n+1)(n—2) ifn> 


-2n'? (nd 


so that 
if we 


1I)>O if n>1 
no restriction 


Therefore M(n,kyp) -On?(n- 
also a+4>0. It therefore 
write the ease (6) in the form 


Is 


(5’) 0O>a tand a’?4-16b>0 
We must first consider the case 
(5, ) Oa fandd >a+4 
We have then a?4+- 1660 so that A(r) has two real 
roots; from (35) we see that 7r,<0 while r,>1 
Therefore 
ad | 


r (?—r)dt<0 


Pin) 47] 1M’ 


1 so that (6,) does not yield nonnegative 
We finally have to in- 


for l<z: 
trigonometric polynomials 
vestigate the possibility that 


a’® 
t and——-< b<4+4+-a 


(dy) O>a 

16 
holds Again a’4-166>0 and A(r) has two real 
roots r, and rr, It follows then from (35) that 
1>rm>1>-—!1 


Assume first that r,>0. 
By a simple computation it is seen that (7) has 


the two minima-+ yr, and— vt, inside the interval 
(—I,41). P(e) is’ therefore nonnegative for 
1<2<1 if and only if P( yr) >0 and P(—y7,) > 90, 


According to (34) and (35) this means that 


*1 
r,dt— 


(1—t)*- "ee 


. as 








~ | 
rpdt>0 | 





(1 —t)"~ 420? ret 
vri 
*1 
T? (1—t*)"*"'(f®#—7,)dt 

. vT1 

or, using the notation of (36) 
Ri(V7r2)— 7, Rv) > 0 

and simultaneously R,(—y7))—7,R\(— yr) >0. 


If +, <0 then there is only one mmimum for Pt, 
and therefore only the first condition remai 
This establishes condition (6) and completes thy: 
proof of statement (F). 

In this section we give proof of lemma 3, which 
was stated on page 141. We write 


Par): (1 —t?)"~"t* dt (s=0, 1) 49) 


se 


and want to determine the conditions which ¢ has 
to satisfy in order that the polynomial 
On t)=P(r)— PP (2) (40 
l<zs+1. 
l<z< +1. There- 
If for some value 


should be nonnegative for 
Clearly 0< P,(2)< Po(z) for 
fore Q(2,¢)<0 for |xz|<1 if ¢>1. 
ft} such that |[f/<1 the function Q(z,%))>0 for 
zi <1 then also Q(z,¢)>0 for |z/<1 and {¢ <|¢ 
To obtain a criterion for the nonnegativity of 
Q(z,¢) in |z/ <1 we have to determine the greatest 
possible ¢ such that |f}< 1 and Q(z,¢)>0 for jz} <! 
An elementary computation shows that the func- 
tion Q(z,¢) has exactly one minimum in the interval 
l<z: 1 which is located at z t. 
We consider next the equation R(z)=Q(—z,2)=0 
or, Written in greater detail, 
(We—2%\dt 2 P(—2z)=0. (4! 


t*)* P\(—2z) 


Clearly R(0)>0 while R(1)=P,(—1)—Po(—1) <0 so 
that the equation (41) has at least one root in (0,1) 
Moreover 

dR 


“ 
— * _g2\n-1 
= 22 ja t?)"~'dt 


so that dR/dz vanishes in the interval 0<2<1 only 
at the point z=0, therefore the equation (41) has 
exactly one root z, in the interval (0,1). 


Since min Q(z, z,)=R(z,)=—0 we see that Q(z, z, 
! 


ral 
in |z|<1 and therefore Q(z, )>0 in /|2/<! 


>0 1 
if f<2,. 
If however {>-, then 
OA—z., D=—Pi(—2z.)— f*P ol —z,) 
(f?— 22)Po(—2.)<0 
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at C(r, ©) assumes also negative values. This 
jlishes the first part of lemma 3. We still have 
cuss the asymptotic behaviour of the solution 
iation (41) 

derive first a useful inequality. 
a root of equation (41), therefore 


By definition 


we see that 


We see therefore that 


30 that 


order to derive an asymptotic estimate for w, 
we rewrite (41) by dividing the interval of integra- 


on into two parts >0) and (0,1). The integrals 


over the second interval may be easily expressed in 


erms of gamma functions A simple computation 


vields 
fb hen 
a Ma sad ' 2) 
2 r (n+) 2n+ 1 
*1 
—2'} (1 {*)" 1—tdt 
write 
2(n+5)r(n+5) 
} - _ Ri 
Tin) 
*1 
I,(z) | (1 —27t2)*-"(1 —#*) dt, 44) 





we have 


20 (n+5) 


~ (n+5) 2", 


I'(n) - 


l 
Liz) \ 71 2n+1 


” 


From the asymptotic formula for log ['(x) it is seen 
that 


r (n+) 2 

COP oo -dte po 

, ») 
» 


hence 


Let 0<r<l 
'* therefore 


We next transform the integral /, 
then (1—z)e7<1 and (1l—z)""'<e 


0<.1—[(1 —z)e7|"~'< (n—1)[1—(1 —z)e’}. 


Since e*7>1+2 we have also (l—z)e7 >1—z’ and 


1—(1—z)e*<_2’* so that finally 0<1—(1—2)"~'e™~* 
(n—1)z? or 0< e~ *~9*—(1 —2z)*-'< (n—1)z*. We 
replace here z by 2*t? and see that 0<e-(@-)2#e—- 
(1 — 2*t?)"-'< (n—1)2't' for 2*#?< 1, so that 
1 — 2747)" ~'!' = e-(a-1) 0? —(n — 1)0, 2 t* for 274t?< 1 


with 0<6,< 1. 


We substitute this into the expression (44) for J, 
and obtain 


ta 25, 
l, e-(n—D (1 —2)dt——" (n—1)z* (46) 
0 oe 
for 2°<.1 with 0<6,<1 


We hav ec 


O< ¢ n—1)27¢? é 


a 3 22 
» — ad | 
= 2°*t we see that ( l—e 2 ) 


Since 0< |] é 


72??? with 0< »<.3/2 and hence 


P n—1)2*t* é at pee ne n—t t 
therefore 
*! *) (n+) ’ 
é 1 1—f?t*)dt é 1 —t*)dt 

2X 

— ¢ n a 
i 
with 0<¢< 1 and 0<X ~. 
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Substituting this into (46) we have 


(w 2 . 
/, ‘ l 1 \dt- 
»r 26,(" | - 
,- ‘ ! ’ 17 
Ld 77) 
with 0<6,<1, 0 a<1, OC Az 
| 
We introduce next a new variable w n+, 
using (45) and (47) we see that w, satisfies the 
equation 
l "I 
vr(; u )—-2u e? | f 1—f\dt 
2X Ww = be °8 w* n | 
lu em ~ > ‘ - = 0 
lS l tO l 
n-+ : ( n-+ ,) 





For each value of n this equation has a root w,, by 
(43) we see that the set a these roots is bounded. 
We consider any accumulation point p of this set 
and a subsequence {w,,\ of the sequence {w, | such 


that lim Wy =P From (48) it is seen that p satisfies 


the equation 
l ' 
F(p) v5 p*)- 2 p* ee (1 t*\dt 0 $9) 


and by transforming the integral we obtain finally 
for p the equation 


Fp) )- | reve "da o (50) 


l 
v"(5—' 


To complete the proof of lemma 3 we have only to 
show that the sequence |w,} converges. We dem- 
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onstrate this by showing that 


Ww, 


F(p)—0 has only one root in (0. ») 


We see easily that F(0)>0 and F( 


l 


» 


has only « 
accumulation point, and this follows if we show t! 


). 0 


that there is at least one root in the interval, 


From (50) we see that 


F'(p) 2o| vs +2 ( de | 


so that F’(p)<0 for p>0 
has exactly one root. 


The ~ _e (41) has been evaluated for certai: 
n, similarly the asymptotic solution p has 


values o 
been found from (50). 


Standards, Washington, D. C., 


values of w, 


2 38197 
3 32114 
' 28244 
5 25505 
6 23434 
7 21708 
S 20463 
” 19346 
10 S305 
15 15122 
20 13141 


p =0.5939157 


Wasnineton, May 25, 1951 


The results are 


60304 
60079 
5991S 
50814 
59745 


5606 
50659 
59630 
50606 
59536 


59500 





This shows that - 


This was done in part at thy 
Computation Laboratory of the National Bureau of 
and in part at the 
computation department of the Institute for Numeri 
cal Analysis in Los Angeles, Calif. 
given in the following table, which shows also the 
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Measurement of Cesium-137 and Cobalt-60 Gamma 
Radiation With a Pressure Ionization Chamber 


G. R. Grove ' 


The characteristics of the pressure ionization chamber at the National Bureau of Stand- 
ards have been studied for the gamma radiation from a cesium-137 source (0.6614 Mev) and 


from a cobalt-60 source (1.1715 and 1.3316 
these sources were measured with this chamber with an accuracy of about 2 percent. 


Mev). 


The dosage rates produced in air by 
Dosage 


rates obtained with an extrapolation chamber and a Victoreen thimble chamber agreed with 
the values measured by the pressure chamber within the limits of the experimental errors. 


l. Introduction 


\ study of the charactristics of the pressure ioni- 

ition chamber at the National Bureau of Standards 
for X-rays up to 400 kv has been discussed by Taylor, 
Singer, and Charlton [1). The results of their study 
showed that it was possible to measure accurately the 
dosage rates of high-voltage X-rays in accordance 
with the international definition of the roentgen. 
lavlor and Singer [2] also used this chamber to 
measure the ionization of the gamma radiation from 
u radium source and determined a value for the 
emission constant of radium. At that time it was 
not feasible to investigate completely the chamber 
response with quantum energies in the million-volt 
region. However, Taylor and Singer believed their 
measurements were accurate to within + 0.4 percent. 

Radioactive materials of sufficient radioactivity in- 
tensity and particularly suited for use with the pres- 
sure chamber have become available only recently. 
The study of the characteristics of the chamber with 
isotopes having gamma radiation in the energy re- 
gions of 0.6 and 1.3 Mev was undertaken to determine 
the reliability of the chamber for measuring the 
dosage rates produced by radiation energies up to 1.3 
Mev. An extrapolation chamber and a Victoreen 
thimble chamber were also tested that 
readings could be compared with that of the pressure 
chamber. 

Cesium-137 and cobalt-60 isotopes were obtained 
from the Oak Ridge National Laboratory for use 
with the chamber. The cesium isotope emits a 
gamma ray of 0.6614 Mev [3], and cobalt-60 emits 2 
quanta in cascade per disintegration with energies of 
1.1715 and 1.3316 Mev [4]. The activities of the 
cesium and cobalt samples were about 8 and 1.5 
curies, respectively. 


so 


2. Apparatus and Procedure 
2.1. Pressure Chamber 


The diaphragm system and chamber are shown 
chematically in figure 1. The source is housed in an 
luminum eylindrical spool, which can be moved 
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axially inside a 1-in. diameter aluminum tube to 
vary the distance between the diaphragm and the 
source. The diameter of the aluminum tube is twice 
as large as the diameter of the diaphragm aperture 
and about three times as large as the diameters of 
the sources used. The aluminum tubing is sup- 
ported in containers that are filled with lead shot to 
reduce the radiation intensity in all directions ex- 
cept along the axis of the tube. A lead wall across 
the back of the tubing confines the radiation to a 
forward direction. 

The gamma-ray beam defined by the diaphragm 
passes through an aluminum window 0.8 mm thick 
and 5.2 cm in diameter into the chamber. The 
cylindrical pressure chamber has an inside length of 
7 ft and a diameter of 30 in. Collecting electrodes 
are arranged in the chamber on a pair of rails in 
such a way that they can be moved from the outside 
along the axis of the tank over a range of about 125 
The spacing of the electrodes also can be 

7.2cm up to about 39 em. 
These electrodes consist of aluminum plates 75 em 
long and 40 cm in height. An insulated aluminum col- 
lecting plate 25.4 em wide and 35 em high is located 
at the center of the grounded electrode. The sec- 
tions of the grounded electrode around the collecting 
plate act as guard plates. Equal potential planes are 
provided in the sensitive ionization volume by means 
0.4-mm aluminum guard wires located in planes 
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parallel to the electrodes. Each of these wires com- 
pletely surrounds the volume defined by the plates 
These guard wires, which are 2 em apart with ad- 
jacent wires connected by 2-megohm metalized re- 
sistors, produce a uniformly graded potential between 
the high-voltage electrode and the grounded guard 
plates. 

Air is pumped into the chamber through a filter 
and a calcium chloride drier, and in addition, several 
trays of drier are kept inside the pressure tank to 
keep the moisture content of the air at a minimum 
The pressure of the gas can be read directly from a 
Bourdon gage with an accuracy of about 0.2 percent, 
and the temperature of the gas can be determined 
from a thermometer fastened to the electrode assem- 
bly, which can be read through a thick glass window 
in the rear of the tank with an accuracy of about 0.4 
pereent 

The potential applied to the high-voltage electrode 
is supplied by a group of about 500 B batteries (45 
v per battery). Potentials from 500 to 20,000 vy 
are available in steps varying from about 500 to 
2,500 v. Using a 100-ya meter with appropriate 
multipliers, the plate potential can be measured 
with an accuracy of about 0.5 percent. A General 
Electric Co. FP-54 electrometer tube is used in a 
conventional detector circuit with a galvanometer of 
approximately 107° amp/mm sensitivity. A Rubi- 
con potentiometer connected to the collecting plate 
through a S.S. White resistor supplies the necessary 
current through the resistor to buck out the current 
from the collecting plate and maintain it at ground 
potential, With the grid of the electrometer tube 
also connected to the collecting electrode, the com- 
pensation current can be adjusted to exactly neu- 
tralize the ionization current as indicated by a bal 
ance of the galvanometer. 

The dosage rate supplied to the ionization chamber 
under saturation conditions can then be determined 
from a know ledge of the effective volume of air de- 
fined by the collecting plate and the diaphragm 
aperture along with the values of the grid resistor 
and the potentiometer voltage necessary to balance 
the detector. 

If the response of the pressure chamber with plate 
voltage and plate separation at various pressures is 
internally consistent and is in agreement with the 
expected characteristic response of a parallel plate 
chamber, the values measured can be relied upon as 
representing the true dosage rates. 

Using the cesium-137 source and a constant plate 
separation of 35.14 em, the ionization current was 
measured as a function of the plate voltage for several 
different pressures ranging from 4.81 to 8.90 atm 
corrected to standard temperature and pressure (see 
fig. 2). The values of the ionization currents per 
unit volume of standard air obtained by the extra- 
polation of the saturation curves to an infinite 
electric field by the technique of Jaffe and Zanstra 
[5] and making the corrections due to air absorption, 

rreed within | percent with their average value. 
With a constant pressure of 8.93 atm, the ionization 
current response to plate voltage was measured for 
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Figure 2 Voltage saturation curves for different chamber 
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Voltage saturation curves for different plate 
separations. 


Fiaure 3. 


@, 9.98 em; ©, 23.96 em; A, 35.14 em, 

different plate separations ranging from 7.18 to 37.94 
em. The saturation curves for three of the plate 
separations are shown in figures 3 and 4. With the 
smaller plate separations, the ionization current in- 
creases more rapidly in the lower plate voltage range 
than does the current from a larger plate separation, 
due to the greater field strength. When the plate 
separation is not large enough, as is the case of 9.98 
cm at 8.93 atm (fig. 3), some electrons are collected 
by the electrodes before they expend all of their recoil 
energies by ionization. From the complete family 
of saturation curves at a constant pressure, a curve 
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the ionization current versus the plate separation 

th constant electric field was determined. The 
onization current was found to fall off for plate sep- 
arations less than about 150 cm of equivalent stand- 
rd air. When corrected for the diameter of the 
amma-ray beam and the chamber pressure, a mini- 
mum sufficient plate separation of about 120 em for 
0.661-Mev radiation was found 

\ similar procedure was followed in measuring 
the radiation intensity from the cobalt-60 source. 
For these quantum energies a minimum plate 
separation of about 400 em of standard air was 
found necessary to obtain the maximum collecting 
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source-diaphragm Figure 6. 


corresponds to a plate separation of about 335 em 
of equivalent standard air when corrected for the 
beam diameter. 

The dosage rates at the limiting diaphragm surface 
for the particular cesium and cobalt sources measured 
in the manner just described were 6.91+0.07 and 
5.10+0.05 r/hr, respectively. These values need 
to be corrected only for the air and aluminum 
absorption and for the diaphragm leakage. 

A test of the inverse-square response was made 
with the pressure chamber for the two 
The distance, F, between the center of the source 
and the limiting diaphragm surface was varied, and 
the concomitant ionization current was measured, 
These data are shown in figure 5, from which it is 
seen that F' is very nearly equal to the effective 
source-diaphragm distance. 


sources, 


2.2. Extrapolation Chamber 


A paraliel-plate extrapolation chamber [6] also was 
used to measure the dosage rates of the two sources. 
This chamber consists of two Lucite disks, the inner 
surfaces of which are painted with a conducting 
aquadag dispersion. Around each disk there is a 
concentric guard ring (fig. 6). The Lucite disks 
are supported in, but insulated from, brass rings, 
one of which is threaded into the other so that the 
plate separation can be varied. A complete revo- 
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Sectional view of extrapolation chamber showing 
approximate size of gamma-ray beam. 
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Extrapolation chamber curve for the cesium-13) 


source, 


Fiounn 7. 


lution of one of the plates changes the plate separa 
tion by 0.150 em, which is about 8 percent of the 
maximum plate separation 

By the use of this chamber with a shielded con 
denser-potentiometer circuit and a vibrating-reed 
electrometer as an indicator, the ionization current 
was measured as a function of the plate separation 
A curve for the eesium source is shown in figure 7 
From the slope of the straight line determined by 


the most pomts, the average ionization current per 
This quantity, together 


centimeter can be found 
with the values of the temperature, pressure, and 
the area of the diaphragm aperture, determines a 
dosage rate at the limiting diaphragm surface 

the cesium and cobalt sources, dosage rates of 6.72 
+£0.13 and 5.20+0.10 r/hr, respectively, were cal- 
culated, 
air and Lucite absorption and for the diaphragm 
leakage 


2.3. Thimble Chamber 


A 25-r Viectoreen thimble chamber was also used 
to measure the dosage rates of the two isotopes 
This chamber and charger unit had been previously 
calibrated for radium radiation, and a scale cor- 
rection was determined with an accuracy of about 
2 percent. With this chamber the dosage rates were 
measured for the two sources at a sufficient distance 
from the diaphragm, such that the largest dimension 
of the sensitive volume of the chamber was less 
than one-half the diameter of the umbra portion of 
the radiation. With the same geometry, readings 
were taken with different thicknesses of Lucite 
shells around the sensitive volume. The rate of 
ionization in the chamber at first increased with 
shell thickness due to the approach to electronic 
equilibrium, and then decreased with additional 


These values need to be corrected only for 


shell thickness due to absorption in the Lucite I} 
extrapolating the slope of the line due to absorpti 
to zero shell thickness, a dosage rate was determin, 

This value is taken to be a good approximation 

the ionization rate with electronic equilibrium a) 
no Laicite absorption. The shell thicknesses neces 

sary for electronic equilibrium were estimated to | 
about 1.0 mm for 0.661-Mev radiation and 1.3 m 

for 1.3-Mev radiation. The extrapolated valu 
were about 8.8 percent for cesium and 4.5 percent fo 
cobalt: higher than the respective values with ele: 
tronic equilibrium shell thicknesses. Extrapolat, 
values of 7.0440.21 and 5.20+0.15 r/hr were ea! 
culated for the cesium and cobalt sources, respe: 
tively. These values are corrected for temperatur 


and pressure and are for the same source-limitin 
aperture distance used with the other two chambers 
These values can then be compared directly with 
the pressure-chamber results 


3. Corrections and Results 


The sensitive volume of the pressure ionization 
chamber is determined by the length of the collector 
plates and the area of the limiting aperture. How 
ever, there is an error involved in considering only 
the radiation that passes through this aperture, as 
the thickness of the diaphragm is not small compared 
to its distance from the source of radiation, and thus 
some of the radiation will pass through the edge of 
the diaphragm and enter the chamber. To estimate 
the magnitude of this error it was necessary to de 
termine what fraction of the total radiation produc 
ing ionization in the chamber was transmitted 
through the edge of the diaphragm. This was 
done by considering the activity of the source to be 
distributed uniformly over the surface of a thin 
disk located at the center of the source (fig. 8). If 
p ahd @ are the polar coordinates of a point on the 
surface of the disk and r is the distance from the 
source-diaphragm-chamber axis in a fixed direction 
at a distance /) along the axis from the source, where 
D is equal to or greater than the distance d between 
the source and the outer diaphragm surface, then 
the total amount of radiation received by the ioniza 
tion chamber is given by 


. prdr pdé, 


where 4 is the thickness of lead traversed by gamma 
rays emitted at p, 6, which pass through the point 
r, and yw is the linear absorption coefficient of lead 
for the radiation in question. The maximum valu 
of p is equal to the radius of the source, and the maxi 
mum value of r is taken to be that radius at which 
the radiation intensity is negligible for all values of 
pand @ The absorber thickness 4 is given by 


bed 


p)+veXp—r cos 6f + (a? — p*)(r? +p? — 2rp cos 6) 


r?+ p’—2rp cos @ 


D} p(r cos 6 
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re S Schematic indication of the variables re quir d for 
letermination of the diaphragm leakage correction. 


ill positive values of the right-hand term. When 
right-hand term becomes negative, 6 is equal to 
This triple integral was evaluated by numeri- 
ntegration, and the resulting intensity was com- 
ed with the intensity that was passed through 
inobstructed aperture for the cesium and cobalt 

, In the case of cesium the correction was 
1) percent, for cobalt it was 4.85 percent. These 
ies are believed to be accurate to about 10 percent. 
he resistance of the S. S. White resistor used in 

he grid circuit of the electrometer was measured 
during the course of the experiment and found to 
value of 3.00310" ohms at 26° C and 
a The precise temperature 


rees 


( nu 
161 10° ohms at 23°C 
of the resistor during a set of measurements was not 
known but was taken to be the average value of the 
temperature of the air in the tank and room tempera- 
The maximum error due to this resistor is 
‘ timated to be about l percent 
hy the application of the corrections due to air, 


| re 


aluminum, and Lucite absorptions [7], and that due 
to the diaphragm leakage, the correct dosage rates 
measured with the pressure and the extrapolation 


chambers were calculated. These values are tabu- 
lated in table 1, along with the thimble-chamber val- 
ues. The cobalt source had been calibrated on June 

1949. The calibrated value corrected to the ap- 
proximate date of this experiment was 5.05+40.15 r/hr 
at a distance equal to the source-diaphragm distance 
The errors assigned to the various values have been 
estimated and consist largely of the errors already 


' 
clint ussed 


Comparative values of the dosage rates at the limiting 
diaphragm surface 


ure chamber 
rapolation cham ber 
ble cham ber 


(sing the pressure chamber, the absorption coef- 
ents for the gamma rays of cesium and cobalt 
were measured for lead, copper, and tin. Sets of 
standard absorbers ranging in thickness from 0.050 
to about 3 em were used in this determination. The 
ues of the measured linear absorption coefficients 
listed in table 2, along with the best available 
oretical values [8]. The theoretical values are 
eved to be accurate within a few percent. The 
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lower energy values are the most inaccurate because 
of the error in determining the more prominent pho- 
toelectric absorption coefficient 


TABLE 2 Linear absorption coe ficients 


Measured 


Calculated 


4. Conclusion 


The response of the pressure ionization chamber 
for radiation energies below 1.3 Mev indicates that 
absolute dosage rates can be determined within a few 
percent in accordance with the international defini- 
tion of the roentgen. The accuracy of the values 
obtained with this pressure chamber are limited 
mainly by three factors: (1) the error in the measure- 
ment of the ionization current, (2) the determination 
of the saturation current obtained by a Jaffe- 
Zanstra extrapolation, and (3) the accuracy of the 
determination of the area of the diaphragm aperture, 
with the correction factor due to leakage of radiation 
through the edge of the diaphragm. A consideration 
of these factors, together with the inherent tolerances 
of the chamber, led to an estimated error of 2 percent 
for the over-all accuracy of the results. 

It appears that it is possible to get an accuracy that 
will allow dosage-rate measurements to be made with 
somewhat less than | percent error by minimizing 
the uncertainties arising from the sources listed 
above. 

The dosage rates measured with the extrapolation 
and thimble chambers are in good agreement with 
the values determined with the pressure chamber. 
This is especially gratifying in the case of the thimble 
chamber, as it is used extensively and there has been 
some question regarding energy dependence for radi- 
ation in the million-electren-volt region. The Vic- 
toreen thimble chamber was not designed to measure 
the dose resulting from high energy radiation; 
however, it was found that the dose due to collimated 
cesium-137 and cobalt-60 gamma-rays could be 
determined relatively accurately by placing Lucite 
shells around the thimble to provide electronic 
equilibrium for the radiation in question 
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Techniques for Growing and Mounting Small Single 
Crystals of Refractory Compounds 


Fred Ordway ' 


An apparatus is described for growing single crystals in a small droplet of melt at tem- 


peratures up to 1,600° ©, 
can be 


The droplet 
observed during the 


process through a binocular microscope, The 


, held by eapillarity at the junction of a thermocouple, 


heating is by 


means of a high-frequency current that can be separated by a filter from the thermocouple 


electromotive force, allowing the 


electromotive 


force to be measured continuously. The 


apparatus is therefore suitable also for determining melting points, 
The process of mounting and grinding the single-crystal specimens for X-ray diffraction 


and a simple 


1. Introduction 


A necessary preliminary 
diffraction studies is the procurement of a suitable 
single-crystal specimen. The ideal specimen should 
be pure, free from flaws or twinning, and well enough 
developed that the crystallographic axes can be found 
from the orientation of edges and faces. In spite of 
of the strictness of these requirements, usable crys- 
tals of many substances can be found in natural 
minerals, or in products of manufacture, because 
large size is not necessary. Even a crystal that is 
only 0.05 to 0.1 mm in least dimension can be used 
with ordinary X-ray cameras. Furthermore, if the 
external development is poor, the crystal can be 
oriented by X-ray patterns, 

It sometimes is found, however, that all available 
specimens of a mineral are either too impure or too 
finely crystallized and that synthetic melts of the 
proper composition, heated in a furnace, likewise 
form crystals that are too small. This is the situ- 
ation among the calcium silicates amd aluminates of 
portland cement clinker, which are being studied at 
the Bureau 

Previous experience in preparing X-ray specimens, 
from melts at low temperatures [1, 13, 18],? indicated 
that the most convenient method is manually con- 
trolled cooling of the smallest droplet that will \ ield 
a crystal large enough to be used. The small size 
of the droplet minimizes the number of nuclei and 
the accidental temperature inequalities that lead to 
growth of unwanted erystals. More important, the 
droplets can be observed by means of a stereoscopic 
microscope during the entire process of crystallization. 
Twins or other extraneous ery stals can be detected 
as soon as they form, and melted away so that growth 
of the desired crystal can be recommenced immedi- 
ately. A much larger body of melt, which would be 
necessary if the crystals were to be grown for optical 
or piezoelectric use, is actually an inconvenience when 
they need be only large enough for X-ray diffraction. 
To take advantage of this fact, a special apparatus 
was assembled for reerystallizing small samples at 


to single-crystal X-ray 


' Research Associate at the National Bureau of Standards, 
Portland Cement Association 
* Figures in brackets indi¢ate the literature references at the end of this paper 


representing the 


but versatile micromanipulator to facilitate these procedures are described, 


high under the stereoscopic micro 
scope. 

In order that the subsequent operations of mount 
ing the specimen and grinding off excessive glass o1 
crystalline material might also be carried out con- 
veniently under the microscope, a micromanipulator 
and certain accessories were constructed. The meth- 
ods described are not claimed to be unique solutions 
of the problem, but they are felt to be neers in 


certain details over the techniques seen elsewhere 


temperatures 


2. Apparatus for Crystal Growing 


2.1 Hot-Wire Apparatus 


Numerous publications have described small fur- 
naces in Which a sample can be placed [7, 15, 28, 30] 
and electric heaters on which the sample can be 
mounted directly [6, 12, 17, 21] for microscopic obser- 
vation. In the latter form of apparatus, elaborate 
insulation is not necessary to protect the microscope 
objectives if the hot zone is made small enough. 
This fact makes possible simpler construction and 
less restriction on the angle from which the sample 
can be viewed. An apparatus of the directly heated 
type was therefore constructed, be heater con- 
sisted simply of a 5-em length of 0.3-mm platinum 
wire bent to a sharp V and supporte d at the ends by 
suitable binding posts. The wire was heated by 
current from a 6.3-v transformer, the primary voltage 
of which was controlled by a continuously variable 
autotransformer or a foot-operated rheostat intended 
for use with sewing machines. The pa t of melt, 
held by capillarity at the point of the V, was ob- 
served through a stereoscopic microscope. A vari- 
able-density filter, originally intended for the eye- 
pieces of a pair of binoculars, served to protect the 
observer from excessive glare. 

This hot-wire apparatus, when protected from 
drafts, was convenient for controlled crystallization 
at temperatures up to 1,500° C. It provided no 
indication of the temperature, however, except the 
reading on an ammeter or voltmeter in the circuit 
and the behavior of the sample itself. The auto- 
transformer settings or meter readings corresponding 
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e melting range could be found by experiment, 
when the temperature was high these settings 
ially changed as the platinum wire volatilized. 
thermore, it was impossible without an inde- 
lent measure of temperature to detect small flue- 
ions due to air currents or changes in line voltage. 


2.2. Hot-Thermocouple Apparatus 
a. Construction 


, considering the addition of a device for meas- 
the sample temperature, a thermocouple was 
erred over a resistance thermometer or radiation 
meter because reasonable accuracy could be 
ined without special calibration of the sensitive 
ent, and suitable measuring instruments were 
eady available, 

\ttempts to measure the sample temperature with 

thermocouple were discouraging because it con- 

ted heat from the sample. Even when made from 
wire only 0.08 mm in diameter it greatly lowered the 
temperature. The possibility of supplying heat to 
the thermocouple from a separate source in order to 
counteract conduction from the sample was con- 
The more suitable procedure finally dis- 
covered, however, was to use high-frequency alter- 
nating current for the heating, so that the direct- 
current electromotive force produced by the thermo- 
couple could be separated out by a filter and measured 
independently. If this is done, obviously the hot- 
wire arrangement need not be a four-terminal device. 
lhe thermocouple itself can also serve as a heating 
element and support for the sample. 

The essential parts that were assembled for the 
hot-thermocouple apparatus are indicated in figure 1. 
rhe heater current is obtained from a 50-w audio 
amplifier of the type used in public address systems. 
The amplifier is fed by a Hewlett-Packard model 


sliit red 
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Circuit diagram for hot-thermocouple apparatus. 
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200BR audio oscillator through a 10-turn helical 
potentiometer, which is used as a voltage divider for 
fine control. This equipment, having been chosen 
primarily because it was readily available, is much 
more elaborate than necessary. As the output need 
be only 25 to 30 w and freedom from harmonics is 
not important, a power supply built specially for the 
purpose could be simple and inexpensive. 

The power for heating the thermocouple is fed 
from the output terminals of the amplifier through 
an additional 50-w audio output transformer, whose 
secondary is used alone as an autotransformer. To 
prevent the transformer from acting as a short cir- 
cuit for the thermocouple electromotive force, a 
blocking condenser must be inserted in series with 
the winding. The blocking condenser in use con- 
sists of three standard 10-yf “oil-filled units connected 
in parallel. 

The blocking condenser adds a considerable reac- 
tive component to the impedance of the amplifier load. 
For most efficient power transfer it is desirable to 
balance the capacitance by means of an equal in- 
ductive reactance in series with the load or reflected 
through the output transformer; that is, the load 
should be a series resonant circuit at the frequency 
used. In addition, the output transformer should of 
course have the proper turns ratio for the usual 
match of resistive impedances. These adjustments 
have not been accurately made in the present appa- 
ratus because it is found that sufficient heating can 
be obtained by the amplifier and transformer im- 
pedances shown in figure 1, with an oscillator fre- 
quency of about 4.7 ke. The resistance of the 
thermocouple is 0.1 to 0.4 ohm, depending on its 
temperature, and the reactance of the condenser is 
1.1 ohms at the operating frequency. More careful 
matching to these values would be necessary if an 
amplifier of lower power were used. If a low-fre- 
quency power supply such as a 400-cycle rotary con- 
verter or a 60-cycle transformer were to be used, the 
impedance match would be extremely important 
because of the high reactance of the condenser at 
low frequencies. The use of 60-cycle power would 
probably be rendered impractical by the large size of 
the components required. 

The thermocouple itself is mounted on an inter- 
changeable four-prong plug-in assembly. The four- 
prong plug is made by cementing together two 
standard double banana plugs which fit into a pair 
of double banana jacks fastened permanently in posi- 
tion under the stereoscopic microscope, as shown in 
figure 2. In order to replac e the thermocouple with 
a clean one for preparing a new melt, it is necessary 
merely to plug in a new assembly. 

The portion of the thermecouple that is actually 
heated, at the junction, consists of 2-cm lengths of 
0.25-mm platinum and platinum—10-percent rhodium 
wire. To reduce heating of the connectors these 
short lengths of fine wire are welded to 0.6-mm 
extension leads of the same metals, which in turn 
make contact with the banana plugs. Several of 
the interchangeable assemblies can be made up in 
advance by welding new lengths of fine thermo- 
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Fieurre 2 Thermoc ouple asse mbly 


couple wire onto the heavier leads. The leads need 
not be replaced, as they are not contaminated by 
contact with the melt; the amount of platinum used 
per sample is therefore only about 0.06 ¢ 

The thermocouple is protected from drafts by two 
sheets of mica supported by a brass frame, which 


clips onto the Bakelite body of one of the double 


banana plugs, as shown by dotted lines in figure 2. 
Without the mica cover no special protection for the 
microscope is necessary if the vertical convection 
current from the hot wire is avoided, but with it 
even this danger to the objectives is eliminated. 

The heater current reaches the thermocouple 
through the nearer pair of banana jacks shown in 
figure 2. The jacks farther from the thermocouple, 
which are less affected by the heat, can be used for 
the measuring circuit if convenient interchange of 
thermocouple assemblies is desired. For the most 
careful work two extension leads made of the thermo- 
couple metals are connected directly to the thermo- 
couple leads by means of the clamping screws in the 
second pair of banana plugs. This arrangement, 
which eliminates contact potentials between the leads 
and the banana plugs, is the one indicated in the 
figure. 

The measuring circuit is conventional except for 
the inclusion of the low-pass filter (fig. 1) to protect 
the measuring instrument from the a-c potential 
across the thermocouple, which may be as high as 3 v. 
The filter must not be interposed between the thermo- 
couple and the cold junction, of course, if the exten- 
sion leads are to serve their purpose. In order that 
a millivoltmeter pyrometer that requires a low d-c 
resistance may be used for temperature measure- 
ments, the primary of a 100-w radio power trans- 
former, having an inductance of 0.6 h and a resistance 
of 13 ohms, is employed as the filter choke. The 


filter condenser consists of two 2-uf oil-filled units j 
parallel. 

All condensers in the heating and measuring cireui 
must be of the nonelectrolytic type with essential] 
zero conductance in either direction. The slig} 
rectifying action of an electrolytic condenser cou), 
lead to a comparatively large error in the measur: 
ment of temperature. 


b. Calibration 


The convenience of the hot-thermocouple appara 
tus for its original purpose of growing small sing| 
crystals was immediately obvious when it was put 
into operation. In one of the first experiments 
lime-alumina melt in the primary phase region for 
tricalcium aluminate produced a crystal of that com 
pound better in size and perfection than had been 
found possible with any other method. The main 
value of the temperature measurement in growing 
crystals is for precise control during one experiment 
lasting up to 2 hr. Reproducibility of readings with 
duplicate samples on the same thermocouple or with 
duplicate thermocouples is less important. In view 
of the other possible uses for the apparatus, however, 
the absolute accuracy of the temperature measure- 
ments was tested as well as the precision of control 

For these tests the oscillator and amplifier were 
powered by a Sola constant-voltage transformer 
A potentiometer whose smallest scale division is 0.001 
mv and a Leeds & Northrup Speedomax recorder 
whose smallest scale division is 0.05 mv were con- 
nected so that either might be used at will. The 
recorder, which prints one point in 3 sec, follows more 
rapid fluctuations than the galvanometer used with 
the potentiometer. When the apparatus was care- 
fully protected from drafts, with the thermocouple 
at about 1,000° C, the minimum fluctuation in the 
recorded electromotive force over 10-min periods was 
£0.02 mv. As the average potentiometer readings 
agreed with those of the recorder to within less than 
0.02 mv at this and several lower temperatures, the 
thermocouple electromotive forces were read there- 
after from the recorder chart. 

Three of the interchangeable thermocouple assem- 
blies were made from single lengths of wire 0.25 mm 
in diameter. These lengths were taken from 100- 
foot rolls that had been purchased by the Chemistry 
Division of the Bureau under strict specifications for 
homogeneity. Calibration certificates were avail- 
able for couples made from each end of the original 
lengths of wire; their agreement within the limits of 
error of the calibration indicated the corrections to 
be valid for all couples made from this wire. The 
corrections amounted at most to +8.5 deg, at an elee- 
tromotive force of 15 mv, over the temperatures read 
from a standard reference table of thermocouple 
potentials [23]. 

The three thermocouples were employed, with a 
correction chart, to determine the melting points of 
five inorganic compounds. The first two couples 
were used with barium disilicate and the third with 
potassium sulfate, sodium sulfate, potassium chloride, 
and sodium chloride, in that order. The barium 





icate Was a sample synthesized as a melting-point 
idard by Eubank [8], and the other salts were of 
mnt grade. It was found most convenient in pre- 
y for each determination to grind the substance 
rT my tee mortar, moisten the powder slightly with 
ent grade xylene, and transfer a small amount of 
paste to the the ‘rmocouple junction with a micro- 
tula. All the substances that have been tested 
ere satisfactorily when treated in this way except 
issium sulfate, which decrepitates at the Ba 
sition point. 
lhe melting points were determined with the appa- 
is in a closed room but without the use of any 
ld against drafts except the two sheets of mica in 
r holder. Under these conditions the recorder 
wed random fluctuations up to +0.05 my, or about 
1,420° C. The results of the experiments, 
1, indicate a maximum error of +5degC 


leg al 
ed in table 

Velting point determinations with the hot-thermo- 
couple apparatus 


Melting point 
given in the lit 
erature [24] 


Melting point 


Compound found 


The value listed for barium disilicate is an average 
of three determinations--two made with different 
hermocouples on samples with no preliminary heat 
treatment, and the third after holding the sample for 
threat 1,420° C. All three estimates agreed within 
0.02 mv, or less than 2 deg. 


c. Other Applications 


The performance of the apparatus in these experi- 
ments indicated it to be satisfactory in the present 
form for determining melting points over the working 
range of the platinum and platinum-rhodium thermo- 
couple wire, with an accuracy of +5 deg even at 
1420° C. If greater precision is desired, the high- 
frequency power supply can be equipped with auto- 
matic regulation to decrease the amplitude of the 
random fluctuations in temperature. Any of the 
well-known methods of regulating heater voltage or 
current, or the resistance of the heating element [20], 
is obviously applicable; perhaps the most desirable 
regulator would be one sensitive to the thermocouple 
electromotive force, using the methods of converting 
the direet to alternating current, and of amplifying 
minute a-e voltages, which are well established in 
industrial controller designs. 

In considering the use of the hot-thermocouple 
apparatus for precise determination of melting 
point or other transition temperature, several factors 
must be taken into account. Volatility of the sub- 
stance renders the method useless if it causes the 

mple to disappear or to change significantly in 

mposition during an experiment. Sodium chlo- 
for instance, vaporizes rapidly enough to 


necessitate some haste in the melting-point deter- 
mination. Corrosion of the thermocouple by the 
molten sample also renders the method useless if the 
composition of the sample changes significantly. 
The effect of slight contamination on the electro- 
motive force of the thermocouple is probably less 
important as there is little temperature gradient 
within the contaminated portion of the wire near 
the junction. 

The thermocouple electromotive force is an ac- 
curate measure of the temperature at the center of 
the junction, but a comparatively large temperature 
gradient exists between the surface of the wire and 
the outer surface of the sample. This gradient is 
actually a help in growing single crystals, and causes 
no error in the determination of true melting points 
if the behavior of the material nearest the junction 
is observed. It may lead to segregation, however, 
if the solid and liquid have different compositions. 
For this reason the apparatus probably should be 
used in the investigation of multicomponent phase 
diagrams only for a rapid survey preliminary to 
detailed investigations by other methods. Bending 
the wire in such a way as to form a radiation shield 
[19] might well eliminate the temperature gradient 
within the sample, although probably at the cost of 
visibility. 

A temperature gradient is to be expected from the 
center of the thermocouple wire to its surface. The 
maximum temperature difference within the wire may 
be calculated approximately by the formula [25] 


AT=qr"/4k, 


where q is the rate of heat generation per unit volume 
of the wire, r is the radius of the wire, and & is its 
heat conductivity. The least favorable values of 
these quantities in the present apparatus may be 
taken as g=3.9< 10° cm (corresponding to 
the 8 w dissipated per centimeter of wire at 1,450° C), 
r=0.0125 em, and k=0.07 cal/see em deg (reported 
by Barratt [3] for platinum-—10-percent rhodium at 
room temperature). The temperature difference cal- 
culated from these values is 2.2 deg C. As this is 
less than the experimental error, the temperature 
gradient within the wire has been ignored in practice. 
The skin effect of alternating current tends to reduce 
the temperature gradient, but that effect is negligible 
at the frequency used in the present apparatus 

The principle used in this apparatus, heating a 
thermocouple by an alternating current that is fil- 
tered out so that the thermocouple clectromotive 
force can be measured alone, has been described in 
the literature [5, 14] as long ago as 1919. The idea 
still seems to offer possibilities in untried applications, 
such as the measurement of surface temperatures 
22], optical pyrometry [9], and elimination of errors 
caused by the heat conductivity of thermocouples 
{29}. 


cal sec 


3. Mounting Specimens 


When a single crystal has been grown in a droplet 
of melt, the preparation is cooled to room tempera- 
ture as rapidly as possible in order to minimize fur- 


155 





ther crystallization The bead of glass is readily 
broken off the wire by bending the latter away from 
the bead with fine pliers Then it is necessary to 
mount the specimen inh proper omentation on the 
usual glass fiber attached to a brass pin that fits the 
goniometer head [$2] of the X-ray apparatus 
Orienting of the erystal is facilitated by the visi 
bility of faces and edges during erystallization 
These often become more difficult to see after the 
but, if necessary, the interfering 
the surface of the glass sur- 


speciinen cools, 
effect of refraction at 
rounding the crystal can be overcome bry Hnmersion 
ina suitable index liquid 

Before the erystal is mounted, the glass fiber is 
attached to the goniometer pin. The adhesive best 
suited for this purpose and for mounting crystals 
stable at fairly high temperatures is jewelers’ shellac, 
which wets the surfaces readily when melted and 
hardens immediately on cooling. The fiber is ce- 
mented in place with the aid of a hand-held hot 
wire [26, p. 198] 

It is desirable that the glass fiber be capable of 
some tilting relative to the goniometer pin in case 
the accidental misorientation of the crystal is not 
fully correctible by means of the movable ares on the 
goniometer head, Frequently this has been done 
by softening the cement holding the fiber on the pin; 
another method has been the use of a goniometer pin 
with a thin, flexible section near the end [18, 27] 
The method used is to solder on the goniometer pin 
a I5-mm length of 0.4-mm (0.016-in.) copper wire, 
as shown in figure 3. The wire is amply rigid to 
support the fiber and specimen, vet it can easily be 
bent with fine pliers or tweezers when necessary 
The advantage of this arrangement is that tilting of 
the erystal can be done by displacing the lower end 
of the fiber, which is attached to the wire. Thus 
the ery stal is not decentered while being tilted, as it 
is when the goniometer pin itself bends 


GLASS FIBER 
. SHELLAC 


COPPER WIRE 


0.125-INCH 
BRASS ROD 





L 


Cioniometer pin with glass fiber 


Ficure 3 


When the fiber has been attached to the goniomet: 
pin, a small amount of shellac is melted onto t| 
upper end. The next step is then to hold this e 
of the fiber against the specimen in the proper orie) 
tation, heat until the shellac adheres, and allow | 
cool in the same position. The heating is done wit 
a tiny hot-wire tool made from a 5-mm length « 
0.015-mm platinum wire, which is connected to thy 
6.5-v transformer in series with a 22-ohm resistor 


4. Micromanipulator 


To facilitate the procedure of mounting the spec; 
men, some form of micromanipulator is desirab| 
Most of those described in the literature are either 
designed for use with a mono-objective microscope at 
high magnifications [2; 4; 10; 16, p. 62-72), and 
therefore delicate and limited in range, or adapted to 
coarser work but difficult to construct [11; 26, p 
202-4]. An attempt was therefore made to design 
an apparatus that would serve the purpose and yet 
require as little precise machining as possible in con 
struction. The apparatus was intended for use with 
a stereoscopic microscope on a vertical-pillar stand 
with horizontal extension arm and weighted base 
It will be deseribed in detail in the hope that the 
description may be useful to others and stimulate 
improvements in the kinematic design [26, 31] of a 
simple manipulator 

The complete apparatus is shown in figure 4. It 
consists of a base plate, A; two manipulators, B; and 
a three-legged platform, C. The base plate and the 
other flat pieces are cut from a ‘-inch steel plate, so 
that their weight insures stability 

Kach of the manipulators, B, is in principle simply 
an ordinary laboratory tripod stand whose base is a 
15° right triangle with leveling screws at the vertices 
The ends of these adjusting screws are positioned by 
six small blocks bearing 90° V-grooves, mounted on 
the base plate. In each set of three blocks, two are 
turned so that their grooves are parallel, and the 
third has its groove perpendicular to the direction 
of the other two. Each manipulator is thereby 
constrained to a single definite position at all times 

Platform C is an equilateral triangle. When set 
in place on the base slate its legs are positioned by 
the three larger V-blocks so that the center of the 
triangle is almost directly above one of the adjusting 
screws of each manipulator, Objects can be placed 
on the platform and operated upon with microtools 
attached to the jointed column that is mounted on 
each of the manipulators. 

The tools are usually moved by means of the ad 
justing screws at the sides and at the rear. Each 
screw causes its manipulator base to tilt about a line 
between the other two screws. The geometry of the 
arrangement is such that the tilts produced by the 
screws at the side and at the rear correspond to hor- 
izontal motions of the microtools suspended above 
the center of the triangular platform. As the serew 
below the center of the platform is at the 90° vertex 
of each manipulator base, the horizontal motions 
are at right angles. The screw at the rear moves 























Fieure 4. Mi 


A, Base plat 


romantp siator, 


B, manipulator; C, platform 


the tool back and forth, and the one at the side moves 
it from side to side, 

\s the height of the tool above the base plate Is 
ipproximately equal to each leg of the 45° right tri- 
angle, the motion of the tool is equal in sensitivity 
and maximum travel to the motion of the corre- 
sponding corner of the manipulator base along the 
adjusting serew. The present model has a sensitiv- 
ity of 0.OS em per turn and a maximum travel of 
25em 

\lthough the two horizontal motions are independ- 
ent, a vertical movement requires simultaneous 
turning of all three adjusting screws. Therefore, 
it is ordinarily most convenient to adjust the height 
beforehand and use only the horizontal motions if 
possible 

The column for supporting the tools, which is 
mounted on each manipulator, has two ball-and 
socket joints (tripod tilt tops, which may be ob- 
tained from most photographic supply houses). The 
two joints are connected by a length of l-in. round 
aluminum red, suitably threaded. At the upper end 
is a shorter length of the aluminum rod, drilled and 
equipped with setserews to hold the \-in. goniometer 
All the tools used with the micromanipulator 
are mounted on \-in. brass rods so that they can be 
attached in the same way. 

(he use of the ball-and-socket joints has the great 

antage that tightening them does not displace 


ibout 
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the object being adjusted, and the chance of acci- 
dental damage to delicate specimens is minimized 
As they are loosened during adjustment, there is no 
tendency to spring back out of position as with fric- 
tion joints. It is desirable that all four ball joints 
be capable of turning to any direction within a com- 
plete hemisphere. ‘Two of those in the present model 
are not, and the resulting limitation in flexibility is 
occasionally inconvenient, even though they are 
used in the lower positions. 

The adjusting serews, which are \ in. in diameter 
and have 32 threads per inch, are made with a con- 
ical recess in the lower end. This 90° cone rests on 
a steel ball 4 in. in diameter, which is positioned by 
the V-block. The combination provides a kinematic 
design in which the smooth, hard surface of the 
steel ball minimizes chattering when the screws are 
turned. Furthermore, if this surface becomes worn 
or damaged the steel ball can be replaced very 
cheaply. <A short piece of flexible plastic tubing 1s 
slipped over each V-block to prevent the steel balls 
from rolling off when one of the manipulators is 
lifted from the base plate 

The adjusting screws were cut on a lathe to insure 
that the screw threads would be coaxial with the 
conical surface; this was the only precise machining 
used in constructing the assembly. The corre- 
sponding threads in the manipulator bases were 
made by hand with a standard tap. 

The screw arrangement has disadvantages, being 
kinematically over-specified and therefore tending 
to allow lateral movement unless made to a very 
fine tolerance. This tendency was overcome to a 
considerable extent by preloading each screw with 
the simple device shown in figure 5, consisting of a 
leaf of Me-in. phosphor bronze with two tapped 
holes and a pointed setscrew. The tool can still be 
moved laterally by almost 0.001 in. if a horizontal 
force is deliberately exerted on one of the adjusting 
knobs, but the force required is great enough so 
that in practice no accidental motion occurs.  Pos- 


_— 














— 


Ficure 5. Spring loading of adjusting screw. 
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sibly other methods can be devised to eliminate 
the lateral movement entirely 

The manipulator has been used for almost a vear 
as & permanent adjunct to the stereoscopic micro 
scope. It has proved versatile and convenient for 
magnifications from “54 to «90. On 
sion, a standard petrographic microscope was set 
on the triangular table, and the manipulator was 
used to operate an implement on its stage at about 

150. Even under these conditions the precision 
and freedom from vibration were satisfactory, al 
though the increased height made the arrangement 
cumbersome 

The present model is less sensitive to vibration 
than the manipulator described by Strong [26, p 
202-4], probably because of the rigidity of the 1-in 
aluminum rod used for the jointed column, Fur 
thermore, the strength of the ball-and-socket joint 
when clamped makes it capable of supporting much 
heavier objects; even a small electric-motor grinder 
weighing more than 200 ¢ has been used as an 
implement 

It is believed that the simplicity, ruggedness, and 
cheapness of this apparatus make it valuable in 
many uses for which more elaborate equipment is 
not absolutely in addition to the mount 
ing of crystals 


one occu 


hecessary, 


5. Grinding 


The specimen containing a desired crystal may 
also contain extra erystals or excessive amounts of 
glass. If so, it is necessary to remove the unde 
sirable portions of the specimen and leave the ervs 
tal undamaged, This can best be done by mount 
ing the specimen as usual, in such a position that 
the undesirable portion is exposed, and then care 
fully removing that portion. A similar procedure 
of mounting and then removing the excess is nee 
essary if cylindrical specimens are desired to sim 
plifv absorption corrections 

The slow removal of a part of the specimen is 
often done by hand with a solvent applied on a 
camel’s-hair brush. In some cases no suitable sol 
vent can be found, and mechanical removal becomes 
necessary. A small grinder to be used for this pur 
pose on the micromanipulator was made from a 
miniature 27.5-v Alnico field motor, by fitting a 
simple adapter on the shaft so that standard abra 
sive tips with \-in. shanks could be used. The 
motor operates satisfactorily on 120-v 60-cyele alter- 
nating current, with a 25-watt bulb and selenium 
rectifier in and a 1,000-microfarad electro 
lytic condenser in parallel. The best grinding sur 
face was obtained by stirring a pinch of 5-4 diamond 
powder into a few drops of thinned lacquer and 
coating a \-in. polished steel shaft with the mixture 
A surface thus prepared was found to cut many 
times faster than any other that has been tried 
In addition, its relative smoothness caused much 
leas vibration of the specimen than had previously 
been considered unavoidable, reducing considerably 


Series 


the chance of accidental breakage or loss. The 
of diamond powder in other types of erystal-grind 
apparatus, such as that of Sturdivant [27] for au 
matically producing a cylinder, would probably ha 
similar advantages; its added convenience mal 
the effort of procurement worthwhile, and its « 
per specimen ts slight 
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Crystal Forms of Chromium Orthophosphate 
Barbara M. Sullivan and Howard F. McMurdie 


There is no evidence for the existence of any crystalline hydrate of chromium ortho- 


phosphate other than the hexahydrate 
material 


to an unstable, fine-grained compound (8CrPO* 
aCrPO*') at 972° C 
forms, 


strained, it will form a stable compound 
patterns are given for the three crystalline 


If the unstrained and coarse hexahydrate is heated, it recrystallizes at 


This breaks down near 130°C (to an amorphous 


Lon © 
If the original hydrate is ground and thus 
X-ray diffraction powder 


and the general relation of chromium 


orthophosphate hexahydrate to other hydrates is discussed 


l. Introduction 


\ study of the forms of chromium orthophosphate 
suggested by the National Institutes of Health 
y needed a radioactive material that was insolu- 
nd nontoxic. Chromium orthophosphate with 
ulioactive phosphorus had been used in similar work 
They found, however, that it was difficult: to 
prepare a uniform anhydrous compound by heating 
the hexahydrate, inasmuch as the temperature of 
dehydration was very critical 
There are only a few previous studies of the com- 
pound reported in the literature. Joseph and Rae 
2}in 1917 reported a violet hexahydrate, as well as a 
rystalline tetra- and dihydrate. In 1948 Clark and 
lai [3] made X-ray powder patterns of the hexa- and 
tetrahydrates and found them identical, Whereas pat 
terns of the dihydrate and anhydrous forms indicated 
amorphous material, 


2. Materials and Equipment 


lhe samples studied here were given to the Bureau 
by the National Institutes of Health. They were 
prepared as hexahydrate (CrPO,.6H,O), using a mod- 
ification of the method of Joseph and Rae [2] de- 
scribed in detail by Ness, Smith, and Evans [4]. 
Differential thermal-analysis studies were made by 
standard methods [5], using a program controller, a 
direct-current amplifier, and automatic recording. 
X-ray patterns were made on a North American 
Philips High-Angle Geiger Counter Spectrometer, 
sing the techniques described by Swanson and 
latge [6] to obtain accurate spacings and intensities. 
ln some cases patterns were rade with the sample at 
high temperatures in the furnace described by Van 
Valkenburg and MecMurdie [7]. CuKa radiation 
was used in all cases (A= 1.5405) 


3. Chromium Orthophosphate Hexahydrate 


lhe hexahydrate, as shown in figure 1, is violet in 

color, prismatic and twinned in habit, biaxial nega- 

with a=1.568, B=1.591, y=1.599, 2V=13° 

d spacings and relative intensities of the powder 
tern are shown in table 1. 


s in brackets indicate the literature references at the end of this paper 


Figure 1. Chromium orthophosphate hexahydrate; 150 


Panter 1. d-Spacings and relative intensities of chrominm 


orthophosphate herahydrate 


The chromium orthophosphate hexahydrate, when 
heated, begins to lose water at about 100° C and 
loses approximately 95 percent of the total between 
100° and 400° C. It is not completely dehydrated 

’. The details of the losses have 


until about 800° C, 
been discussed by Ness, Smith, and Evans [4]. 
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They showed that while most of the loss is near 130° 
C, the vapor-pressure curve indicates no sharp 
breaks, and there is considerable loss, even up to 
800° C. With the initial loss of water, the material 
changes first to a bright-green color and then becomes 
dark brown. The X-ray powder pattern of material 








Differential thermal ana »} CrP 6H,O 
heated at 12 2 


aCrPO,t 


Figure 2 


a 


Ficure 3. 
orthophosphate 


Amor phous material obtained by heating chromium 
hezahydrate at 600° C for 66% hours; 150. 


heated up to 130° C becomes gradually weaker } 
is otherwise unchanged. Above 130° C the patt: 
indicates amorphous material. No evidence w; 
found here for the existence of the crystalline tetr; 
hydrate reported by others [2, 3]. It was report, 
as being produced by boiling the hexhydrate in H,0 
for half an hour. Material treated in this mann 
was bright green and showed a weak X-ray patter: 
of the hexahydrate. Under the microscope it ap 
peared as a mixture of the crystalline hexahydrat; 
and amorphous material. Weight-loss determina 
tions of the boiled material indicated the loss of 
approximately 1 molecule of water. 

The differential thermal analysis patterns of 
chromium orthophosphate hexahydrate, both ground 
and unground, heated at 12% deg/min, are given in 
figure 2. The large endothermic break near 130° C 
is caused by the loss of water and the breakdown of the 
crystal structure. Above that temperature the mate- 
rial appears dark brown and almost opaque in the 
microscope,and the X-ray diffraction pattern indicates 
that it is amorphous. Figure 3 shows a sample that 
has been held at 600° C for 66% hours. It may be 
noted that the external shape of the original hexahy- 
drate is retained. This is typical of material heated 
in this temperature range. 


4. Anhydrous Chromium Orthophosphate 


The exothermic breaks on the differential thermal- 
analysis curves (fig. 2) near 1,000° C indicate the 
crystallization of anhydrous chromium orthophos- 
phate. Two different crystal forms exist, one of 
which is stable. The stable high-temperature form 
is designated a, and the lower, unstable form, 8 
The 8 form is produced when unground material 
is heated and its formation produces the sharp break 
at 1,014° C. At higher temperatures this form grad- 
ually changes to a chromium orthophosphate. Al- 
though the X-ray pattern does not show the «@ form 
until it has been heated to 1,400° C, microscopic 
examination indicates its formation in small amounts 
as low as 1,200° C. Fairly pure 8 chromium ortho- 
phosphate can be formed when the unground hexa- 
hydrate is heated a short time between 1,000° and 
1,100° C, but prolonged heating near 1,100° C results 
in the formation of a noticeable amount of the a form 
When the hexahydrate is ground before heating, the « 
form crystallizes with no intermediate step. At a 
heating rate of 12% deg/min this takes place with 
an exothermic break at 972° C. In the high-tem- 
perature X-ray furnace [7] only the a form was 
found, as it was necessary to grind the sample to 
hold it in position. 

Tables 2 and 3 list the d spacings and relative 
intensities of the two anhydrous forms. <A typical 
example of the 8 form is shown in figure 4. The 
very fine highly birefringent particles show no extinc- 
tion positions and have an average index of 1.908 
The a@ form (fig. 5) is strongly pleochroic, changing 
from bright green to deep blue. It is biaxial positive 
with a very large 2), and with a=1.761 (blue), and 
y=1.844 (green). The presence of small amounts 
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TaBLe 2. d-Spacings and relative intensities of anhydrous 
chromium orthophosphate—8 form 


d-Spacings and relative intensities of anhydrous 
chromium orthophosphate a form 


5. Discussion 

Anhydrous chromiun orthophosphate 8S torm) 

heating the unground heaahudrate to 1.100° C 
ential thermal-analys furnace; 1450. 


No evidence was found in this investigation for 
the existence of any crystalline hydrate of chromium 
Er, orthophosphate other than the hexahydrate. The 
dehydration is not completed below 800° C because 
of the formation of an amorphous material that has 
an extremely large absorbing surface. The dehydra- 


tion curves of Ness, Smith, and Evans [4] do not 
show sharp breaks as do compounds that form new 
crystal phases immediately following dehydration, 
but are much more gradual, which, it is believed, is 
typical of conditions when amorphous phases exist 
for scme time before recrystallization. 

The formation of an amorphous phase that exists 
over a wide temperature range is of special interest 
here. It is known that kaolinite forms such a phase 
between its decomposition at about 550° C and the 
formation of yAl,O; at 980° C [8], and that certain 
other silicates, such as pyrophyllite [9], pass through 
an amorphous state. These compounds differ in 
two distinct ways from chromium orthophosphate 
hexahydrate. In these silicates the ‘“‘water’’ is in the 
form of OH groups and is lost above 500° C, and the 
silicate groups are linked together with shared O 
atoms. In the case of CrPO,. 6H,O, however, it is 
reasonable to expect that the water is present as 
molecular water surrounding the Cr atoms, and the 
PO, groups are discrete unlinked groups. 

Aluminum and iron orthophosphates have inver- 

‘ sions and structures similar to those of silica [10]. 
URES. Anhydrous chromium orthophosphate (a form) | Tt was at first thought that chromium orthophosphate 
tained by heating the hexahydrate at 1,500° C for 2 hours; oye 
140, ' , . was related to silica in an analogous manner, but no 
evidence was found to support this. 
of the a form is readily detected by the microscope The formation of an unstable form of the anhy- 

‘cause of its unusual pleochroism. drous compound, when unground material is used, 

The a form was stable and unmelted up to 1,800° C, | is evidence that the 8 form is more closely allied to 
the limit of available furnaces the original structure than is the a form. When 





the starting material is ground with an attendant 
increase in surface, it approaches equilibrium more 
rapidly and the 8 form is not obtained. It is ap- 
parent that the 8 form is not stable above 973° C. 

The formation of a preliminary unstable phase is 
not an unusual occurrence; the formation of yAl,O, 
in kaolinite is such a case [8]. When 3CaQ.Al,O,.6H,O 
is dehydrated, the first crystal phases formed are 
12C'aO.7Al,O, and CaO, whereas the final equilibrium 
compound is the single phase 3CaQ.Al,O, [11| 
CGirim and Bradley [9] have noted the occurrence of 
various unstable phases in dehydration of clays and 
other silicates 


6. Summary 


There is a single crystalline hydrate of chromium 
orthophosphate (CrPO,6H,O) that is violet, platy, 
biaxial negative with a= 1.568, B=1.591, y—1.599, 
and 2V=13°. 

This hydrate breaks down near 130° C to an 
amorphous material that retains some water up to 
so0e © 

If the original material is unstrained and coarse, 
it will reerystallize on heating at about 1,020° C 
to an unstable fine-grained BCrPO, with an average 
index of refraction of 1.908. 

If the original hydrate is ground and thus strained, 
the material will form a stable compound (aCrPO,) 
that is highly pleochroic, biaxial negative with 
a—1.761, y~ 1.844, and 2V large. 


O 


7. References 

{1} H. Allen, L. H. Hempelmann, Jr., and N. A. Womae! 

The effect of insoluble radiophosphorus (chromiy 

phosphate) when applied interstitially in the trea 

ment of adenocarcinoma of the mamma in mic 

Cancer Research §, 239 (1945). 

[2] A. F. Joseph and W. N. Rae, Chromium phosphat: 
J. Chem. Soc. UM, 196 (1917). 

{3} G. L. Clark and A. P. Tai, Phosphate coating of alum 
num and polymorphism of chromium phosphat 
Science 107, 505 (1048). 

[4] A. T. Ness, R. E. Smith, and R. L. Evans, The prey 

aration and properties of normal and radix 

active chromic phosphate, J. Am. Chem. Soc. (publ 
cation pending 

Speil, L. H. Berkelhamer, J. Pask, 

Differential thermal analysis, U. 8. 

Tech. Pap. 664 (1945). 

[6] H. BE. Swanson and EF. Tatge, Standard X-ray diffractio: 

patterns, J. Research NBS 46, 318 (1951) RP2202 

Van Valkenburg, Jr., and H. F. MeMurdie, High 

temperature X-ray diffraction apparatus, J. Researc! 

NBS 38, 415 (1947) RP1782. 

H. Insley and R. H. Ewell, Thermal behavior of kaolin 
minerals, J. Research NBS 14, 615 (1935) RP792. 

[9] W. F. Bradley and R. E. Grim, High-temperature 
thermal effects of clay and related materials, Am 
Mineral, 36, 182 (1951 

M. J. Buerger, Crystals based on the silica structure, 
Am. Mineral. 33, 751 (1948). 

{11} W. G. Sehneider and T. Thorvaldson, The dehydration 

of tricalcium aluminate hexahydrate, Can. J. Research 
19 (B), 123 (1941 


some 


and B. Davis 


[5] 8. 
Bureau of Ming 


[7] A 


[8] 


{10} 


WasuHiIncton, October 10, 1951. 





